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The Chairman’s Column 


To MEMBERS OF THE PGEC 


This year I hope to keep all of you informed on cur- 
rent and future happenings in the PGEC. I hope to do 
it principally through this column, which will appear 
(when there is something to be said) in each issue of 
these TRANSACTIONS. 

On April 29, C. W. Rosenthal of the Administrative 
Committee attended, in my behalf, a meeting in New 
York of the national Professional Groups Committee. 
The following summary is based on his report to me. 

1) Joint PG chapters will be approved where neces- 

sary and advisable. 

2) PG members who are 4 months delinquent in dues 
will be dropped from membership. 

3) There is activity under way to standardize the 
covers of the various Group Transactions and 
also the forms used in reviewing submitted papers. 

4) The financial status of national Professional 
Groups was discussed. An Ad Hoc Committee 
has been formed to study this problem and mem- 
ber’s suggestions to cure it (e.g., to accept adver- 
tising in the TRANSACTIONS). Dr. R. M. Emberson, 
rR. Lo MeFarian, Dr, L.. GC; Van Atta and: the 
undersigned have been appointed to this com- 
mittee. The chairman will be Dr. John T. Hender- 
son, past president of the IRE and chairman of 
the board of directors. 

5) A survey of the financial affairs of Professional 
Group chapters indicated that the principal prob- 
lems were not financial, but rather in personnel 
and organizational matters. The personnel prob- 
lem is one of strong leadership and an adequate 
carry-through from year to year, while organiza- 
tional difficulties are best helped by the Pro- 
fessional Groups Chapter Coordinator of the cog- 
nizant Section. 

6) There is some feeling among Professional Groups 
which do not hold any national meetings of their 
own, that they should get preferred treatment and 
time at the national IRE Convention, and that 
Groups which hold national meetings should re- 
strict their national convention papers to tutorial 
ones or ones of broad interest; specialized papers 
should be presented at group national meetings. 

During the 1958 Western Joint Computer Conference 
the Administrative Committee held an informal meet- 
ing. Since there had been a meeting about six weeks 
previously in New York there was no pressing business 
to transact and consequently, no formal actions were 
taken. 

However, a number of items of business were dis- 
cussed; significant findings are outlined here. 

1) With respect to the problem of obtaining an ab- 

stract service to replace one previously conducted 


by Dr. Harry Huskey on a volunteer basis, the 
consensus of those present was that the PGEC 
should spearhead this effort and not wait for the 
cooperation of other technical societies. Frank 
Heart of Lincoln Laboratories and Larry Jones 
of Westinghouse-Baltimore, the Ad Hoc Committee 
for this problem, are preparing a specific proposal 
for consideration by the Administrative Commit- 
tee meeting at its August 21 meeting during 
WESCON in Los Angeles. 

2) Members present felt that we have no need at the 
moment of financial return from advertising in 
these TRANSACTIONS; for the moment, therefore, 
we will not accept advertising. 

3) The National Joint Computer Conference re- 
cently amended its proposed charter to meet the 
objections of the PGEC Administrative Commit- 
tee. The matter will be voted upon by the Ad- 
ministrative Committee shortly. 

4) Those present were in favor of rerunning the 
membership survey which was conducted in the 
fall of 1956 since it was felt that significant changes 
have occurred in the last 18 months. As before, 
to protect the anonymity of the information, a 
large public accounting firm will be used as the 
mail-drop, and will carefully dissociate any identi- 
fying information. Since this matter is subject 
to the vote of the Administrative Committee, no 
action has been taken yet. 

5) The problem of PGEC participation in an Inter- 
national Federation of Computer Societies was 
discussed. This is a matter which requires JRE 
attention, and Dr. Arnold Cohen is investigating. 

6) As of May 1, the PGEC had 7338 members. 

7) Chairman Melville of the Student Activities Com- 
mittee indicated that he is not in favor of con- 
tinuing our postdoctoral fellowship beyond this 
year, primarily because fellowship money is 
readily available elsewhere, and the PGEC does 
not receive a significant amount of publicity. An 
alternate plan will be submitted by Melville for 
Administrative Committee consideration. 

8) Editor Nash indicated that he will recommend 
some changes in the membership and structure of 
the editorial board and Publications Committee. 

Plans are now firm for an International Conference 

on Information Processing to be held June 13-21, 1959. 
Sponsored by UNESCO, the conference will be in Paris. 
Only invited papers will be used by the U. S. Committee 
for the ICIP. Further information may be obtained from 
the committee at Box 4999, Washington 8, D. C. All 
travel arrangements must be coordinated through this 


committee. 
WILuis H. WARE 
Chairman 
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Willis H. Ware (A’43—SM’49) was born in At- 
lantic City, N. J., on August 31, 1920. He received 
the B.S. degree from the University of Pennsyl- 
vania in 1941, and the S.M. degree in 1942 from 
Massachusetts Institute of Technology, where he 
was a Tau Beta Pi Fellow. 

From 1942 to 1946 he was employed by the 
Hazeltine Electronics Corporation for research 
and development in radar and IFF. In 1946 he be- 
came one of the original members of the staff of 
the Electronic Computer Project at the Institute 
for Advanced Study, Princeton, N. J. There he 
worked on the design and development of the 
large scale general purpose electronic digital com- 
puter which was later to set the pattern for the 
construction of several other “Princeton class ma- 
chines.” At the same time he continued his gradu- 
ate work at Princeton University, from which he 
received the Ph.D. degree in 1951. In 1952, Dr. 
Ware joined the RAND Corporation, where he is 
continuing his work with digital machines. 

Dr. Ware has held several national and local 
offices in the PGEC. He is National Chairman for 
1958-1959 and has held posts as National Vice- 


Willis H. Ware 


Chairman, 1957-1958; National Secretary-Treas- 
urer, 1954-1955; National Chairman of the Stu- 
dent Activities Committee, 1956-1957; Chairman 
of the Los Angeles Chapter, 1954-1955; Treasurer 
of the Los Angeles Chapter, 1953-1954. He also 
has been a member of IRE Technical Committees 
8.4 and 8.5 (computer definitions) from 1948 to 
the present time, and served as Chairman of Com- 
mittee 8.5 from 1951 to 1958. 

Conference Chairman of the 1958 Western 
Joint Computer Conference, Dr. Ware is cur- 


- rently serving as a member of the IRE national Ad 


Hoc Committee to review Professional Group Fi- 
nances. He is also a member of the Technical Pro- 
gram Committee of the International Conference 
on Information Processing. In 1957 he was Fi- 
nance Chairman of the Western Joint Computer 
Conference and also the Professional Groups Co- 
ordinator of the Los Angeles section. In 1956 he 
was the alternate Chairman of the Technical Pro- 
gram Committee of WESCON. 

Dr. Ware is a member of Tau Beta Pi, Sigma Xi, 
Eta Kappa Nu, Pi Mu Epsilon, and the American 
Association for the Advancement of Science. 
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Design of AC Computing Amplifiers 
Using Transistors* 
C. A. KRAUSE} anp R. R. LOWEt 


Summary—A design philosophy for transistorized analog com- 
puting amplific-s is presented. A design procedure for a summing 
amplifier to drive a specific resolver in a 400 cps system is described, 
and the performance of the resulting circuit is evaluated. Whereas 
experience with vacuum tube amplifiers in similar applications has 
led to the conclusion that the amplifier input impedance should be as 
large as possible, the inverse is true in the transistor amplifier. 


I. INTRODUCTION 


T HAS been felt that the transistor had inherent dis- 

| advantages in analog summing amplifiers because 
of its relatively low input impedance. Investigation 
shows that this is not the case and that, in general, a 
current amplifying device is most valuable when its 
input impedance approaches zero.! This and subsequent 
conclusions resulted from the development of a silicon 
transistor summing amplifier to drive a resolver of 
15 kilohms tuned resistance at a computing frequency 
of 400 cps. It was desired that the summing resistors 
be as large as 500 kilohms (this will be seen to have a 
direct bearing upon the required current gain), and that 
the computing error due to finite amplifier gain not 
exceed 0.05 per cent when the closed loop gain is unity. 
The various configurations in which the amplifier- 
resolver combination are to be used are shown in Fig. 1. 
Parallel and series compensator feedback provide multi- 
plication of the input by the sine and/or cosine of the 
resolver angle, series feedback having potentially higher 
loop gain but restricted to a single input. Parallel rotor 
feedback provides for division by the same trigonomet- 
ric functions. Since the rotor feedback case can never 
have any more loop gain than the compensator feed- 
‘back cases, it is the less critical application (ignoring 
resolver idiosyncrasies). The approach taken was to 
provide adequate loop gain in the parallel compensator 
feedback case, so that the series feedback configuration 
had reduced significance. Thus provision for adequate 
loop gain and stability margins in the case of parallel 
compensator feedback is essentially a complete solution. 


Il. SUMMING AMPLIFIER INPUT IMPEDANCE 


The impedance at the summing node of a good com- 
puting amplifier is very nearly zero. The input imped- 
ance seen from a summing input to the amplifier is, 
then, essentially the value of the summing resistor, 


* Manuscript received by the PGEC, December 18, 1956; revised 
manuscript received, April 18, 1958. ; 

+ Norden Division, United Aircraft Corp., Gardena, Calif. 

1R, L. Wallace, Jr., and G. Reisbeck, “Duality as a guide in tran- 
sistor circuit design,” Bell Sys. Tech. J., vol. 30, pp. 381-418; April, 
1951. 
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Fig. 1—Various configurations of the amplifier-resolver 
combination. 


regardless of the amplifier’s own internal input imped- 
ance. In a vacuum tube summing amplifier, the inter- 
nal input impedance should be as high as possible to 
prevent deterioration of loop voltage gain by its shunt- 
ing effect. In a transistor summing amplifier, the internal 
input impedance should be as small as possible so that a 
voltage error in the relationships at the extremities of 
the feedback and summing resistors will produce the 
maximum possible current signal into the transistor. 
Consider the amplifiers of Fig. 2. In each circuit, the 
finite gain is evidenced by the fact that e;/R,; is not 
equal to éo/Ry, 7.e., a non-zero current flows into the 
summing point. Assuming that the voltage at the 
amplifier input is small in each case, the error term can 
be approximated as 


Ry 
tetera 


Expressing & in terms of each of the two circuits, 


eg 1 
222 eee 
vt Ry Ry 
1+—4+— 
Roe 
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Fig. 2—Approximate equivalent circuits for vacuum tube and tran- 
sistor summing amplifiers. (a) Vacuum tube amplifier. (b) Tran- 
sistor amplifier. 


and 
Cb il 
—- ) 
Ry Ry 
1+—+— 
ae R, 


V5 


In the vacuum tube case, e,/E,, should clearly be 
maximized, leading to the conclusion that R, should 
also be maximized. Because the two expressions are 
similar, the same conclusion will be reached regarding 
the R, of the transistor unless care is exercised. As the 
transistor amplifies 7; and not ¢, and e,=7;7, it can be 
shown that 
1; 1 


Ew a ; 
11 1 + R R; 


Thus 7; should be minimized in order to maximize 
4i/ Ftp. 

Since this conclusion is directly opposed to the 
vacuum tube case and to some earlier opinions regarding 
transistor amplifiers, it must be justified. Consider the 
approximate equivalent circuit of a summing amplifier 
shown in Fig. 3. Assuming that e is very small for the 
purpose of computing the power in Ry, and defining 
the power gain, Ap, as 


€0 (A 
A» = ed ) 
R,/ UG 
where 
P RR; 
L= , 
Rr + Ry 
it can be shown that 
€0 Ry 1 


Ay Ry /1; + Ry/Rs + 1 
V ApR1/ 1; 
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Fig. 3—Approximate equivalent circuit of a transistor summing 
amplifier. (Note: For small-signal, low-frequency application, 
any degree of rigor may be obtained from this circuit by properly 
defining 7; and A;.) 
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Fig. 4—Required and developed power gain as a function of input 
resistance for vacuum tube and transistor summing amplifiers. 


If the desired ratio between e) and e, is —R;/R,, the 
magnitude of the error term, £, is 


Mp (Ee :) 
AR rt RS ; 


Fig. 4 shows A, as a function of 7; for various values of 
R,/R;, (the desired closed loop gain) with | E| =0.05 
per cent. 

Referring to the A,-7r; plot, the dashed lines labeled 
TR and VT show the developed power gain input im- 
pedance for a transistor amplifier and a vacuum tube 
amplifier, respectively, both of which were designed 
for the “matched” condition with 25 computing in- 
puts. For the vacuum tube case, an otherwise identi- 
cal amplifier develops more power gain as the input 
impedance (grid resistor) is increased, even providing 
slightly more gain than required for the specified ac- 
curacy. This is intuitively satisfying and indicates de- 
creasing power requirement in the input circuit (at 
constant output) as the grid resistor is increased. For 
the transistor case, an otherwise identical amplifier 
develops more power gain as the input impedance 
(characteristic of the first stage, or intentionally added 
series resistor) is decreased, again providing slightly 
more gain than is required for the specified accuracy, 
and decreased input power requirement at constant 
output. 

Another inference from the A,-R: plot is that tran- 
sistor amplifier computing accuracy is relatively in- 
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Fig. 5—Accuracy as a function of the number of summing resistors 
for amplifier current gain of 100,000 and different values of input 
resistance. 


dependent of the number of inputs. To a first approxi- 
mation, the number of inputs (at equal summing resist- 
ance) required to reduce the amplifier accuracy to one- 
half of the single-input value is equal to the ratio of the 
summing resistance to the amplifier input resistance. 
For example, an amplifier having 0.05 per cent accuracy 
with a 500K feedback resistor and one 500K summing 
resistor, and with an input resistance of 10K, will have 
an accuracy of about 0.1 per cent with fifty 500K in- 
puts. Alternatively, multiplication of a single input by 
variation of an input summing resistor can be carried 
out with considerably more accuracy than in the com- 
parable vacuum-tube circumstance. The amplifier hav- 
ing characteristics described above can maintain half 
the single-input accuracy in multiplications up to 50:1. 
This point is shown graphically in Fig. 5, which shows 
the accuracy of such an amplifier having a current gain 
of 100,000 for various values of Rin. The curve plotted 
for an Rj, of 7.5 kilohms is of special significance, since 
this corresponds approximately to the final design. 


III. Loop Gain 


It has been shown that multiple summing resistors 
have little effect upon a properly designed transistor 
amplifier. It is also true that the magnitude of the load 
resistor has little effect upon the vacuum tube amplifier, 
but a large effect upon the transistor version. As 7; ap- 


proaches zero, 
€0 a Ry 
lee oe Rp 2 Rs 
Ry 


Thus a larger value of current gain is required if Rz is 
reduced. For instance, if the tuned resistance of the 
resolver is reduced from 15 kilohms to 10 kilohms, the 
gain requirement is increased almost 50 per cent. The 
percentage of the amplifier current which is fed back 
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Fig. 6—Final amplifier configuration. 
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Fig. 7—Equations for collector current stabilization. B is the dc 
current gain at the operating point and 7, re, 7c, and J,. are used 
in the conventional sense. 


is approximately inversely proportional to the ratio 
of the feedback resistor to the load resistor. 


IV. AMPLIFIER CONFIGURATION 
A. Operating Point Stabilization 


The methods of operating point stabilization should 
be established before the basic amplifier configuration is 
finalized, since such stabilization can be achieved only 
at the expense of dynamic current gain. Thus the degree 
of required stabilization affects the choice of transistor 
types and the manner in which these transistors must 
be connected and interconnected. The first two stages 
of Fig. 6 have been stabilized by means of the collector- 
to-base feedback resistors, the base-to-ground shunt re- 
sistor, and the emitter resistors. 

Equations can be developed? which express the degree 
of operating point stability as a function of transistor 
parameters and external resistors (see Fig. 7). These 
equations give quite realistic results in the case of 
germanium transistors, but the current-amplifying prop- 
erties peculiar to silicon reduce their usefulness. Since 


2A. W. Lo, R. O. Endres, e¢ al., “Transistor Electronics,” Pren- 
tice-Hall, Inc., New York, N. Y., pp. 131-154; 1955. 
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the current gain of the grown junction silicon unit is 
about zero at zero emitter current, increases quite 
rapidly to some maximum and then decreases more slow- 
ly as emitter current increases, the small-signal value 
of current gain gives only qualitative results. These 
equations were used only to establish the order of mag- 
nitude of the stabilizing components, and extensive 
laboratory measurements were employed to optimize 
their values. 

It would have been possible to improve the operating 
point stability of the amplifier with the use of fewer 
components if the second stage had been direct-coupled 
to the first in the same manner as the third and fourth 
stages are each coupled to the preceding stage. As will 
be shown in the discussion of low-frequency stabiliza- 
tion, this can be accomplished only at the expense of an 
undesirably large by-pass capacitor. In the existing 
amplifier, the third and fourth stages have extremely 
stable operating points, established almost entirely by 
the base input voltages derived from the previous 
stage, and relatively large emitter resistors. 

The proper use of temperature-sensitive resistors in 
the compensating networks can improve operating point 
stability. This technique was not necessary in this 
amplifier, since the operating currents are quite narrow- 
ly confined in the regions of maximum current gain by 
the existing networks. No change of maximum un- 
distorted output with temperature was observed, and 
Fig. 8 shows that the gain decrease with temperature 
is no more than expected due to current gain decrease 
with temperature, 7.e., current level shift has not con- 
tributed significantly to the effect. As will be discussed 
later, temperature sensitive resistors might be properly 
incorporated to partially compensate for this gain de- 
crease. The fact that the gain of the amplifier did not 
significantly increase at high temperatures where the 
individual transistor current gains have increased 20 
per cent or more is due to increased insertion losses 
caused by higher input and lower output impedances of 
the individual stages. 


B. Individual Stages 


The method of connecting the individual stages is 
largely determined by the desire for maximum current 
gain from a given number of transistors without the use 
of transformers. The common base connection is im- 
mediately eliminated because it has a current gain of 
less than unity. The common collector connection is 
generally less efficient than the common emitter because 
of the large insertion loss caused by its relatively large 
input impedance and low output impedance. However, 
where four stages must produce a net reversal of phase, 
the use of a single non-phase-reversing common col- 
lector makes it unnecessary to obtain a phase reversal in 
the resolver load, thereby minimizing the adverse 
affects of interwinding capacities. Since the input im- 
pedance of a common collector stage is minimized when 
its load resistance is smallest, and since the input resist- 
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Fig. 8—Relative open loop current gain as a 
function of temperature. 


ance of the Type 953 transistor is significantly less than 
the Type 905, making the third stage the single common 
collector stage minimizes the insertion loss. 


C. Transistor Types 


To the knowledge of the authors, the single manufac- 
turer of production quantities of high quality silicon 
junction transistors at the time this amplifier was de- 
signed was Texas Instruments, Incorporated. Since the 
required output voltage was to be a minimum of 20 
volts rms, only the Types 952 (80 v) and 953 (120 v) 
have adequate collector voltage maxima. The 953 was 
chosen because collector voltage maxima are reduced 
at high temperatures, because an appreciable voltage 
must be maintained across the transistor for class A 
operation, because the maximum undistorted output 
was desirable, and because the price and availability 
of the two types were identical. 

Since the minimum guaranteed common-emitter cur- 
rent gain of the Type 953 is less than 10, and since the 
required total current gain must be about 100,000, it is 
desirable that the other three stages have a large guaran- 
teed minimum current gain. This is also desirable from 


‘the standpoints of allowing maximum ac and de de- 


generation and optimizing low temperature perform- 
ance. Choice of the Type 905 with a minimum common 
emitter current gain of about 40 (this is the best avail- 
able) gives performance, under the most adverse con- 
ditions, that is marginally adequate. The current gain 
thus achieved and its effect upon computing accuracy 
will be discussed in Section VI. 


V. DYNAMIC STABILIZATION 


The dynamic stabilization problem can logically be 
divided into two parts—e.g., instabilities above and 
below the nominal 400 cps operating frequency. The 
high-frequency problem is of little concern in this dis- 
cussion since its origin and its solution both stem from 
the specific resolver load. Because the particular resolver 
was tuned and because its dissipation is low, the gain 
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ecreases rapidly with high frequencies. It was possible 
to place the 180° phase shift frequency well beyond the 
nity gain frequency by eliminating the effect of the 
resolver leakage reactance. This can be accomplished 
py a single capacitor in the case of compensating winding 
feedback, and by a form of parallel tee network for the 
other configurations. The gain-phase characteristics of 
the transistors had little to do with the problem because 
f the sharply tuned load. 
_ Low-frequency oscillations are primarily caused by 
leading phase shifts due to coupling and by-pass capaci- 
tors. The philosophy employed in the solution of this 
problem was the attempt to create a single “dominant 
lead,” wherein one resistance-capacitance time con- 
stant has reduced the open loop gain below unity before 
other time constants make significant contributions to 
tthe phase shift. The dominant lead of the circuit of 
Fig. 7 results from the 40 microfarad emitter by-pass 
capacitor of the fourth stage, the input resistance and 
‘current gain of this stage, and the output resistance of 
‘the third stage. This time constant, about 0.001 second, 
lis at least an order of magnitude smaller than all others. 
If the second stage had been direct-coupled to the 
first, the by-pass capacitor for the emitter resistor 
‘would have been much greater than the 40 microfarad 
‘coupling capacitor, as impedances from emitter to 
ground are multiplied by the stage current gain in their 
effect upon input impedance. If the large time constant 
‘created by the 40 microfarad coupling capacitor were to 
be maintained, the by-pass capacitor required for direct 
‘coupling would have to be about 2000 microfarads, 
‘regardless of the size of the emitter resistor. 


| VI. DrEstGN MopIFICATIONS FOR IMPROVED 
ACCURACY 


_ The computing accuracy of the described type of 
summing amplifier can be improved by several meth- 
‘ods. Such methods as reducing the emitter degeneration 
and/or reducing the size of the summing resistors have 
the disadvantage of reducing the gain margins. Emitter 
resistors with a positive temperature coefficient might 
be used to give the same high-temperature margins 
with reduced low-temperature degeneration. The com- 
puting accuracy might also be improved by obtaining 
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the transistors to a more limited current gain tolerance, 
but there are inherent disadvantages of selection tech- 
niques. 

It may also be practical in some cases to employ 
temperature control such that the lowered current 
gains at very low temperature are eliminated. 

In this particular case, it was originally assumed that 
the ambient temperature could be maintained at or 
above 25°C, if necessary. The resultant amplifier under 
the most adverse conditions maintained the desired 
accuracy to about —30°C. Reduction of the 100 ohm 
emitter resistors to 47 ohms, or of the 500 kilohm sum- 
ming resistors to 300 kilohms maintains the 0.05 per 
cent accuracy to —55°C, but lowers the gain margins by 
about 4 db. A minimum “alpha” specification for the 
type 953 transistor of about 0.94, rather than 0.90, would 
allow the amplifier to operate within the specified ac- 
curacy throughout the temperature range. 

An excellent method for improving the performance 
of such an amplifier, where feasible, is the addition of 
another transistor stage. This would not only achieve 
increased computing accuracy, but also allow greater 
degeneration within the amplifier. The increased de- 
generation would increase the gain margins by reducing 
gain dependence on transistor characteristics. 


VII. CoNncLUSION 


As originally conceived, the purpose of the amplifier 
described was to duplicate the characteristics and per- 
formance of previous vacuum-tube devices. Early in 
the course of the preliminary design, fundamental con- 
siderations of transistor operation precipitated an ap- 
proach which resulted in an amplifier which was ex- 
tremely satisfactory in external performance, but which 
had internal characteristics quite different from those 
initially envisioned. Some preconceived notions about 
input impedance and impedance matching were dis- 
carded, and some natural advantages of the transistor 
in an ac computing amplifier application were dis- 
covered. Methods of division can be devised which do 
not suffer severely from loss of loop gain, and very 
large numbers of summing inputs are practical. The 
outstanding disadvantage is that loop gain varies di- 
rectly with load resistance. 
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A Note on Contact Networks for Switching 
Functions of Four Variables* 
RODERICK GOULDT 


Summary—Several corrections are given to a recent tabulation 
of two-terminal contact networks realizing the switching functions of 
four variables. Nineteen new networks which are more economical 
in contacts than those previously tabulated are also presented and 
certain four-variable functions possessing a useful complementary 
relationship are listed. In conclusion this paper indicates an area 
where further work on four-variable contact network synthesis is 
needed. 


SWITCHING function of four binary variables 
A is definied by mapping the values 0 and 1 onto 

the 24 configurations of values of the arguments. 
It follows that there are 2", or 65,536, different func- 
tions of four variables. These may be divided into 402 
classes such that all of the functions of each class are 
equivalent under argument transformations (7.e., com- 
plementation or rearrangement of the independent 
variables). A list of 402 functions, comprising one repre- 
sentative from each of the 402 classes, was published 
by Aiken and his Harvard associates in 1951 [1]. The 
utility of this list is demonstrated by the fact that a 
vacuum-tube switching circuit for any four-variable 
function may be obtained from the table of 402 cir- 
cuits given in [1]. The expression “function of four 
variables” in this paper refers only to those included in 
the Harvard list. 

In a recent book [2], Higgonet and Gréa included a 
table of two-terminal contact network realizations for 
238 of the 402 four-variable functions; their networks 
were designed to use as few contacts as possible. A 
check of this table has revealed a number of errors. The 
circuits listed for the following functions! are incorrect: 
49(5-19), 68(6-11), 83(6-26), 84(6-27), 158(7-51), 185(8- 
22), 211(8-48), 231(8-68). Correct networks for these 
functions are given in Fig. 1. These new networks were 
obtained by a synthesis method, described by Gould [3], 
which involves using the given switching function to 
obtain a vector space; the latter in turn determines the 
desired network. 

In addition to the errors already mentioned, the fol- 


* Manuscript received by the PGEC, February 3, 1958. A portion 
of the work reported in this paper was done while the author held a 
Bell Telephone Labs. Fellowship at the Computation Lab., Harvard 
University, Cambridge, Mass. 

{+ Deceased. Formerly at Comité d’Etude et d’Exploitation des 
Calculateurs Electroniques, Brussels, Belgium. 

1 In the original Harvard tabulation of the 402 functions, each 
was identified by an integer from the range 0, 1, - - - , 401. The iden- 
tical list of functions was used by Higgonet and Gréa, but different 
identification numbers were assigned. In this paper, each function is 
identified by the Harvard number, followed by the Higgonet-Gréa 
number in parentheses. Functions having a Harvard number greater 
than 237 do not appear in the Higgonet-Gréa table, so for these only 
one number is given. 


lowing misprints occur in circuits of Higgonet and Gréa 
[2]: 78(6-21)—at bottom of network, z’ should be 3; 
96(6-39)—the unmarked contact connected between x 
and gz’ should be labeled y; 107(6-50)—in second x con- 
tact down from top, x should be x’; 136(7-29)—at top 
of network, x should be x’; 199(8-36)—at top of network, 
x’ should be x; 208(8-45)—at bottom of network, w 
should be w’; 209(8-46)—in contact connected between 
z’ and y, x should be x’. It is believed that all of the 
remaining four-variable circuits of [2] are correct. 

For any switching function f, a parameter m, 0<m 
<16, is defined as the number of configurations of the 
argument values for which f=1. Each function f; hav- 
ing m,<8 is the complement of another function fe 
having m,=(16—m). If a planar contact network for 
fi is known, a network for f. can be found directly by 
means of a well-known construction [4], and con- 
versely. Apparently, it is for this reason that Higgonet 
and Gréa gave networks only for functions having 
m<8, although a number of their networks are non- 
planar. 

If the function f; has m=8, the complementary func- 
tion fi’ is in most cases equivalent to f; under trans- 
formations of the arguments. However, fi’ is sometimes 
equivalent to another function fp having m=8. Of the 
74 four-variable functions for which m=8, there are 
16 pairs such that either member of each pair is equiva- 
lent to the complement of the other member. These 
16 pairs of functions are listed in Table I (p. 198), and 
with each pair is given the transformation which must 
be applied to the first function in order to make it com- 
plementary to the second [1]. The facts summarized in 
Table I are apparently not generally realized. For exam- 


- ple, the networks given by Higgonet and Gréa in [2] for 


the functions 175(8-12) and 181(8-18) use 8 and 6 con- 
tacts respectively. Since the networks are planar, it is 
evident that a 6-contact network for function 175 may 
be obtained from that for 181. The networks listed in 
[2] for the following other functions might similarly 
be simplified: 176(8-13), 192(8-29), 193(8-30), 206(8-43), 
213(8-50), 214(8-51), 216(8-53), 218(8-55), 230(8-67), 
232(8-69). 

A second table of two-terminal contact networks, 
covering all of the 402 functions of four variables, has 
been compiled independently by Moore [5]. Moore’s 
circuits are in many cases more economical than the 
correct circuits. However, application of the vector 
space synthesis method previously mentioned has pro- 
duced nineteen networks which require in each case one 
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Fig. 1—Contact networks realizing various functions of four variables. 
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TABLE I 


Distinct Four-VARIABLE FUNCTIONS WHICH ARE COMPLEMENTARY 
UNDER ARGUMENT TRANSFORMATIONS 


Functions Transformation 
172(8- 9) 186(8-23) w's'y'x 
173(8-10) 217(8-54) sw'x'y’ 
175(8-12) 181 (8-18) waxy" 
176(8-13) 194(8-31) yy x'w' 2! 
178(8-15) 188(8-25) w'y x’Z 
183(8-20) 218(8-55) zw'x’y! 
187(8-24) 207 (8-44) x'w'y 2! 
191(8-28) 231 (8-68) w'x y'2! 
192(8-29) 211(8-48) w'y'2'x 
193 (8-30) 221(8-58) 2 y'w's! 
196(8-33) 213(8-50) w'y’Z x 
197 (8-34) 230(8-67) x'2'y w 
199(8-36) 232(8-69) w'x y Z 
200(8-37) 216(8-53) wy 2x 
206(8-43) 223 (8-60) y w' x's! 
214(8-51) 229 (8-66) wexy 


fewer contact than the corresponding network of Moore 
and one or two fewer than that of Higgonet and Gréa 
[2]. These new networks, which are shown in Fig. 1, 
realize the following functions: 68(6-11), 84(6-27), 
96(6-39), 138(7-31), 141(7-34), 147(7-40), 154(7-47), 
158(7-51), 171(8-8), 182(8-19), 268, 271, 277, 284, 304, 
320, 324, 332 and 339.2 Moore [5] has established lower 
bounds for the number of contacts required to realize 
each of the 402 functions. It follows from these bounds 
that each of the circuits of Fig. 1 uses the minimum pos- 
sible number of contacts. 

Exhaustive application of the vector space synthesis 
method has indicated that a number of Moore’s bounds 
may be raised. However, for nine of the 402 functions, 


2 The author has been informed that similar networks for eleven 
of these functions have been found previously: for 182 by S. H. 
Caldwell; for 96, 138, 141, 147, 154, 171, 182, 268, 271, 284, 332 by 
D. A. Huffman; for 96, 141, 171, 271, 332 by E. J. McCluskey, Jr. 
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TABLE II 


Lower BounpDs* ON THE NUMBER OF CONTACTS REQUIRED TO 
REALIZE CERTAIN FouR-VARIABLE FUNCTIONS 


Function Lower Bound 
92(6-35) 10 

160(7-53) 11 

163(7-56) 12 

205 (8-42) 11 

212(8-49) 10 

264 11 

290 11 

293 12 

340 11 


* Bounds due to E. F. Moore. 


the most economical networks now known require a 
number of contacts greater than the presently known 
lower bounds. 

These nine functions with their present bounds are 
listed in Table II. Other circuit designers may be in- 
terested in attempting either to synthesize networks us- 
ing the number of contacts shown in Table II, or to 
raise the bounds. 
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/ Summary—The number of integrators in an analog-computer 
etup should be equal to the order of the differential equation de- 
cribing the system. This paper presents a new procedure for tracing 
ne loop currents which results in one-to-one correspondence be- 
ween the number of integrators in the stimulation setup and the 
ount of independent energy-storing elements in the network, i.e., 
ne degree of the system’s characteristic equation. The generality of 
ne procedure proves that it is always possible to trace the loop cur- 
ents in such a way that excess integrators are avoided. The loop- 
malysis and the branch-variables-analysis approaches are discussed 
md examples given. 


INTRODUCTION 


C NIMULATION of an electrical network on an analog 
S computer can be unsatisfactory because of hidden 
regenerative loops. Such regenerative loops may 
‘rise when the computer setup contains more inte- 
‘rrators than the order of the differential equation de- 
cribing the network. The regenerative loops are re- 
erred to as hidden loops, since the unwanted feedback 
irises because the actual computer components depart 
rom exact, designated values. In other words, the 
-omputer solves a differential equation of a higher order 
than the physical system that is to be simulated; and 
‘he extraneous roots may cause considerable departures 
rom the correct solutions. This problem has been 
yointed out and discussed by Walters! and more re- 
sently has been the subject of a communication by 
Scott.” 
| The danger of excess integrators has been avoided by 
_arrowe,® who developed the so-called “direct simula- 
ion” method. It has recently been pointed out‘ that 
=xCeSS integrators can be eliminated in the loop-analy- 
sis approach and in many cases in the node-analysis ap- 
sroach when the loops, or the reference node, are prop- 
rly chosen. This paper presents a new procedure for 
‘racing the loop currents in such a way that there is one- 


* Manuscript received by the PGEC, February 3, 1958; revised 
manuscript received, June 30, 1958. This work was in part conducted 
yy Project MICHIGAN under Dept. of the Army Contract (DA- 
36-039-SC-52654), administered by the U. S. Army Signal Corps. 

+ Willow Run Labs., University of Mich., Ypsilanti, Mich. 

1L. G. Walters, “Hidden regenerative loops in electronic analog 
omputers,” IRE Trans. oN ELECTRONIC CompuTERS, vol. EC-2, 
»p. 1-4; June, 1953. : 

2N. R. Scott, “On the use of redundant integrators in analog 
omputers,” IRE TRANS. ON ELECTRONIC CompuTERs, vol. EC-6, 
yp. 287-288; December, 1957. ; 

8 V, L. Larrowe, “Direct simulation bypasses mathematics,” Con- 
rol Engineering, vol. 1, pp. 25-31; November, 1954. i f 

4 J. Otterman, “How to avoid extra integrators when simulating 
2LC networks,” Control Engineering, vol. 4, pp. 111-114; November, 
1957. 
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On the Loop- and Node-Analysis Approaches to the 
Simulation of Electrical Networks: 
JOSEPH OTTERMAN+ 


to-one correspondence between the number of integra- 
tors in the simulation setup and the count of independ- 
ent energy-storing elements in the network, 7.e., the 
order of the differential equation describing the net- 
work.’ The procedure is rather easy to apply, even to 
complicated networks or, through mechanical-electrical 
analogy, to mechanical systems. The generality of the 
procedure proves that it is always possible to trace the 
loop currents in such a way that excess integrators will 
be avoided. The node-analysis and the branch-variables- 
analysis approaches are discussed. 


Loop-ANALYSIS APPROACH 


For an example of a computer setup with an excess 
integrator, consider the simulation by the loop-analysis 
approach of the network discussed by Walters.! The 
network is shown in Fig. 1. If the loop currents are 
traced as in Fig. 1(a), the loop equations are: 
pL 


(ox +Ri+ | fete In = Ex), (1) 


1 

pCi 
if 

—pL-h+ (s+ R+—) i= 0 (2) 


where p is the differentiation operator and where the 
initial conditions are disregarded. The computer setup, 
the possible instability of which has been analyzed by 
Walters, is shown in Fig. 2. This setup requires four 
integrators although there are only three energy-storing 
elements in the network. 

When only one loop current is traced through the 
inductance L, as shown in Fig. 1(b), the loop equations 
are: 


1 il 1 
R — )I R R — + — )I,=E; 3 
( oe +( Le a) 2 (p) (3) 
1 
— pie ae (+=) I, = 0. (4) 


The corresponding computer setup, shown in Fig. 3, 
requires only three integrators. The total number of 
amplifiers is reduced from 8 to 6. 

The presentation which now will be given for tracing 
the loop currents to avoid excess integrators is compli- 


5 J. Otterman, “On the order of the differential equation describ- 
ing an electrical network,” Proc. IRE, vol. 45, pp. 1024-1025; July, 
1957. 
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(a) The loop equations contain terms with , that range from 
—1 to +1 in powers of p; and terms with J2 that range from —1 to 
+1 in powers of p. 


(b) The loop equations contain terms with J; that range from —1 
to +1 in powers of p; and terms with Jz that range from —1 to 0 
in powers of p. 


Fig. 1—Different tracing of loop currents will result in eliminating 
an excess integrator in the computer setup. 


cated by the fact that its aim is to demonstrate the 
generality of the procedure. The actual rules for tracing 
loop currents in such a way that excess integrators will 
be avoided [as was done in Fig. 1(b)] can be gleaned 
from this presentation by skipping all discussion per- 
taining to the count of nodes, elements, and separate 
parts. 

In the following, a node is defined asa junction of two 


or more elements, a junction of two ends of the same . 


elements, or an isolated end of a single element. All 
the current sources in the network are assumed to have 
been converted into voltage sources. The appearance of 
a voltage source in any loop does not affect the problem, 
and voltage sources are disregarded in the following dis- 
cussion. Mutual inductances are not counted as elements 
and do not require any special consideration. 

The procedure may best be described as an instruc- 
tion to redraw the network according to certain direc- 
tives. First draw all the resistances in the network. If 
any closed loops have been formed, trace loop currents 
through them. The number M, of the loops required 
will be given by the formula 


M, = B,— Nr +3; (5) 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


September 


Fig. 2—Four integrators are required if loop currents 
are traced as in Fig. 1(a). 


Fig. 3—Only three integrators are required if loop currents 
are traced as in Fig. 1(b). 


where B,, N, and S, are, respectively, the number of 
elements, the number of nodes, and the number of 
separate parts of the resistive network.® Since the loops 
contain resistances only, the loop equations are alge- 
braic and no integrators are necessary to represent them. 

Then, element by element, draw all the inductances 
and capacitances in an arbitrary order. During this 
process of increasing the partial drawing of the network, 
trace a loop current whenever a loop has been formed. 
An element or elements through which the loop current 
is traced for the first time, z.e., an element or elements 
which are being added to the loop structure, will be 
referred to as “new” elements. 

The formation of a new loop may occur in several 
topologically different ways, as is shown in Fig. 4. 
However, in each case the count of independent loops 
M of the partially drawn network, defined by 


M=B-N+S, (6) 


where B is the number of elements or branches, N the 
number of nodes, and S the number of separate parts 
in the partially drawn network, increases by only one. 
Thus, if we start with M, loops, given by (5), and trace 
one loop current each time a new loop is formed, we 
have always traced in the partial drawing of the net- 


5 The simplest counting procedure is to disregard all the elements 
already drawn which do not form part of closed loops. The count of 
M, by this procedure is identical to M, arrived at by counting the 
number of resistances, nodes, and separate parts in the resistive 
network, even though the counts of B,, N, and S, are different. 
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Fig. 4—The addition of new elements can occur in topologically dif- 
ferent ways. The solid lines denote the previously drawn loop 
structure, the dotted lines denote the new elements. The 
rectangles denote an element or series of elements. In each case 
the count of the loops M=B—WN-+-S is increased by one. 


work the number of loops given by (6). At the time we 
trace the loop current through the last new element 
(the addition of which completed the network), the 
right side of (6) becomes B,—N.+S., where B., NV. and 
S, are, respectively, the number of elements or branches, 
the number of nodes, and the number of separate parts 
of the complete network; and we have traced the neces- 
sary number of M, loop currents. Moreover, the loop- 
current equations form an independent set. This is so 
because each loop current is passed through new ele- 
ments; the equation for the voltages in this loop is there- 
fore independent of equations written for the previously 
drawn loops structure. 

The loop current is traced through new elements and 
lexcept in the case where new elements by themselves 
form a closed loop, as shown in Fig. 4(b) and Fig. 4(c) | 
through elements through which other loop currents 
have been traced. The completion of the current through 
the previously existing loop structure is done in such a 
way as to minimize the range in powers of p of voltage 
terms associated with the new loop current. When the 
loop current passes through L and C, or through L, R, 
and C, the expressions for voltages associated with the 
current include terms divided by p and multiplied by 
p; and simulation requires two integrators. When the 
loop current passes through L and R, or R and C, the 
range in powers of p of the associated terms for voltages 
is one; simulation requires one integrator. When the 
loop current passes through like elements, the terms for 


Otterman: Simulation of Electrical Networks 


201 


associated voltages exhibit the same power in p (p71 
when the loop consists of capacitance, p® when the 
loop consists of resistances, and p! when the loop con- 
sists of inductances), and simulation requires no inte- 
grators. 

The procedure for minimizing the range in powers of 
p for voltage terms associated with each new loop cur- 
rent treats separately seven cases, depending in part on 
what new elements are being added and in part on the 
possible choices of completing the loop through the 
previously traced loop structure. Examination of con- 
secutive cases is required. It can be shown that in each 
case the additional number of integrators required is 
equal to the increase in the number of independent 
energy-storing elements. 


Case I. Addition of L Elements Only 


The new elements include inductances only. If it is 
possible to complete the loop through the previous 
loop structure in such a way that only inductances are 
traversed, then trace the loop structure through this 
path. The terms for voltages associated with this loop 
current exhibit the same power in p. No additional inte- 
grator is required in the computer analog of the network. 
This is in accordance with the fact that the added induct- 
ance forms part of an inductances loop and is therefore 
dependent in terms of Kirchhoff’s law of voltages.® Its 
addition, therefore, does not increase the order of the 
differential equation of the drawn loop structure. If 
several inductances are being added, all but one are 
dependent in terms of Kirchhoff’s law of currents as 
applied to the partially drawn network, so the increase 
in the order of the differential equation is still zero.® 

If a closed loop through inductances only cannot be 
found, consider Case II and Case III. 


Case II. Addition of L or Rand L Elements 


The new elements are inductances, or inductances and 
resistances. If a closed path through the previously 
drawn loop structure can be found that does not traverse 
any capacitances, then trace the loop current in such 
a way. Simulation of the current will require one inte- 
grator. This is in accordance with the fact that the 
addition of an independent inductance increases the 
order of the differential equation of the loop structure 
by one. If several inductances are being added, all but 
one are dependent in terms of Kirchhoff’s law of currents 
as applied to the partially drawn network, so the in- 
crease in the differential equation is still one. 

If a path that does not traverse any capacitance 
cannot be found, Case III applies. 


Case III. Addition of L or Rand L Elements 


If a path that does not traverse any capacitance 
cannot be found, then trace a loop that includes a 
capacitance or capacitances. The new loop passes 
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through at least one inductance and one capacitance; 
thus the simulation requires two integrators. The capac- 
itances in the previously drawn loop structure were not 
independent in terms of Kirchhoff’s law of currents as 
applied to the partially drawn network, since a cut-set’ 
of capacitances existed.’ The addition of the new ele- 
ments provides an additional current path, eliminating 
a cut-set of capacitances and converting one dependent 
capacitance to an independent capacitance. The in- 
crease of the order of the differential equation is two, 
since the counts of independent inductances and inde- 
pendent capacitances are increased each by one. 


Case IV. Addition of C Elements Only 


Read Case I, substituting the word “capacitance” for 
“inductance.” 
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If the loop currents have been traced in accordance 
with the above procedure, the loops for stating Kirch- 
hoff’s law of voltages can be chosen in a different way 
than the loops of currents. In other words, adding and 
subtracting of loop equations does not affect the prob- 
lem of minimizing the number of integrators. 

Examples of the procedure for tracing the loop cur- 
rents are given in Figs. 5 and 6 (pp. 203-205). 


NopE-ANALYsIS APPROACH 


The node-analysis approach, too, quite often results 
in a higher number of integrators than the degree of 
the system’s characteristic equation. Walters! analyzes 
the network shown in Fig. 1, with the node-voltage 
variables as shown in Fig. 7(a). The node equations 
are: 


1 Ex(p) 
—+9C,)Vi- Cr Vet Ds (7) 
Ry Ky 
1 
ea AG tie Vit (So + ati + C2) Y2— pCe2 V3=0 (8) 
1 
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Case V. Addition of C or Rand C Elements 


Read Case II, substituting the word “capacitance” 
for “inductance,” and vice-versa. 


Case VI. Addition of C or Rand C Elements 


Read Case III, substituting the word “capacitance” 
for “inductance,” and vice-versa. 


Case VII. Addition of L and C or L, R, and C Elements 


When both Z and C are added, their addition cannot 
increase the number of loops of capacitances or induct- 
ances, as happens in Cases I and IV. The order of the 
differential equation is increased by two, since the 
counts of independent inductances and independent 


capacitances are increased each by one. Two integrators — 


are required. The loop can be completed through the 
loop structure, which was previously traced, in an ar- 
bitrary manner. 


The corresponding computer setup, shown in Fig. 8, 
contains four integrators. 

Excess integrators can be avoided by using the node- 
analysis approach if the network is an LC continuous 
network,‘ that is, a line can be traced from a node to 
every capacitance in the network by traversing capaci- 
tances only, and a line can be traced from the same 
node to every inductance in the network by traversing 
inductances only. Such a node should be chosen as the 
reference node. In an LC continuous network with 
only capacitances (only inductances) and resistances, 
any node adjoining a capacitance (an inductance) and a 
resistance can be chosen as the reference node. 

The network analyzed by Walters is actually an LC 
continuous network. To avoid the danger of excess inte- 
grators, the node adjoining inductance Z and capaci- 
tances C; and C2 should be chosen as the reference 
node. If the node-voltage variables are chosen as in Fig. 
7(b), the node equations are: 


1 1 Exp) 
& re ) } Ry, yr 0 R, (10) 
1 1 1 1 Exp) 
—-— -Vit (— — =) V — =— 
1 : L a Ry 4 Rs : Rs» : Ry (11) 
1 
Ones Vee (=+ pcs) V3\= 0. (12) 


7 A cut-set is a set of elements such that when each element in the 


set is cut into two, the network is separated into two parts, 


The computer setup, shown in Fig. 9, contains only 
three integrators. 


(a) Simulation of the loop currents requires four integrators if 
the loops are traced as above. 
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RI R2 


R3 


(b) The resistive network. No loops formed. 


(c) Addition of Z;, Re, R3 to the loop structure. The new elements 
by themselves form a closed loop. 


Ri 


(e) Addition of Z2. Case II applies. The simpler of two possible 


loop currents was chosen. 


(d) Addition of C, Ri. Case V applies. 


(f) Simulation of the loop currents requires three integrators if 
the loops are traced as above. 


Fig. 5—Tracing of loop currents through a simple network. 


The network shown in Fig. 5 is an LC continuous 
network, and in the node analysis the node common to 
C, Li, L2 and R3 should be taken as the reference node, 
The network shown in Fig. 6 is not LC continuous. 


BRANCH-VARIABLES-ANALYSIS APPROACH 


Both the node analysis and the loop analysis are 
systematical approaches based on the fundamental 
relations in electrical networks. Those fundamental 
relations consist of 2B equations, where B is the number 
of branches or elements. The equations contain 2B un- 
knowns, which are the current through each branch and 
the voltage across each branch. B of the equations relate 
voltage across a branch to the current through the 
branch. M equations relate voltages around M inde- 


pendent loops. N-S equations, where N is the number 
of nodes and S the number of separate parts, relate the 
currents through N-S independent node-pairs, or, more 
generally, through N-S independent cut-sets. 

Analyzing the network in terms of those basic vari- 
ables is more laborious, because the number of vari- 
ables is larger; but it is, in some respects, much more 
flexible. 

For an example of the branch-variables-analysis ap- 
proach, consider the simulation of the network shown 
in Fig. 1. There are five elements in the network, and 
the number of variables is therefore 10. The positive 
direction of currents and voltage-drops is taken down- 
wards and from left to right. The five equations relating 
voltages to currents are: 
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(a) The network to be simulated contains 14 energy storing ele- (b) The resistive network. No loops formed. 
ments. There is one cut-set of inductances, L3, L4, one cut-set of capaci- 
tances, C2, C3, C1, Cs, Cz, one loop of inductances Ls, Ls and one loop 
of capacitances Cs, Cs, Cz, Cs. Order of the differential equation is 10. 


Ri 


(c) Addition of Cs, R; to the loop structure. The new elements by 
themselves form a closed loop. 


Ri 


(e) Addition of Z;. Case I applies. (f) Addition of Cs, Cz, Re. Case V applies. The current could be 
traced through C3 instead of Rs. 


Ri 
(g) Addition of C;. Case IV applies. (h) Addition of Ls, Cs, Ri, Ls. Case VII applies. Loop current 
can be traced in an arbitrary way through the previously drawn loop 
structure. 


Fig. 6—Continued on next page. 
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i@ Addition of Lz, C3. Case VII applies. Loop current can be traced in (j) Addition of Cx, Cy. Case V applies 
| 5 Cite ; 


an arbitrary way through the previously drawn loop structure. 


L4 G7. L6é 


(kk) Addition}of Z;. Case III applies. Loop current can be traced in 
an arbitraryfway through the previously drawn loop structure. 


Fig. 6—Tracing the loop currents through a more complicated network. 


V1 V2 V3 V1 ) V3 


0 
(a) The node equations contain terms with V; that range from 0 (b) The node equations contain terms with Vi that range from 
to 1 in powers of p; terms with V2 that range from —1to1inpowersof 0 to 1 in powers of p; terms with V2 that range from —1 to 0 in powers 
p;.and terms with V; that range from 0 to 1 in powers of p. of p; aud terms with V; that range from 0 to 1 in powers of p. 


Fig. 7—Different choice of reference node will result in eliminating an excess integrator in the computer setup. 


pas 2 Vea (13) The three equations relating currents through cut-sets 
Pai jes are: 
if = = 
ee for (14) giao 0 (18) 
PCr Ice — Ir2 = 0 (19) 
I 
Veo = ee (15) fp ples Ten= "0: (20) 
pC2 
Vre The two equations relating voltages in loops are: 
In = = (16) 
V, Ex’) + Veit Vi+ Ver = 0 (21) 
L 
i= (17) 


pL Voo t+ Vaz — Vi = 0. (22) 
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Fig. 8—Four integrators are required if the reference node 
is chosen as in Fig. 7(a). 


-V2/pk 


Fig. 9—Only three integrators are required if the reference node 
is chosen as in Fig. 7(b). 


The corresponding computer setup is shown in Fig. 10. 
This setup could also be obtained by applying the tech- 
nique of Larrowe.? 

To avoid extra integrators when a cut-set of induct- 
ances or capacitances occurs, the equation expressing 
the voltage across one element in the set in terms of the 
current through this element should be replaced by an 
equation expressing the same voltage in terms of volt- 
ages across other elements in the set. For instance, in 
the cut-set of Z3; and LZ, in the network shown in Fig. 
6(a), the element JL; is arbitrarily chosen as a dependent 
element and the equation for the voltage across this 


element 
Viz = plat zs (23) 


is replaced by 
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Fig. 10—The network of Fig. 1 can be simulated by this setup in 
which all the branch-variables are available for recording p urposes. 


L3 
Vi3= —- pPLsl ra = —— Vy. (24) 
Ly 


Conversely, when a loop of inductances or capaci- 
tances occurs, the equation expressing the current 
through one element in the set should be replaced by an 
equation expressing the same current in terms of cur- 
rents through other elements in the set. For instance 
[Fig. 6(a) ] the equation for the current 


Vie 
Its = — (25) 
L6 Jae 
is replaced by 
V Lek 
ype foe (26) 
ple Le 


Using this technique, no integrators are associated 
with the energy-storing elements chosen as the depend- 
ent elements. 


CONCLUSION 


The loop-analysis, node-analysis, and branch-vari- 
ables-analysis approaches have been reviewed. The 
choice of the most suitable approach depends on many 
factors, such as the relative number of loops and nodes, 
the presence or absence of mutual inductances, and the 


desirability of representing voltages across elements or 


voltages relative to a certain point in the network. 
Therefore, familiarity with the different approaches is 
distinctly advantageous. 


CoS 


| Summary—The usual definition given for the parity check is 
| unwieldy and not particularly suited for the analysis or the study of 
| the arithmetic properties of parity checking. The definition of parity 
by means of congruences provides a convenient mathematical basis 
‘for the concepts of the parity check. In this paper congruence nota- 
| tion is used to generalize the concepts of parity to include nonbase 
| two number systems. Consideration is given to the cases where the 
check base is equal to the number base, and where it is not equal to 
the number base. The arithmetic properties of each case are con- 
| sidered by means of congruences. 


INTRODUCTION 
He ee CONCEPT of the parity check! is widely 


known and has frequently been used to detect and 

in some cases to correct errors in digital systems. 
The parity check is obtained by associating with each 
_ representation of information an auxiliary digit or group 
of digits known as the parity check digit(s). The proper 
choice of the relationship between the check digits 
and the information digits permits the establishment of 

_ mathematical operations, which may indicate erroneous 
digits or provide correction data. Parity checking 
usually is defined for the binary system in terms of the 
_ odd or evenness of the number of ONE digits in a speci- 
fied block of binary digits. If the number of ONE digits 
is even the parity check digit is a zero. If the number of 
ONE digits is odd the parity check digit is a one. This 
definition is sufficient for most practical applications of 
the parity concept, but gives no insight into the funda- 
mental principles of parity. A restatement of parity in 
terms of linear congruences provides a mathematical 
basis for the parity concept. In terms of congruences, 
the parity check may be defined for the general case 
where the check base does not equal the number base. 

_ The parity check which employs a check base equal to 
the number base (as in the usual definition given above) 
is a special case of the general parity check. The con- 
gruence notation also facilitates the study of the arith- 
metic properties of the parity check. The properties of 
the parity check with regard to arithmetic operations 
has received little attention in the literature of error 
correcting and detecting codes. Previous research in 
this area is due to Block et al.? and Robertson.‘ 


* Manuscript received by the PGEC, March 10, 1958; revised 
manuscript received, June 23, 1958. 

+ University of Michigan, Ann Arbor, Mich. 

1R, W. Hamming, “Error detecting and correcting codes,” Bell 
Sys. Tech. J., vol. 26, pp. 147-160; April, 1950. 4 

2R. M. Block, R. V. D. Cambell, and M. Ellis, “The logical 
design of the raytheon computer,” Math. Tables and Aids to Comp., 
vol. 3, pp. 286-296, 317-323; October, 1948. 

3 R. Block, U. S. Patent No. 2,634,052. : - 

4J. E. Robertson, “Error Detection and Correction in Binary 
Parallel Digital Computers,” Digital Computer Lab., University of 
Illinois, Urbana, Ill., “Electronic Digital Computer” Internal Rep. 
No. 37, pp. 70-81; 1952. 
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Robertson studied the application of a Hamming error 
correcting code in a general-purpose digital computer. 
The papers by Block describe a particular error detect- 
ing scheme employed in the RAYDAC computer. 

The primary contribution of this paper is the unifica- 
tion of the parity check concepts within the framework 
of congruences. 


CONGRUENCES 


The remainder of this paper makes considerable use of 
congruence notation. A short presentation of pertinent 
concepts and relationships is included at this point. 
Additional material on the subject may be found in any 
standard text on number theory. 

The congruence 


A =a Mod db (1) 


is read “A is congruent to a modulo b.” The congruence 
relationship states that 


A=a-+bi (2) 


is valid for some value of ¢, where A, a, b and ¢ are inte- 
gers. a is called the residue and 6 the base or modulus of 
the number A. 

As examples of congruences, consider 


10 = 7 Mod 3 (3) 
10 = 4 Mod 3 (4) 
10 = 1 Mod3. (5) 


In these examples the integers 7, 4, and 1 form a resi- 
due class of 10 Mod 3. Of particular importance is the 
least positive residue of the class which in this example 
is one. The least positive residue is that residue for 
which 0<a<d.§ 

Unless otherwise specified, the term residue as used 
in this paper means the least positive residue. 

For the most part, congruences may be treated in the 
same manner as the equality sign in ordinary arithmetic. 
The main exception pertains to division. The properties 
of congruences used here are presented below without 
proof. 


1) Congruences to the same modulus may be added 
and the result is a valid congruence. 


n 


DAs = (2 as) Mod 6 


i= 


(6) 


5 G. H. Hardy, and E. M. Wright, “An Introduction to the Theory 
of Numbers,” Oxford University Press, London, Eng.; 1956. 
6 The equality sign may exist on only one side of the inequality. 
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2) Congruences to the same modulus may be multi- 
plied and the result is a valid congruence. 


n 


IT 4: = (TI as) Mod. (7) 


pil 


3) Congruences are transitive. If d=B and B=C 


then A=C. 
4) The following congruence relationships provide 
a basis for numerical checking procedures 


ob” = a Mod (6 — 1) (8) 
ob” = +o Mod(d+1) neven (9) 
= —oMod(d+1) xn odd. (10) 


Il 


GENERALIZATION OF PARITY 


A parity check need not necessarily include every 
digit of the number. In the general case, multiple parity 
checks may be employed, each check concerning only 
specific digits. However, for purposes of simplification, 
only the case of a single parity check over all digits is 
considered here. The parity of a number from a con- 
sistently weighted number system is defined now. Con- 
sider the number 


(11) 


This is the usual shorthand notation for the polynomial 
representation of a number N in a consistently weighted 
number system, base b: 


N = ondb™ 3 + on_1b”™? + ++ + + 0261 + 010°. 


OnOn—-1 ° © * 6201 


(12) 


A general parity check digit p associated with number 
N is defined by 


F(N) = p Mod m (13) 


where p is the least positive residue. 
This paper considers two different functions of JN. 
The first function is 


F(N) = N (14) 
and the parity check digit p is given by 
N = p Modm m # b. (15) 


This type of check is based directly on the numerical 
properties of linear congruences and the check digit is 
some function of the magnitude of the number JN. This 
type of check is called a numerical parity check. 

If the check base m is identical to the number base 
b, then 


N =o; Mod 6. (16) 


The check is unsatisfactory because it is a function 
of only the low order digit of the number representa- 
tion. This is due to the fact that 


a0” = 0 Mod b n> 0. C7) 
A check for the case m=b is obtained if 
TEGO) =O grab Gpetoateetes 4 abe Cee a Gls (18) 
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The parity check digit p for this type of check, which 
is termed digital checking, is defined by the relation 


> 0: = p Mod b. 
t=1 


An equivalent definition may be given in terms of ring 
addition. Ring addition modulo 0 is specified by the 
symbol @. Ring addition is simply normal addition 
without carry, for example, 


(6—1)@1=0. (20) 


In terms of ring addition modulo 0 the parity digit p is 
given as 


P= On @ Ona OY Dor O or. (21) 


If the number system is binary, and the check base 
is two, then the previous definition for the digital check 
provides a parity digit which makes the total number 
of ONE digits of the number plus check representation 
even in correspondence with the usual parity proced- 
ures. The definition is extended easily to an odd parity 
check. For this case the parity check digit is given by 
p, the diminished radix complement of p defined as 


P@p=b-1. (22) 


As an example of the parity definition for m=), con- 
sider the following. 
Given 12894, a number in base 10, then 


(1+2+8+9-+ 4) = p Mod 10 (23) 


Or 


P=1620869604=4. (24) 


Both types of parity checks are restricted in appli- 
cation to the detection of certain classes of errors in 
the number representation. Error correction is only 
possible if multiple parity checks are employed.7:$ The 
digital parity check, m=), detects all alterations except 
the alterations in which the sum of the digit perturba- 
tions is congruent to 0 modulo b. Let the decimal num- 
ber of the previous example be corrupted by +5 in 
the first column and by —5 in the third column. The 
result is an erroneous representation but the error is 
obviously undetectable. 


1620309@9=p=4. (25) 


A numerical parity check Mod m, m#b, is insensi- 
tive to an error if the magnitude of the error is con- 
gruent to zero Mod m. Thus the error is detectable only 
if the check base m is not a divisor of the magnitude of 
the error. The magnitude of the error in the previous 
example is 495. The error is not detected by Mod 3, 5, 
9 or 11, checks. 

Parity checks are particularly effective in detecting 


’M. J. E. Golay, “Notes on digital coding,” Proc. IRE, vol. 37 
p. 657; June, 1949. : es 

8M. J. E. Golay, “Binary coding,” IRE TRANS. ON INFORMATION 
Tueory, No. PGIT-4, pp. 23-28; September, 1954. 
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rrors for which only one digit has been modified. For 
his class of errors a parity check m=6 is always effec- 
ive. If m is less than 6 then perturbations congruent 
0 zero modulo m are undetectable. 


ARITHMETIC INVARIANCE OF THE DIGITAL 
PARITY CHECK 


In the arithmetic system employed in most digital 
machines, a one-to-one correspondence is established 
between the machine digit representation and a set of 
positive and negative fractions. The rules of machine 
arithmetic operations are set up in such a way as to 
preserve this one-to-one correspondence under the oper- 
ation of addition. Addition is the fundamental arith- 


metic operation of the machine and the other arith- 
‘metic operations, subtraction, multiplication and di- 
vision, are defined in terms of addition and the auxiliary 
‘operations, shift and complement. 

It is obvious that the digital parity check is invariant 
‘to the shift operation if no digits are deleted or added 
‘to the representation. This is true because the digital 
parity check is a function of the modulo sum of the indi- 
vidual digits rather than a function of the magnitude of 
the representation. Also, the parity check is independent 
‘of the position of the radix point. For this reason the 
‘remaining discussion does not explicitly consider wheth- 
er the representations are in correspondence with frac- 
tions or integers. 

- Consider now the addition of two consistently 
weighted numbers in base } representation. The check 
base, of course, is also equal to b. 


| AtD=S 

: nase 

ee g 

he Aan Og 

Now 

| S=aO0d, OG (26) 

and 
(a; + di + cs) = sit e410 (27) 
a1 =0 (28) 
c¢=O0ore;=1if1A1; (29) 


c; is the ith carry digit from the ith minus one column. 
The parity check digit p for the sum is given as 


Cpr sn © Snr D Een = ps (30) 

or 
Cnz1 B On Ban DnB ++: Ou Pd = ps, (31) 

but 
Qn @>:: Oa = pa (32) 
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and 
dn, B+++ Odi = pp; (33) 
therefore, 
Cn4i BD +++ 62 B pa @ pp = pg. (34) 


The last expression provides a means of checking the 
addition operation in terms of the parity check digits of 
the sum, the addend, and the augend and the generated 
carries. Alternatively, the parity check may be in terms 


of p’, the radix complement or additive inverse of p 
defined as 


py ep=d. (35) 
If p’ rather than # is used as the check digit, the check 
procedure for addition is given as 


Cn +++ Oa @ fs’ = pa' © pp’. (36) 


The class of addition errors detected by the digital 
parity check is the same as the class of representation 
errors defined in the previous section. However, it must 
be remembered that the logic of the usual addition 
circuits is such that the various carry digits are not in- 
dependent. The parity checking procedure is valid only 
if the carry digits are correct. 

Consider the decimal addition of the following num- 
bers: 


Olw an 


Oe? 
2 4 
See) 


Due to the occurrence of one carry, the check of addition 
is 


685@1=2 


2=2. (37) 


The check of the complement operation is straight- 
forward. The diminished radix complement of a repre- 
sentation A is obtained by substituting for each digit 
a; the diminished radix complement, 4. 


a@ag=b-1 (38) 
given 
A= Gy 2204 and A = G, = = 'ay (39) 
pa = an @D Ga, and pa=a,@ -@ dG (40) 
and 
On Din B+: OuyGH= ps pa=f. (Al) 
Then 
f =n(b — 1) Mod b (42) 
(f + n) = nb Mod b (43) 
(f + n) = 0 Mod b (44) 
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and 

pa t+ patn=0 Moda, (45) 
or 

pa © pa B n Mod d = 0. (46) 


A check of the diminished radix complement consists 
of the modulo b sum of the parity check of the normal 
representation, the parity check of the complemented 
representation and the number of digits comprising 
the representation. In the special case where b=2, the 
check of the complemented representation is equal to 
the check of the normal representation if 2 is even, and 
opposite if 2 is odd. 

A check of the radix complement operation is some- 
what more complex than that required for the dimin- 
ished radix complement. The complications are due to 
the definition of the radix complement rather than to 
the nature of the parity check procedure. Two checking 
procedures are possible. The first is based on the fact 
that the radix complement A’ is obtained from the 
diminished radix complement A by the addition of one 
to the lowest order digit. This method requires cogni- 
zance of the carries produced when the one is added to the 
lowest order digit. The check is given by the relation- 
ship 

n+1 


pa = 10 p= ® Dc: Mod d. 


1=2 


(47) 


The second procedure for obtaining the radix com- 
plement check is based upon the well-known rule for 
obtaining the digitwise radix complement: “Starting 
with the low order digit, radix complement the first 
non-zero digit; the remaining higher order digits are 
changed by diminished radix complementation.” This 
method requires a knowledge of the number of zeros 
preceding the first nonzero digit. Let this number be gq. 

Given 


A= Witn Te o04a10g “8 101 (48) 
PUES. 2s), wigs Mtge (49) 
a, = O17 = @ (50) 
so 
PA = On OD Ags O Gg (51) 
Pan = Oyun @ Oat. Did o11 (52) 
bi Ope as Ones ae 
®@ dare B aqr1 @ a’ q41, (53) 
pa ® par = [(n —g — 1) — 1) +8] Mod (54) 
or 
pa © par © (n — g) Mod db = 1. (55) 


Multiplication consists of repetitive addition with 
appropriate shift operations. The result is a double- 
length product. The digital parity of the double-length 
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product is found from consideration of the addition 
operations since parity is not affected by the shift opera- 
tion. The linear congruence is multiplicative. The check - 
procedure for the product P=A XD is 


Po = (pa © po) ® YY ci Mod b (56) 


or 


(pa X po) + do cs = pp Mod 6. (57) 


The symbol © indicates multiplication, modulo }b. It 
is observed that the digital parity check provides no 
protection against faulty shift operations. Parity pro- 
tection against this type of error is obtained if multiple 
checks are employed. 

The division algorithm commonly employed is easily 
explained in terms of the remainder theorem. 

Consider 


dE R 


BO eh re (58) 
or 

X=QV+R (59) 
Ria OY (60) 

where 
QO = qb + --- + gb. (61) 

Then 
R= X — (qb +--+ + ob 7)Y. (62) 


The subtraction of Y from X, q; times is accomplished 
by the addition of Y, q; times. Carries which are gener- 
ated must be considered. The parity check expression 
is given as 


pr = px ® (po © pr) ® Dc Mod b. (63) 


This check may be employed for each division opera- 
tion or only as a check of the over-all division process. 
In any case the check is insensitive to shift errors. 


‘Examples of Arithmetic Checking (m=b= 10) 


Addition 
3 FG TES) © the circled digits 
+4 5 8 _® are the check digits 
BO): Phot ae) 
Check 
pa ® po ® Dc; Mod b = 4. (64) 
3@07@G2=2 (65) 
Complement 


nines complement 
= 0345@ 
A = 9654@ 


1958 


Check 
pa B pa @ n Modd = O 
46264=0 
ten’s complement 
A = 034500 @ 
A’ = 965500 © 


(66) 
(67) 


/ Check 
par © pa © (n — g) Modb = 1 
9@2@4=1 


(68) 
(69) 
|. Product 


346 ® 
213 © 
346 
346 
692 
346 
1038 

346 

“4498 

346 

39098 

346 

73698 @ 


5 carries 
generated 


Check 
bp = (ba © Pv) © Docs Mod B 
3=(306005 
3=8@5=3 


(70) 
(71) 


Division 
. 
: 073699 ® xX 
) ~~ 0346 mes ‘3 
nis ten’s complement of 0346 © is 9654 @ 


Carries 
073699 2 
9654 Carries 

—— 

039099 2 0347 3 
9654 9654 
004499 0001 
9654 9654 
969899 9655 
04499 3 001 
01039 Fs 
9654 fe a 
97579 De co= 
1039 1 
9654 
0693 2 
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Check 
pr = px ® (pe © Pr) ® Dc; Mod bb 
1=40(604)63=46463=1. (72) 


NUMERICAL CHECKING, Mop m(m#b), oF ARITH- 
METIC OPERATIONS 


The modulus of the numerical parity check is not the 
same magnitude as the radix of the number representa- 
tion. Congruence relationships are valid for any modu- 
lus. However, only certain moduli yield procedures 
which permit straightforward computation of the parity 
digits. There exist in particular two interesting and 
useful checks based on procedures congruent modulo 
(6—1) and congruent modulo (6+1). These checks are 
called the diminished base numerical check and the 
augmented base numerical check, respectively. The 
simplest checks of the numerical class depend directly 
on certain congruence relationships. Numerical check- 
ing procedures have one very desirable feature not 
associated with digital checking procedures; they do not 
require a tally or cognizance of the carries generated by 
the addition process. 


THE DIMINISHED BASE NUMERICAL CHECK 


The diminished base numerical check is dependent 
on the following properties of linear congruences: 


a0" = o; Mod (6 — 1). :(73) 


This relationship permits parity calculation independent 
of the digit position or weight. Consider the sum 


xXx+Y=S 


X= gb) Se ab? (74) 
VY = y,b"1 + -+- + 10° (75) 
S = Cn41b" + 5,671 + - +--+ 516%. (76) 
Since congruences are additive, 
X = >> 2b! = DS x; Mod (6 — 1). (77) 


t=1 t=s1 


The parity check digit px associated with the represen- 
tation X is defined by the term 


n 
Di 
t=1 


reduced to the least positive residue of the class. The 
parity check digit px also may be defined by 


xt in De OHO (78) 
for either definition 
X= s xsb*-! = px Mod (6 — 1). (79) 
i=l 
The check digit py is given as 
Y= = yib*-! = px Mod (6 — 1). (80) 


i=l 
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Linear congruences are additive, therefore, 
X+V =D) (e+ yb"? = (bx + py) Mod (6 — 1) (81) 
u—\ 


and 


ps = px © pr. (82) 
The modulo (6—1) check is numerical and hence 
does not require cognizance of the generated carries. 
Nevertheless, the check is sensitive to errors in carry 
generation and propagation. The check detects all digit 
alterations for which the modulo (6—1) sum of the 
changes does not equal zero. This is equivalent to the 
statement that an alteration is-detectable only if b—1 
is not a factor of the magnitude of the arithmetic error. 
A single carry failure causes an alteration of minus one 
and is therefore detectable. However, the error is not 
correctable on the basis of the parity information. 

A check of the complement operation is obtained 
from the basic definition of the complemented represen- 
tation. The definitions which follow are for representa- 
tions of n+m-+1 digits. The representation has n digits 
to the left of the radix point, m digits to the right of 
the radix point and one sign digit. The diminished 
radix complement of X is X where 


X+X = pti —p™, (83) 
The radix complement X’ is given as 
X +X = 3, (84) 


The usual congruence relationships are defined for 
integer values. It is, therefore, convenient to consider 
the radix point on the extreme right of the representa- 
tion. The check is obtained by observing that 
b”("t! — 6-”) = 0 Mod (6 — 1) (85) 

and 
b”b"+l = 1 Mod (6 — 1). (86) 


Then the check for the diminished radix complement is 


px ® px = 0 (87) 
and the check for the radix complement is 
pre@ pros. (88) 


Congruences are multiplicative. Therefore, a check of 
the product 
Pi ROS 
is given simply as 
pe = prO Ps (89) 


with no cognizance of carries required. 
A check involving the division of congruences is not 
practical since the correct division procedure is depend- 
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ent upon whether the divisor and (b—1) have any com- 
mon factors. However, this is of little consequence since 
machine division is usually obtained by a process con- - 
sisting of repetitive add and shift operations. 


xX a R (90) 
yee ee 
X=OQOV+R. (91) 


The check of the division process is independent of 
carries and must be true at every step in the division 
process. The check is 


px = (po © py) @ hr. (92) 


THE AUGMENTED BASE NUMERICAL CHECK 


The diminished base numerical check avoids the 
need for carry cognizance and for that reason offers a 
simpler check procedure than is obtainable by means of 
digital checking. 

However, for the binary number system the dimin- 
ished base check is meaningless. One solution is the 
augmented base check. This system is also carry inde- 
pendent, and parity check digits are relatively easy to 
obtain. The augmented base check is dependent on the 
following property of congruences. 


a:b? = + o; Mod (6 + 1) if7 is even (93) 
=! sg; Mod (baits atas odd aeenOn) 
The check digit for representation A is given as 
(+++ as — ag + a3 — a2 + a1) = p Mod (641) (95) 
where 
A = +++ + a5b4 + a4b? + a3b? + aed! + a1b°. (96) 


The checking procedures employed for the arith- 
metic processes are similar to the procedures for the 
diminished radix except that ring addition @ and ring 
multiplication © are defined modulo (6+1). 

In the binary case the number base is two and the 
check base is three. The parity of a representation A 
is given as 


(+++ + as + 2a, + a3 + 2a. + a1) = p Mod3 (97) 


or 


7+ + ® as © 2as © ag © 2ae © ay = P. (98) 
Direct subtraction is avoided by using the radix com- 
plement of one with respect to the base three. 

Note that for the binary case the augmented base 
check is identical to the diminished base check, base 
four. 


EXAMPLES OF THE AUGMENTED BASE CHECK 


Let b=2; then the check is modulo 3. 
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Addition X+Y=S 
1 
1 


Or0e0 ste te1-. 0 Check 
Ont O20 510 bx®Opy=ps 
Se 10601 =00 
OMe tat TO" 0. f° 0 
Complement 
A =00011101. n=7 
Diminished Radix Com- 
plement 
A=11100010. @ Check 
J AAAS gee 


= (b"t!—1) Mod (6+1) 
since n+1 is even 
10601 =00 
Radix Complement 
A’=11100011. Check 
pa@Ppa =b"™ Mod (b+1) 
since n+1 is even 


10610=01 
Multiplication RK S=P 
1011 
x 1001 
— Check 
1011 PrObs=Pp 
101100 
10©00=00 


1100011 @ 
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Division X/Y=QO+Vpr/V 


Ore 10011 tox Check 

O10 ey px =(boOphy) Bbre 
ror TN 6s Shel 00 =(01©10) @01 

0 000111 —R ~10@10=00 


The above is a check of 
the first step in the division 
process. 


—sign column 


CONCLUSION 


The congruence notation used provides a fundamental 
basis for the concept of both the digital and the numeri- 
cal type of parity check. Congruence notation also 
provides the mathematical tools needed to consider 
the arithmetic properties of the parity checks and the 
extension of the parity check concept to nonbinary 
systems. 

In the final analysis, the effectiveness of a parity 
check is dependent upon the error structure of the 
carry process. It is necessary to consider in detail the 
mechanism of binary addition before the effectiveness 
of the parity check may be evaluated. A subsequent 
paper by the author considers the error structure of the 
binary adder and the effectiveness of the various parity 
checks. 


Investigations of Magnetic Amplifiers with Feedback« 


HARRY J. GRAY, JR.+ 


Summary—Sine wave carrier excited magnetic amplifiers have 
been investigated to determine if the figure of merit can be im- 
proved through the use of feedback techniques. It is shown that the 
power gain can be made unlimited but that a finite rise time is pre- 
served. Hence the figure of merit as ordinarily defined becomes 
meaningless. It is shown, however, that voltage gain divided by rise 
time remains nearly constant under feedback and is a more suitably 
defined figure of merit. 


INTRODUCTION 


AGNETIC amplifiers are attractive devices 
for computer applications because of their po- 
tentially unlimited life, high reliability, and 


* Manuscript received by the PGEC, March 14, 1958. This work 
was performed while the author was at the Remington Rand Univac 
Div. of Sperry-Rand Corp. and was supported by the AF Cambridge 
Res. Ctr., Air Res. and Dev. Command, under Contract AF19(604)- 
1376. 

+ Moore School of Elec. Eng., University of Pennsylvania, Phila- 
delphia, Pa. Formerly at Remington Rand Univac Div. of Sperry- 
Rand Corp., Philadelphia, Pa. 


relative insensitivity to environment, when compared 
with electron tubes. Carrier magnetic amplifiers have 
not been used in megacycle computer applications in 
the past because of limited power gains and speed. 
However, it has been suggested that feedback tech- 
niques might improve the power gain with little loss of 
speed. 

There are essentially two possibilities: application 
of sufficient positive feedback to produce a regenerative 
pulse amplifier, and use of feedback techniques to pro- 
duce a stable but fast-rise high-gain amplifier having 
many potential applications. This paper deals with the 
problem of understanding the behavior of amplifiers of 
the second type. The first type is a special case of the 
second. 


A Heuristic APPROACH 


Consider the self-saturating magnetic amplifier shown 
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Re 


Fig. 2—Equivalent circuits for self-saturating 
magnetic amplifier. 


in Fig. 1. Pulse tests indicate that the terminal be- 
havior is approximately equivalent to the terminal be- 
havior of the circuit seen in Fig. 2. (Z is resistive.) 

Consider the application of pure negative voltage 
feedback and pure positive current feedback to the 
equivalent circuit of Fig. 2. Practical considerations are 
dealt with later. Then the equivalent circuit becomes 
as shown in Fig. 3. 

The transformed input impedance is 


KpeRLZ 
y Ie or <4 
L 
ZA(s) = Re apy See ee 1 
(s) KZ (1) 
R-+ Rz 
The transformed voltage gain is 
BIR 6 
R-+ Ry 
K(s)°= PotD) 
KZ KpyR1tZ 
R [1 ule Ghee aes 
R-+ Rr R= R; 
The power gain is given by 
ss 
jG pe EN a 
Ry 


with s=0. The power gain is infinite when Z, or K, 
is infinite. Examination of the conditions for which K, 
is infinite shows that Z, is then zero. In the limit, K, 
is infinite. This corresponds to a current actuated am- 
plifier. From (1), Z, is infinite when 

R+ Ry 


KC 3 
j z (3) 


that is, when the current feedback is adjusted to the 
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K¢ilo 


Fig. 3—Equivalent circuits with pure negative voltage feedback 
and pure positive current feedback applied. 


Fig. 4—Magnetic amplifier with positive and negative feedback. 


point of instability in the absence of negative voltage 
feedback. Using this result in (2) yields for the voltage 
gain and 10-90 per cent rise time, K,=1/K,;, and 
PRE RSS 
DO) a O 
Owain 


(> 


The quantity K,/T, becomes 
ZR 


OS ee ee 
2.IE(REF RD 


(4) 
which is a constant of the amplifier. From (2) it is seen 
that (4) is valid also for values of Ky; other than given 
by -@); 

This last case is similar to vacuum-tube operation. 
Hence, at present, it is the most interesting case. 

‘The following amplifier was tested in the circuit of 
of Fig. 4. 

Cores =4 wraps of $ mil by 1/32 inch 4-79 

Mo-Permalloy on 0.t-inch-diameter form. 

Output winding =111 T No. 46 HF. 

Control winding=111 T No. 46 HF (links both 

cores). 

Diodes = IN64. 

Re == 240: 

R,= 100. 

RF supply = 12-volt peak at 2.25 mc. 

Some early data are summarized in Table I for 0.1- 
volt pulse input. The voltage gain for Ry=470Q (ap- 
proximate point of infinite input impedance) was 25 
and 5.8 for Ky, equal to 0.02 and 0.2, respectively. The 
agreement with theory is sufficiently good to lend weight 
to the argument. 
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TABLE I 
K./T.(us7) 
Ry 
470 | 560 | 680 | 820 | 1000 | 1200 | 1500 
Kyy 
O02 \GalgMles Odie? | sty eleaol 1s \t 1.6 
0.2 | MoM 10.1) “a7 M245 


The implications of (4) are that the rise time has a 
limit of zero as K, approaches zero. This evidently is 
not true with a finite RF supply frequency. Further- 
more the analysis does not give the insight into the 
performance of the amplifier that is necessary for de- 
sign. A more rigorous approach requires a cycle-by- 
cycle analysis in the manner of Johannessen! whose re- 
markable paper suggested some of the items discussed 
here. 


DIFFERENCE EQUATION ANALYSIS 


The reader is referred to the Johannessen paper for 
the mechanics of the analysis. The same assumptions 
are made here as in that paper, except that: 


1) The reactors are assumed to have a finite mag- 
netizing impedance Z,, which is essentially resistive 
in the megacycle range. 

2) The winding resistance is negligible. However, 
a loss of output volt-seconds occurs due to diode 
drops and winding reactance when the core is 
saturated. The effect is assumed to be equivalent 
to the loss of volt-seconds caused by a fictitious 
winding resistance, Ry. 


No External Feedback 
A nodal analysis of the circuit of Fig. 1 yields 


Eo(n + 1) = KyoE.(n) + KsEo(n) (5) 
where 
Ry, Ls Lila 
Kyu = N : = 
Rea Ry NAR Z Ze Rs 
A Zm\|N?R. 2Rr+ Rw 
Ky = =a : (6) 
Zm + N?R||Re Rz+ Znl|[N?R. Ri + Rw 


The solution of the difference equation (5) for a step 
in e, of amplitude £, applied at the start of interval zero 
is 


Mae 


E,(n) ae hae Sie (7) 


sy 

Avoiding the quagmire of precise definition, one ob- 
tains from (7) as a measure of the rise time and maxi- 
mum voltage gain 


1 P. R. Johannessen, “Analysis of magnetic amplifiers by the use 
of difference equations,” Trans. AIEE (Commun. and Electronics), 
vol. 73, pp. 700-711; January, 1955. 
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T BD ni 
etapa 
) In Uy a) 
eae (8) 
Get tose RY 
The ratio K,/T, becomes 
K, WK yo 
~ (9) 


A a2 
m(1 — Ky) (1 + anoere: - 
which is not independent of Ky as in the above section. 
However, the quantity 
Kay or 
Tisai 


is plotted vs K; and in Fig. 5, demonstrating that this 
quantity is fairly constant over a wide range. 
Hence for K;>0, a good approximation is 


ISG wK yo 


— = Ree 
Ee Qa f 


(10) 


Thus it is reasonable to take fK,, as a figure of merit. 


Properties of Kyo 


Eq. (5) describes the behavior of an amplifier which 
possesses a single time constant; hence it is the simplest 
to analyze. Internal feedback is present as shown by the 
term K,E)(m). All amplifiers discussed in this paper 
can be described by (5). External feedback will be used 
to control Ky and modify the input impedance. _ 

Eq. (10) is a direct and fundamental consequence of 
(5) giving the figure of merit, fK,., some general signifi- 
cance. From (6) the reciprocal of K,, becomes 


1 1 (14 ies! ssa 
Ro) RS | ees si 


a linear function of R, or Ry. This result allows an exper- 
imental check to be made. The amplifier described in 
the preceding section was operated with a 10-volt peak 
RF supply at 2.25 mc and with a 100—ohm load. Meas- 
urements of voltage gain and rise time were made 
and K,, computed from (8). The results are plotted vs 
R, in Fig. 6. The amplifier was adjusted in each case 
for maximum voltage gain with a peak filtered output 
of 4 volt corresponding to a signal input utilizing fully 
the dynamic range of the amplifier. 

The linear dependence on R, is evident and the slope 
is consistent with theory if a finite diode reverse impe- 
dance (R,) is assumed. However, characteristically of 
all data of this type, the intercept is too low. The proper 
intercept is obtained when the amplifier is adjusted 
for maximum voltage gain with a small-signal input. 
Plots of this type were made in order to determine 
experimentally the maximum value of K,.(Kyom) for an 
amplifier to compare various amplifier structures before 
feedback is applied. 
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02 0.3 0.4 0.5 0.6 07 08 09 1.0 
K¢ 


Fig. 5—(K,/T;,) (4/wK vo) vs Ky. 


THEORY (Zm=!K, Rw=252., Ry= © 


EXPERIMENT 
0.5 
O.; 
{0} | a ee | It J 1 : 
10) 100 200 ~300 400 500 600 £700 800 900 1000 


Rc OHMS 


Fig. 6—1/K 0 vs Re. 


From (6) Kyom=N, and from (9) the maximum value 
of K,/T, is Nw/x. This result can be derived more 
directly as follows. The input volt-seconds are trans- 
formed by the control to output turns ratio, V. Hence, 
neglecting losses, the maximum voltage gain is V. The 
output occurs in one-half cycle, t/w. Taking this to be 
the rise time, the maximum figure of merit Nw/m is 
yielded. In a like manner, a single-core amplifier has a 
maximum voltage gain of JN, a rise time of one cycle, 
and therefore a figure of merit, Nw/27. 


External Feedback 


Consider the circuit shown in Fig. 7. A nodal analysis 
yields for K,, and Ky, 


‘Ry Mi, 
her = N ————— 
R,+ R, N?R. 
ZL Ry, 
Ge - 
N SG Bose UKeN 
NARI|Ril| = Ry 
M 
Z 
(11) 
NA Re N Ry, 
= 
NK fo M 23 aI Ray 
where 
N 2 
Z = N?R|Zn| Rall (=) R;. 
M 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


September 


Fig. 7—Magnetic amplifier with feedback winding. 


The input impedance is, statically, 
R, 
i KS, 


¢ (12) 
Eq. (12) is sufficiently startling to require an outline of 
its derivation. Fig. 8 shows equivalent circuits for 
Fig. 7 for two cases. 

Solving for 


: Nee + €r = 
LIN SS ae cores unsaturated 
NoR, 
. Nec + er — NKyo 
lip = core II saturated. 
N?R, 


Averaging the currents over a cycle yields (12). 

Note that the input impedance is infinite when 
Ky =1/K,. To determine the necessary circuit condi- 
tions, set Ky,=0 in the second part of (11) and set 
Ky;=1. The resulting equation 


ZL, 


N?R, 


N 
| Rell = Rs 


Ry Z 
East ieee RELY, 


ears 


= =1 (13) 


yea 


z 


is a constraint on the circuit constants. Calculation of 
1— Ky; subject to the constraint of (13) yields 


$ R, NKj 
{Ke Pr (14) 
Ree eve: 
and 
Ks 1 
fs ded one keaear Ks 


follows from the first part of (11). Hence, when (13) is 
satisfied, the input impedance is infinite and K,=1/Ky. 
The following amplifier was constructed. 
Cores =3 wraps of $ mil X1/32 inch 4-79 Mo-Permal- 
loy on 0.1-inch-diameter form. 
Power winding =147T No. 46 HF. 
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CORES UNSATURATED 


CORE Il SATURATED 


Fig. 8—Equivalent circuits for magnetic amplifier 
with feedback winding. 


Control and feedback windings = 


2—147T No. 46 HF 
1—74T No. 46 HF 
3) 1 No. 46 BE 
1—18T No. 46 HF. 


Diodes=EM15 (high conductance, fast recovery 

type). 

Carrier = 10-volt peak at 2.25 mc. 

Rw =50Q by power output vs R,z test. 

Various combinations of turns for the feedback and 
control windings were tried along with various values of 
Ky. Voltage gains, power gains, and rise times were 
measured for a step input; rise times were estimated by 
cycle counting. 

Typical data with R;=2.2K, N=2, M=4, Rr, =1000 
and Z,~ «© areshown in Table II. 


TABLE II 
R. Kyo Ks T; (usec) Ge il KG 
100 0.5 2 1.34 1.49 0.66 
470 0.5 2 2.9 0.69 0.31 
1K 0.5 2, 4.25 0.47 On2it 
1K 0.7 LAS 3.8 0.38 0.17 
1K 1.0 0 2.45 | 0.41 0.18 


The relation K,=1/Ky, is seen to be satisfied by the 
data. Note that if R, is fixed, K,/T;, is very nearly con- 
stant. Also note the dependence of K,, on R,. 

That K,, decreases when R, increases is a severe 
practical handicap because a driving source having 
little power output will have a high internal impedance 
contributing to R, and lowering K,., hence reducing 
the voltage gain bandwidth. Empiricism has shown that 

the addition of frequency sensitive circuit elements in 
the output and feedback circuits limits the reduction 
of K,,. as R, is increased. Such measures when applied 
to this amplifier yielded a voltage gain bandwidth as 
defined by (10) of 1.8 mc. which was independent of 
the turns ratios for 1.1<M<2.7 and surprisingly for 
0.7<N<1.6. The value of R, in these tests was 4.7K. 
The dependence on R, was slight. 


CONCLUSIONS 


It has been demonstrated that feedback techniques 
can be applied successfully to carrier magnetic ampli- 
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fiers to obtain useful results. It is shown that power 
gain bandwidth is not as meaningful as a figure of merit 
as is voltage gain bandwidth. Furthermore, an amplifier 
subjected to combined positive and negative feedback 
can be adjusted to have an infinite input impedance. 
Under these circumstances, its voltage gain is equal to 
the reciprocal of the voltage feedback factor and its rise 
time is given by K,/fKyo. 


DEFINITIONS OF SYMBOLS 


é-=instantaneous control voltage. 
és =instantaneous supply voltage. 
é) =instantaneous output voltage. 
é€r=instantaneous voltage across core I. 
E=average voltages, averaged over a half- 
cycle of the supply voltage. The same 
subscripts as for instantaneous voltages 
are used. 
f=frequency of supply voltage, cycles per 
second w=2rf. 
1=instantaneous current. 
1, =instantaneous control circuit current. 
Ky; =feedback factor. 
Ky, =voltage feedback factor. 
Ky;=current feedback factor. 
K,=voltage gain. 
K,=power gain. 
K,yo= voltage gain before feedback. 
Kyom=Maximum value of Kyo. 
M=turns ratio of feedback winding to con- 
trol winding. 
N=turns ratio of output winding to control 
winding. 
R,=resistance of control circuit. 
R;=feedback resistor. 
Rr=resistance of load. 
Ry =equivalent resistance of output winding. 
R,=equivalent leakage resistance of diode. 
s=complex frequency variable. 
T,,=rise time of output. 
Z-=input impedance. 
Zm= magnetizing impedance of reactor re- 
ferred to output winding. 
Z =an impedance. 
n—1,n,n+1=three consecutive half-cycles of the sup- 
ply voltage. 
| =an operator enabling impedances instead 
of admittances to be written in a nodal 
analysis. Note that 


Ri Rs 1 
ee ee 
Ry + Rs G, aia G2 
Also, 
Ri|| — Re = — (—Ril| Re) 
and 


R,||R.|| — Re = Ra. 
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A New Class of Digital Division Methods* 


JAMES E. ROBERTSON{ 


Summary—This paper describes a class of division methods best 
suited for use in digital computers with facilities for floating point 
arithmetic. The division methods may be contrasted with conven- 
tional division procedures by considering the nature of each quotient 
digit as generated during the division process. In restoring division, 
each quotient digit has one of the values 0, 1,-- - , r—1, for an ar- 
bitrary integer radix r. In nonrestoring division, each quotient digit 
has one of the values —(r—1),---, —1, +1,---, +(r—1). For 
the division methods described here, each quotient digit has one of 
the values —n, —(n—1),---,-—1,0,1,---,n—1, n, where nis an 
integer such that #(r—1)<n<r—1. A method for serial conversion 
of the quotient digits to conventional (restoring) form is given. Exam- 
ples of new division procedures for radix 4 and radix 10 are given. 


INTRODUCTION! 
A DIVISION method can be categorized by listing 


the permissible values of each quotient digit as 

generated during the division process. For an 
arbitrary radix 7, each quotient digit generated during a 
conventional restoring division has one of the values 
0,1, --+,7—1. For nonrestoring division each quotient 
digit has one of the values —(r—1), —(r—2), ---, —1, 
1, 2,---,7—1, with 0 excluded. The purpose of this 
paper is to describe a class of division methods in which 
each quotient digit has one of the values —n, —(n 
—1),---, —1, 0, 1,---, , where m is an integer 
such that 4(r—1) <u <r—1. Conversion of the quotient 
to the conventional restoring form is required for the 
latter two classes of division methods. 

Division, as executed in most digital computers now 
in use, involves a recursive process which may be pre- 
ceded by preliminary operations and followed by termi- 
nal operations. Most of the time required for a digital 
division is spent in the repeated execution of the recur- 
sive process. For nonrestoring division and for the class 
of division methods proposed here, the recursive process 
can be described by 


Lip. = 1Hj — Qitid pi nO lenceria A 


for which the following notation is employed: 


x;=partial remainder resulting from the jth execu- 
tion of the recursive process 
x9 = dividend 


* Manuscript received by the PGEC, March 19, 1958; revised 
manuscript received, July 1, 1958. The investigations leading to this 
paper were part of a computer study program by the staff of the 
Digital Computer Lab., University of Illinois, Urbana, Ill., under the 
joint support of the Office of Naval Res. and the Atomic Energy 

omm. 

+ University of Illinois, Urbana, III. 

1 Descriptions of conventional restoring and nonrestoring meth- 
ods can be found in a number of texts, such as the following. 

R. K. Richards, “Arithmetic Operations in Digital Computers,” 
D. Van Nostrand Co., Inc., New York, N. Y.; 1955. 

M. Phister, Jr., “Logical Design of Digital Computers,” John 
Wiley and Sons, Inc., New York, N. Y.; 1958. 


Xm = remainder 
q; = (for fractions) the jth digit of the quotient to the 
right of the radix point 
m =the number of digits, radix 7, used to represent 
the quotient 
d= divisor. 


Each x;(j=0, 1,---, m) satisfies |x,|<k|d|, and 
the sign of each gji1 is chosen so that | x54: = | (r| x;| 
—| gj41| |d|)]. A straightforward but lengthy analysis 
of all cases that may arise reveals another property of 
each quotient digit, namely | gj <k(r—1). For non- 
restoring division, =1; it is shown that division meth- 
ods exist for certain discrete values of & in the range 
L<R<1, 

From the equation for the recursive process it can 
be shown readily that the division procedure is correct. 

For j=0, 


x1 = rx — quid. 
For7 = 
Xo = 1X1 — god = r°xX0 — (rqi + qa)d. 
For j=m—1, 


Dig ea gr ohne 


> + 1GQm—1 + Gm)d. 


The shifted remainder 7~"x,, is then 


m 
1 "tm, = %9 — d >, rg; 


t=1 


where 
™m 
ae 
4=1 


represents the quotient Q. It follows that Qd+r-™x,, 
=X 9; 1.e., the sum of the shifted remainder and the 
product of the quotient and divisor is the dividend. 

The mechanization of the recursive process requires 
three distinct steps. 


1) The partial remainder x; is shifted, 7.e., multiplied 
by the radix r. 

2) One of several permissible arithmetic procedures is 
selected, such that the maximum absolute value 
of r| x; | , hamely kr| d| , is reduced by the amount 
1/r, so that the result x41 satisfies | x41 <k\d]. 
It should be emphasized that the reduction is in 
the range over which partial remainders may vary, 
and not necessarily in the absolute values of 
specific partial remainders. 

3) A quotient digit is generated corresponding to the 
arithmetic procedure selected. 


PO Fe a ee ee eee eS 
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The key to the study of division methods lies in the 
analysis of arithmetic procedures which reduce the al- 
lowable range of absolute values of the shifted partial 
remainder (7x;) by the amount 1/r. 


ANALYSIS OF ARITHMETIC PROCEDURES 


It is convenient to normalize the partial remainders 
with respect to the absolute value of the divisor. If the 
substitution z;=x,/|d| is made, the range restrictions 
become —k<2;41<k and —rk<rz;<rk. Attention then 
is focused on arithmetic procedures which transform 
rz; into z;41. The easily mechanized procedures involve 
addition or subtraction of integral multiples of the 
divisor from rx; to yield x,4:; after normalization, the 
procedures involve addition or subtraction of integers 
from rz; to yield 2;,:. The arithmetic procedures can 
be represented as a family of straight lines of the form 
8j41=72;—1, where a=—n,--+-, —1, 0, 1, 2,---, 2, 
as shown in Fig. 1. 


cn A 
em! 


Zep att! — 24,22 I 


ii 27417 -F 


2)41 =(Zj—A 


Fig. 1—A normalized graph illustrating arithmetic 
procedures during division. 


In order for one of the proposed class of division 
methods to exist, it must be possible to superimpose a 
rectangle on the family of straight lines of Fig. 1 in 
such a way that the following occur. 


1) The rectangle is centered at the origin, with 
vertices at (+rk, +k). 

2) The projections on the rz; axis of the line segments 
within the rectangle cover that portion of the rz; 
axis within the rectangle. 


It follows from the first condition that the vertices 
of the rectangle lie on the lines through the origin of 
slope +1/r, whose equations are 2;41= +2;. The second 
condition requires that the rectangle be sufficiently 
large to insure that k >}. Two additional considerations 
govern the choice of size of the rectangle, 7.¢., the choice 
of k and of n. 


3) The number of lines of the form 2;41=72;—7 neces- 
sary to satisfy condition 2) should be minimized. 
The number of multiples of the divisor which must 
be formed is proportional to the number of lines 
employed. 

4) The overlap of projections of the line segments on 
the rz; axis should be maximized. The precision 
necessary in the selection process decreases as the 
overlap increases. 
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The two considerations are mutually contradictory 
since 3) requires that the size of the rectangle should be 
decreased and 4) requires that the size be increased. 

The considerations govern the choice of & to the 
extent that k should take one of a discrete set of values 
such that the vertex (rk, k) of the rectangle lies on a line 
2341 =12;—n. Any one division method can then be char- 
acterized by the positive integers chosen for r and for n. 

The value of k as a function of y and 1 can be found 
by solving for the point of intersection of the lines 
2j;41=2; and 2;41=72;—n. The value of 2,41 at the point 
of intersection is k, and is found to be u/(r—1). The 
requirement that k >} becomes n> (r—1)/2. 

The choice of » for some given radix r involves a 
balance between time and equipment costs associated 
with the selection process, on the one hand, and similar 
costs in forming multiples of the divisor, on the other. 
Since the balance is so much a function of design details, 
the choice of ” is discussed further only in connection 
with specific examples. 


QUOTIENT CONVERSION 


The conventional representation of a quotient re- 
quires that each digit be one of the positive integers 
0, 1,--+, r—1. Since nonrestoring division and the 
method described here involve negative digits in the 
quotient, some means of conversion is required. The 
technique employed in conventional nonrestoring di- 
vision, except for the special case of the binary system, 
can be described for a radix complement representation 
by the following rules. 


1) If q.<0, replace qi by qi.’ =r+qi and set the sign 
of the quotient negative; if q>0, gq’ =q, and the 
quotient is positive. 

2) For j=1, 2," , m=1, inspeet’ g;” and "Gyan. 
If qji1<0, replace q;’ by gq;’—1 and replace 
Git by gin’ =r+gin. UH giz >0, 9,’ is left un- 
changed, and qj41’ =qj41. 


The rules require a serial inspection of the quotient 
digits, the most significant digit first. When a negative 
digit is encountered, it is added to the radix, and a unit 
is borrowed from the next most significant digit. 

For the proposed division method, the conversion is 
complicated by the fact that 0’s are permissible quotient 
digits. The inspection of the sign of g;41 (or q in rule 1) 
must be replaced by an inspection of signs of the divisor 
d and the partial remainder x; (xo in rule 1) to determine 
the sign of the next zonzero quotient digit. Agreement of 
signs of x; and d corresponds to q;41>0 in the above 
rules; disagreement corresponds to qj41<0. Alterna- 
tively, signs can be determined in the usual way and 
can be associated with those q;1:1 which are zero. These 
modifications provide for a borrow propagation through 
a sequence of 0’s. 

Both sets of conversion rules require that no quotient 
digit g; be such that | g;| >r—1. In particular, the 
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parameter 7 is one value that |; | can assume, and 
therefore »<r—1. Thus, 7 is restricted to the range 
4(r—1)<n<r-—1, since it was established previously 
that 4235 (F—1). 


PRELIMINARY AND TERMINAL OPERATIONS 


Some of the requirements which necessitate prelimi- 
nary or terminal operations for division methods are 
listed below. 


1) The requirement of standardizing the quotient, in 
a floating point division, and the requirement of 
overflow detection in a fixed point unit. 

2) The requirement for a rounded quotient. 

3) The requirement that a correct remainder be 
generated. 


Procedures vary in computers now in use. For the di- 
vision methods described in this paper, the above re- 
quirements necessitate comparable preliminary or termi- 
nal operations. 

An additional requirement imposed by the division 
methods discussed is that each partial remainder x; 
should satisfy | x,| <k|d| where d is the divisor and 
4<k<1. In particular, the restriction applies to the 
dividend x», and may necessitate additional preliminary 
operations for division methods for which k<1, in 
contrast to conventional procedures for which k=1. 

Preliminary standardization of the divisor simplifies 
the selection of the correct multiple of the divisor during 
the recursive process, since the precision required for 
selection increases as the minimum absolute value of the 
divisor decreases. The proposed methods, therefore, 
are best suited for use in arithmetic units having facili- 
ties for floating point operations. The methods can be 
used for fixed point division if facilities for simultane- 
ously shifting divisor and dividend are available. 

Additional problems are posed when a remainder 
Xm Must be generated such that the division algorithm 
Od+r-"xm=xo (where Q is the quotient represented 
by m digits, radix 7; dis the divisor, and x» the dividend), 
is satisfied. If d and xp are initially shifted left p digital 
positions so that d’=r?d and x0’ =r?xo, then the value 
OL xy “such; thatsOd)--r-"wn, =Xo 1s foundstos be i,’ 
=r?x,. For a fixed point division, it would thus be 
necessary to shift the remainder x,,’ p digital positions 
right to obtain the correct remainder x,,. 

A second difficulty arises when the quotient conver- 
sion rules require the least significant digit gm’ of the 
converted quotient to have the value 7. If, as is often 
the case, facilities for addition are not available for the 
quotient register, g,,’ can be set to the value r—1, and 
the remainder x, must then be replaced by «n’=xn-+d. 
The division algorithm becomes 


(0 — r-™)d + r-™(am + d) = Xo. 


Similar results are obtained when conventional non- 
restoring division methods are employed. 
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Example 1—Radix 4 Division 


Conventional radix 4 division methods require either 
two uses of a binary adder (use of one adder sequentially 
or two adders in parallel) or the formation and storage 
to full precision of 3d, where d is the divisor. The di- 
vision method of the class proposed here with r=4 and 
n= 2 requires a single binary adder, conditional doubling 
and complementing circuits, and a selection circuit to 
compare 7x; with 0.5d and 1.5d to a precision of 7 
binary digits, if }<|d| <1. 

The mechanization of the division scheme is indicated 
diagrammatically in Fig. 2. 


BINARY _ 
ADDER 
COMPLEMENTING 


—>— 


CIRCUIT 


CONDITIONAL 
—=— DOUBLING 
CIRCUIT 


Fig. 2—Radix 4 division. 


The selection circuit performs three functions, based 
upon the seven most significant binary digits of 4x; (or 
x; if desired) and of d. 


1) Compares signs of x; and d. If signs agree, the 
complementing circuit is set to subtract; if signs 
disagree, the complementing circuit is set to add. 

2) Compares absolute values of 4x; and 1.5d. If 
4| x;| > 13] d| , the conditional doubling circuit 
must be set to form 2d. If 4|x,;| <14|d|, the cir- 
cuit must be set to form d. If 14|d| <4|.;| <12|d|, 
the conditional doubling circuit can be set either 
way, depending upon the design details of the 
selection circuit. 

_ 3) Makes a similar comparison of absolute values of 
4x; and ($+4%)d. For the smaller values of 4| x, | ; 
Xj41=4x;; otherwise x;41 is transferred from the 
binary adder. 


After x;41 is formed, 4x,4: is formed by a radix 4 left 
shift to replace 4«;. The values selected for the quotient 
digit g;;1 must correspond to the selections made, as sum- 
marized in Table I. 


TABLE I 
Selection 1 Selection 2 Selection 3 din 
Subtract 2d 934140 +2 
Subtract 1d Gi =0 +0 
Subtract ld 934140 +1 
Add 2d 0 —2 
Add 1d = () 


Add 1d ~0 = 
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Example 2—Radix 10 Division 


From the many possible choices available to the de- 
signer, the radix 10 division method of this example is 
based upon the following. 


1) The excess three representation of decimal digits 


is chosen. 

2) A single decimal adder (excess three code) is used 
sequentially. 

3) Storage is provided for the divisor d, but not for 
any of its multiples. 

4) A complementing circuit and a conditional dou- 
bling and quintupling circuit are employed, with 
+d, +2d, and +5d available as inputs to the 
adder. 

5) The division method is characterized by r=10, 
w= 72 with 

n fi 
nee ae 
oak. 9 


Doubling and quintupling circuits can be described 
by the sets of Boolean equations 


q = b(a\/ c\V de) V/ ac(d\ e) 
r = é ® (ade \/ adé) 
s=a@OdOe for doubling 
t=€é 
b= 
and 
v= 
w= alé\/b\V cd) V éb(c V d) 
x = b @ (ecd \/ ééd) for quintupling 
y=20cOd 
zZ=d 


where, in excess three notation, a, b, c, d; q, 7, s, t; and 
v, W, xX, y represent one decimal digit of the divisor, 
2x(divisor), and 5x(divisor), respectively. 

For doubling, e is the binary carry input and u is the 
binary carry output to the next most significant decimal 
digit. For quintupling, e and z are most easily described 
as incoming and outgoing binary borrow signals in a 
halving circuit, with a wired-in decimal shift of v, w, x, 
and y converting the halving circuit to a quintupling 
circuit. 

If the permutation a'=)b, b’=c, c’=d, d'’=e, e’=4 
is made in the divisor digits, and the permutation 
v=4, w=q, X=, y=S, =, at the output of the circuit, 
the doubling circuit is transformed into a quintupling 
circuit. Thus the same hardware can be used sequen- 
tially for both doubling and quintupling, provided 
permutation and complementation of input and output 
signals is correctly arranged. 

The arrangement of hardware required for the di- 
vision scheme is shown in Fig. 3. The two steps required 
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for the generation of each decimal quotient digit can be 
described as follows. 

Each quotient digit g;;1 can be decomposed into two 
digits Gia; Ghee 4, such that G41 = Qi’ +qj41", where 
Qin1 = —5, 0, or 5 and gj41/’= —2, —1, 0, 1, or 2. Step 
1 corresponds to the determination of q;41' and results 
in the formation of a quantity x,41’ such that 


Kp3 = 102; if Qi+1 = (i) 


or 


DECIMAL 
ADDER 
EXCESS 3 


COMPLEMENTING 
CIRCUIT 


PERMUTATION 
CIRCUITS 


[ SIGN COMPARISON | 
(1) COMPARE 10 Xp, 25d 
(2q COMPARE Xjt1, 15d 
(2b) COMPARE Xj+ 1, O.5d 


QUINTUPLING - 
CIRCUIT 


Fig. 3—A decimal division method. 


Step 1 reduces the range 0<10|x;| <7$3|d| to 
O<|xj41'| <24\d . The result of step 2 is %j41, whose 
range is 0< | Xj <$| d|. Step 1 is not required when- 
ever g;41' =0. The details of operations performed for all 
ranges of 10x; and xj;11’ are summarized in Table II. 
The overlap of $|d| in the ranges is such that three 
decimal digits of d and of 10x; or x;,1’ are required for 


TABLE II 
Range of 10x; Step 1 Ga Range of x41’ 
—7£|d| <10xj<—22|d| Add 5|a| 28 
—22|d|<10x;< 2%|d| Not required 0}—24|d| <xj41'<2$|d| 
22|d|<10x;< 7$\d Subtract 5|d| +5) 
Range of a;44’ Step 2 G+" Range of aj41 
~24|d| <ajy"<—13|d| Add 2/4 9 
12 |d| <xj’<— B|d| Add |d| ay 
— $ld|<xjp'< Fld 0}—-4|d| <x <$ld| 
2\d| <ajy'< 19ld Subtract |d| +1 
12|d| <xju'< 25 |d Subtract 2|d| +2 
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comparison, provided d is standardized to the range 
1/10<|d| <1. 

The average number of operations necessary for the 
division method, assuming all quotient digits —7,---, 
+7 are equally likely, is 14 operations per digit of the 
quotient. This figure is to be compared with 3.4 opera- 
tions per quotient digit for a conventional nonrestoring 
division method employing doubling and quintupling 
circuits.? 


CONCLUSION 


The new methods of division described here, when 
coupled with two rather obvious comments, show 
promise of leading to new developments in digital 
arithmetic. These comments are: 


1) Multiplication is the inverse of division. 
2) The representation of quotient digits is, at least 
in the interesting cases, redundant. 


This paper is the result of an attempt to find a division 
inverse to a multiplication method recently described 
by Lehman.’ The radix 4 division example is the inverse 
of the radix 4 equivalent of the binary multiplication 
described by Lehman, except for the manner in which 
redundancy is employed in the generation of quoteint 
digits in the one case and the recoding of multiplier 
digits in the other. The multiplication method which 
is the inverse to the radix 10 division is similar to that 
used on the IBM 602A computing punch; the fact that 
fewer operations are required per quotient digit for the 
division method indicates that a better multiplier digit 
encoding scheme would reduce the number of opera- 
tions required for multiplication. It is clear that multi- 
plication methods inverse to other division methods of 
the class described here must exist, and it is felt that 
further investigation will lead to a broader understand- 
ing of digital arithmetic. The interrelationships be- 
tween quotient digit redundancy, selection procedures, 
and divisor multiple generation should be studied also, 
particularly in those cases for which several steps per 
quotient digit are desirable (e.g., the radix 10 example). 


2 Richards, op. cit., pp. 274-275. 

3M. Lehman, “High-speed digital multiplication,” IRE TRANs. 
on ELectronic Computers, vol. EC-6, pp. 204-205; September, 
1957. 

4 Richards, op. cit., pp. 262-263. 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


September 


APPENDIX 
PRECISION REQUIRED FOR SELECTION 


An estimate of the precision required of the divisor d 
and the shifted partial remainder rx; for selection of the 
correct arithmetic procedure can be gained by estimat- 
ing the variation in the ratio rx,;/d resulting from trunca- 
tion of rx; and of d. This variation must be less than 
the overlap of projections on the rz; axis (Fig. 1) of two 
successive lines of the form 2;41=72;—-7@ 2;41' =r2;—(4 
+1). 

The overlap is the difference in values of rz; for 
Zj41=n/(r—1), and for gj’ = —n/(r—1), and is 2n/(r 
—1)—1. Truncation errors in rx; and d can be expressed 
by setting a<|x,| <a+Aa, 6<|d| <b+Abd. The varia- 
tion in the estimates of the ratio | rx;/d| is then 


[== a ] 
her (An 


It is required therefore that 
a+ Aa a 
b b+ Ab 


1 a 1 2n 
~—(a0+ as) <a —( — 1). 
b b r\r—-—1 


Assuming Aa =Ad, one obtains 


b[2n — (r — 1)] 
(ry — 1)r(1 + a/d) 
The minimum value of 6, assuming standardization, 


radix 7, is b=1/r; the maximum value of r(a/b) is 
4(2n—1). Therefore, 


z 2[2n — (r — 1)] 
r(r — 1)(2r + 2n — 1) 


For the example r=4, n=2; Aa<1/66, indicating 
that seven binary digits are required for the comparison. 
For r=10, n=7; Aa <1/297, indicating that three deci- 
mal digits are sufficient for the comparison. 
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Magnetic Core Pulse-Switching Circuits for 
Standard Packages” 


JACK L. ROSENFELD} 


Summary—A new method for the logical design of magnetic core 
pulse-switching circuits is presented. This method has features 
which make it excellent for use in standard packages. These features 
are the absence of spurious noise signals at the output; the fact that 
outputs are independent of the order of arrival of input pulses; the 
fact that interchanging components does not affect circuit behavior; 
and the fact that moderate changes in clock pulse amplitude and 
duration do not cause false operation. Furthermore, the number of 
components required by this system compares favorably with the 
numbers required by other systems. A computing system can be built 
by properly interconnecting a few different types of such packages. 

The new system uses advance current drive but performs the 
logical operations with both forward and backward windings in an 
output network. The use of both types of windings permits a reduc- 
tion in the total number of cores required. 

Tests performed on actual circuits yielded very encouraging re- 
sults. The circuits operated as predicted, and the performance was 
most satisfactory. 


INTRODUCTION 
Se methods [1 ]—[8] have been developed for 


using magnetic cores as logical elements in com- 
binational switching circuits. For certain applica- 
tions an all-core computing system (except for transistor 
or vacuum tube driver) may be desired, and to design 
each individual switching circuit would be a prodigious 
task for even a small computing system. No present 
method of constructing magnetic core switching circuits 
seems to lend itself to economical standardization. 
This paper describes a method of magnetic core pulse- 
switching which enables the logical designer to con- 
struct a switching circuit by wiring together a few 
standard building blocks, using a small number of 
cores and diodes. The method is quite economical and 
especially adaptable for application to standard build- 
ing blocks. 


ANALYSIS OF LOGICAL DESIGN 


In this paper magnetic core circuits are drawn with 
modified mirror symbols [5]. The fact that a conductor 
is wound on a particular core is indicated by a short 
line segment (the “mirror”) cutting the conductor 
line at a 45° angle (see Fig. 1). This mirror is labeled 
with the designation of the core on which the winding 
appears. The number of turns in a winding may be 
noted beside the mirror symbol (N., Ms, and JN, in 
Fig. 1). Current through the winding on a core produces 


* Manuscript received by the PGEC, March 20, 1958. This paper 
is based on the author’s thesis of the same title submitted in partial 
fulfillment of the requirements for the M.S. degree at M.I.T., Dept. 
of Elec. Eng., June, 1957. The work described was done while the 
author was employed by Bell Telephone Labs. 

j M.I.T., Cambridge, Mass. 


flux in either of two directions, depending upon the 
direction of the winding. This direction is specified in 
mirror symbols by “reflecting” the current in the mirror. 
When the direction of flux in a core is changed by cur- 
rent in one of its windings, a voltage appears on the 
other windings. The polarity of this voltage can be 
found by reflecting the setting current in its mirror, 
inverting the direction of the reflection, and reflecting 
the resultant “ray” in the mirrors for each winding of 
the core. These secondary reflections show the direction 
of the induced voltages. Fig. 1 shows magnetic core A 
set in the “downward” direction by current J. If the 
flux was previously “upward,” voltage is induced in 
the other windings, tending to force current to 
flow “backward” (from right to left) in lead 4 and 
“forward” (from left to right) in lead l,. Lead J3 is not 
connected to a winding on core A. Unless otherwise 
indicated, input currents flow from left to right. 


Windings 


Fig. 1—Modified mirror symbols. 


A randomly selected switching function of four 
variables, F(x1, 2, %3, X4) =%3/X4 +1! x2%3/%4! +04! x0! x3%4 
+2x1%2%3%4', was chosen to illustrate the operation of the 
new logical system. This is shown in Fig. 2. The pulse 
sequence is as follows: inputs x; are pulsed; then ad- 
vance current A is pulsed, producing output F across 
load Zz. When the advance lead is pulsed, those cores 
to which inputs have been previously applied deliver 
voltages across their output windings. Each of the four 
leads 1; corresponds to one of the four terms in a dis- 
junctive expression of the demonstration function. 
Observe that one and only one winding on each lead is 
wound in the forward direction. The resultant voltage 
across a lead (from node a to node 6 for lead J;) can be 
in the forward direction only if the one forward winding 
on that lead, and no other winding, is excited. For 
example, lead /, will have forward current induced if 
core x3’ has been set and core x4’ has not been set. If 


N 


(a) (b) 


Fig. 2—New logical system. (a) Inputs, (b) output windings. 


both cores have been set, their output voltages cancel 
each other. If core x,’ alone has been set, the back volt- 
age is such as to prevent current from flowing forward 
in lead /,. Thus this lead corresponds to the term x3/x4. 
When this lead, or any other lead, has a forward current 
induced, the advance current is forced to pass through 
load Zz. When none of the leads are forward biased, 
the advance current flows through diode d;. Thus the 
output function is realized. 

These output windings may be combined in series- 
parallel fashion to eliminate windings. For example, the 
dashed wire may be soldered in and the dashed cut made 
in lead J3. This eliminates one of the windings on core 
x, and corresponds to the identity 21'x203/x4' +-x1' x0! x3%4 
= x1 (x0%x3'x4' +x2'x3%4). Similar modifications may be 
made as long as every path has one and only one for- 
ward winding, as illustrated in Fig. 3. (Manipulations 
for removing superfluous diodes are discussed in the 
following section.) 

This new system is partly a combination of systems 
developed by others, combining them in such a manner 
as to produce decided advantages. The novel feature is 
the introduction of both forward and backward wind- 
ings in the output network. The forward windings force 
current through the load for F=1; in this respect the 
circuit behaves like an AF switch [5]. When F=0, the 
backward windings prevent current from flowing to the 
load, as in AB-type circuits. The advantages of this 
system will become clear as the discussion proceeds. 


RESTRICTIONS ON LOGICAL DESIGN 


Some of the difficulties encountered in the electrical 
design of magnetic core switching circuits can most 
easily be eliminated by placing restrictions on the 
logical design of the circuitry. In the following discus- 
sion it will be assumed that the loads of the core switch- 
ing circuits are other cores—an appropriate decision for 
an all-core computing system. But first, some general 
considerations should be made. 
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Fig. 3—Output network in minimal form. 


The load core of a switching circuit must switch 
faster than the core which is driving it; if this is not so 
the driving core will stop switching before the load core 
is completely switched. When a logic circuit core is driv- 
ing a load, the resultant drive on the circuit core is 
N.,I—N;I, since there are NJ ampere-turns driving 
the core “down” and Ny;J driving it “up.” (Na is the 
number of turns in an advance winding; N,;, forward 
winding; N;, backward winding; N;, load core input 
winding; and J, the advance current.) Sands [9] has 
shown that the variation of inverse switching time with 
magneto-motive force (drive) is linear. Therefore, NaI 
—N,I must be less than the drive to the load core, N,J, 
if the load core is to switch faster than the driving core. 
However, if the advance current were to split equally 
between two paths, the net drive on the driving cores 
would be (Na—4N,)I. If Na—43N;> Ni, then the driving 
cores may finish switching before the load core is com- 
pletely switched, but if no two paths conduct simul- 
taneously this difficulty does not occur. 

Thus the switching function to be realized is written 
in disjunctive form before translating it to circuit 
form. Each path corresponds to one product term of the 
sum of disjunctive products, so no two paths may ever 
conduct simultaneously. In Fig. 2 the paths have been 
designed to correspond to the disjunctive terms of the 
demonstration function, as described by the Karnaugh 
map [10] of Fig. 4. 

The above discussion also implies certain restrictions 
on the design of multiple-output switching circuits. 
Consider a circuit which is to drive two load cores. 
Load core 1 is to switch for Fi(x1, x2, x3) =1, and load 
core 2 for Fo(x1, x2, x3) =1. If there are some combina- 
tions of inputs for which both Fy and fF, are unity, 
there must be three mutually exclusive (disjunctive) 
outputs. Output 1 carries a current pulse when load 1 
alone is to switch (fi= Fi F2’); output 2 carries a pulse 
when load 2 alone is to switch (fe= Fi’ Fy); and output 
3 is pulsed when both load 1 and load 2 are to switch 
(fs= Fif2). Output 1 is connected to an input winding 
on load core 1; output 2 is connected to an input wind- 
ing on load core 2; output 3 is connected to an input 
winding on load 1 in series with one on load 2. Of 
course, the circuits for each output must be designed 
with all paths disjunctive. This restriction on multiple 
output circuits just described for two load cores may 
be generalized for any number of load cores. Fig. 5 
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Fig. 4—Karnaugh map of the switching function 
in disjunctive form, 
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Loads 


Output Network 


Fig. 5—Output network and loads of binary serial full adder. 


illustrates these principles. It is the output network of 
a binary serial full adder with inputs A and B added 
to Co, the carry digit from the previous time slot. The 
outputs are S, the sum, and C, the carry. The comple- 
ments of the inputs are furnished, and the circuit 
delivers the complements of the outputs. 

A further restriction upon the design of multiple- 
output circuits is that no magnetic core may furnish 
the forward winding for more than one output. Ob- 
-serve that this restriction is satisfied by the binary 
adder circuit of Fig. 5. Output CS’ is driven by forward 
windings on cores Cy and Cy’; CS is driven by B; C’S 
by A’ and A; and C’S’ by B’. (Refer now to Fig. 6.) If 
a driving circuit core (A) should have forward wind- 
ings (A; and A2) in the paths of two disjunctive loads, 
I, and Ly, the following situation might arise: when the 
inputs are such that L; should be switched, the forward 
core (A) in the path of Z; is driven by (Na—Ny)I 
ampere-turns; whereas a core with a backward winding 
(B’) in the path of ZL, is driven by N.J ampere-turns. 
Asa result, the “backward” core (B’) finishes switching 
much sooner than the “forward” core (A). If the 
“forward” core (A) has a forward winding (A2) in the 
path of ZL», then the forward voltage forces some cur- 
rent to flow in ZL, when none should flow, after the 
“backward” core (B’) finishes switching. The circuit 
of Fig. 5 operates properly because forward windings 
of different cores drive each output. 
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Fig. 6—lllustrating a restriction on multiple-output circuits. 


Diodes must be placed in output networks so that 
no circulating currents can exist. Such circulating cur- 
rents, whether occurring during the input phase or dur- 
ing the advance phase, could improperly switch cores 
with windings in their paths. One method of properly 
placing diodes is to assume one diode with the proper 
polarity is placed in each path between every pair of 
nodes. The absence of some of these diodes will permit 
circulating currents or permit current to flow in the 
wrong direction. The other diodes are superfluous. 


ECONOMICAL LOGICAL DESIGN 


An economical realization of a given switching func- 
tion incorporates a minimal number of magnetic cores, 
windings, diodes, and man-hours spent in logical design. 
The problem of minimizing magnetic cores is a fasci- 
nating topic, and investigating core minimization does 
shed some light on other problems. 

In Fig. 2 the demonstration function is realized using 
only five cores. If the output windings were all forward 
or all backward, eight cores would be needed—primed 
and unprimed inputs for all four variables. The follow- 
ing discussion describes a useful technique for minimiza- 
tion. It is possible to form three of the four combinations 
of two variables and their complements using only two 
magnetic cores, in conjunction with other cores repre- 
senting the other variables. In Fig. 7 cores x; and x; are 
used to deliver product terms involving x;,'x;’, x;’x;, and 
x,x;’. The restriction that no more than one of these 
two cores may be wound in the forward direction is the 
reason that x,;x; cannot be produced. If all four combi- 
nations are required by any function, then at least 
three cores must be used. In the light of these facts, a 
switching function may be checked to determine which 
pairs of variables appear with three different combina- 
tions of primed and unprimed terms, and which ones 
appear with all four combinations. This indicates for 
which pairs of variables two cores may suffice and for 
which pairs three are needed. It is then necessary to 
test the various possibilities. One of several different 
expressions (disjunctive) of the given function may 
yield a network using fewer cores. 

The final step in minimization is to reduce the num- 
ber of output windings. This reduction is similar to the 
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Fig. 7—Producing product terms with two cores. 
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Fig. 8—Typical circuit for producing a function 
of three variables. 


reduction of the number of contact pairs in relay con- 
tact networks, with the exception that output windings 
conduct current in one direction only. 


EXPERIMENTAL STUDIES 


Fig. 8 is the schematic of an experimental circuit 
which sets a load core, L, when the inputs satisfy the 
equation F=x1%2'x3' +-x1x%2%3-+%1'x3=1. (Note that this 
function can be produced with three cores—Fig. 9.) 
Arnold Engineering Company permalloy tape cores 
(20 wrap, + mil) and Transitron Electronic Corporation 
T7G diodes were used. The number of turns per winding 
were calculated according to well-established princi- 
ples [11]. Negative current pulses of 100 milliamperes 
amplitude and 2 microseconds duration were applied 
alternately at the terminals labeled ¢; and ¢. During 
¢; one of each pair of cores representing primed and un- 
primed variables was set, depending upon the positions 
of the input switches; also, the load core was reset. Dur- 
ing ¢2 the output was delivered to the load. The switch- 
ing of the load core was observed across 20 turns con- 
nected to an oscilloscope. Fig. 10 is a series of traces tak- 
en from oscilloscope photographs. They represent the 
voltage across the sensing winding on the load core 
(shown in the lower right portion of Fig. 8) during ¢1. 
The traces correspond directly to the Karnaugh map, 
shown to the right. Notice the small “blips” which ap- 
pear, whether or not F=1, at the beginning and end of 
the input pulses of ¢:. These are shuttle voltages (de- 
scribed previously [5]) which cannot be eliminated. 

A circuit was constructed to produce the demonstra- 
tion function of four variables previously discussed. 
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Fig. 10—Output waveforms for core circuit producing 
function of three variables. 


This circuit also performed most satisfactorily. When 
the circuits were built with different cores, or when the 
cores were interchanged, no difference in behavior could 
be observed. This interchangeability is the feature of 
the new system which makes it so ideal for application 
to standard packages. 

One problem encountered was the slight switching of 
load cores which occurred when no output signal was 
supposed to be delivered to the load. The reason for 
this malfunction is that, for some combinations of inputs, 
when F=0 none of the cores corresponding to windings 
in a path may be switching, so no net voltage is gener- 
ated along the path. This occurs for the top path of the 
output network of Fig. 8 when (x1, x2, x3) =(0, 0, 0). 
Therefore a part of the advance current may sneak down 
this forward path and partially switch the load instead 
of flowing through the shunt diode. When an extra 
diode was placed in series with the input winding of the 
load core, the partial switching was greatly limited. 
This diode effectively increases the impedance of the 
load path relative to that of the shunt path. 

It was originally assumed that N, must be greater 
than Ny,, to insure a net backward voltage when a core 
with a forward winding and one with a backward wind- 
ing in the same path are switching simultaneously. Ex- 
perimental studies revealed that making N; equal to Ny; 
did not appreciably increase the noise seen at the out- 
put when F=0. Hence, all experimental circuits were 
constructed with N,=N;=16 turns. This significant 
fact implies that only three types of windings are needed 
on any core: input windings (6 turns in the experimental 
circuits), one advance winding (21 turns), and forward- 
backward windings (16 turns). The only difference 
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between forward and backward windings is the way 
they are connected into circuits; 7.e., a forward winding 
becomes a backward winding if its terminals are inter- 
changed. Since only three types of windings are re- 
quired, the expense of winding the cores is considerably 
reduced. 

Tests were made to see how sensitive the experimental 
circuits were to changes in the characteristics of the 
current pulses furnished by the core driver. Circuit 
behavior was not significantly affected by moderate 
variation in pulse width, pulse amplitude, and pulse 
rise time. An increase in pulse width above the necessary 
minimum did not change the response of the circuits. 
The test circuits also operated properly for a very large 
range of pulse amplitude. Increasing J above the speci- 
fied value merely speeded up the circuit operation; 
cores switched faster, but the circuits still operated ac- 
cording to logical specifications. 


STANDARD PACKAGES 


The actual packaging of core circuits must be done 
with the over-all system in mind. After the designer 
has derived the proper number of turns for the three 
types of windings, he must decide upon convenient 
numbers of cores to include in the packages, and how 
many input windings and forward-backward windings 
to wind on each core. His object is to be satisfied with 
as few different types of packages as are consistent 
with economical use of the components. 

The components might be cast in an epoxy type of 
resin. (Librascope manufactures its “Decision Elements” 
in that form.) An actual 6-core package might consist 
of 6 cores and 5 diodes. There would be an advance 
winding on each core, and the cores would be connected 
in series, with just the two extremities appearing at 
surface terminals. The input windings and forward- 
backward windings would be terminated at the surface. 
Diodes also would be terminated at the surface of the 
package, and there might be several sets of shunted 

_ terminals to serve as nodes. A particular switching func- 
tion would be produced by properly interconnecting 
the terminals on a standard package. 


COMPARISONS AND CONCLUSIONS 


Among existing systems of core logic, the ones most 
suitable for standardization are those of Minnick [4], 
Karnaugh [5], and Andrews [6]. These systems do have 
their drawbacks. Minnick was one of the first (1953) in- 
vestigators of magnetic core pulse-switching. The circuits 
he described generally require several clock pulses per 
cycle of operation; they cycle all cores once every cycle 
of operation, deliver unwanted outputs before the de- 
sired output is delivered, and use a comparatively large 
number of components (cores, diodes, windings). Kar- 
naugh’s “A-type” circuits remedied these shortcomings. 
However, the “A-type” circuits demand that all input 
pulses arrive simultaneously. In devices where input 
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pulses must travel from different physical locations, it 
is likely that the inputs will not arrive simultaneously. 
Andrews, with his “shunt-type” circuit, solved this 
problem by moving the logic function from the input 
windings to the output windings. In a “shunt-type” 
circuit inputs may arrive in any order whatsoever. The 
difficulties with “shunt-type” logic are an unwanted 
pulse of output current which appears at the beginning 
of the advance pulse, and a greater number of cores 
than found in “A-type” circuits. 

These final problems are solved in the new type of 
logic circuits. No noise pulse appears. Furthermore, the 
number of cores, diodes, and windings which are used 
in the new system compare favorably with the numbers 
required in equivalent “A-type” circuits. The primary 
disadvantage of the new system, however, is the rela- 
tively large number of turns which must be wound on 
the cores. 

The salient advantage of the new system is the ease 
with which circuits can be constructed. Once all cores 
are wound with the appropriate number of turns in each 
winding, circuits are built by interconnecting windings 
according to sketches similar to Fig. 2. The fact that any 
core may be chosen to represent any variable—.e., that 
circuit behavior is not dependent upon which cores are 
used where—makes the new type of logic an especially 
good choice for use in a system of standard packages. 

The experimental studies of the new system demon- 
strated the practicability of building combinational 
switching circuits using advance current drive, with 
logic performed by forward and backward windings in 
an output network. It was observed that such circuits 
may be designed to be independent of moderate changes 
in driving pulse amplitude, duration, and rise time. 
Moreover, the behavior is not affected by interchanging 
components. 
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The Switching Characteristics of 4-79 Permalloy 
Cores with Different Anneals* 
T. D. ROSSING}, W. M. OVERN, ann V. J. KORKOWSKIt 


Summary—The magnetic properties of 4-79 Permalloy cores, 
which have been annealed at different temperatures, are discussed. 
Cores annealed at relatively low temperatures are characterized by 
high coercivity, slower switching, insensitivity to strain, mag- 
netization difficult to rotate, and insensitivity to an applied trans- 
verse field. Some cores exhibit a preference to one remanent state 
over the opposite state. 


URING the attempted development of a memory 
D with nondestructive readout, investigations 
have been made of the switching characteristics 
of 4-79 Moly-Permalloy, ferrites, and vacuum-deposited 
films. For reasons which will become apparent in this 
paper, Permalloy tape cores are not too promising in 
this application. However, some of the data gathered in 
these investigations sheds new light on the remagnetiza- 
tion process in ferromagnetics. 

Several experiments! have shown that the wall- 
motion switching model describes the remagnetization 
of Permalloy tape cores. For switching times of 10 
microseconds or less, the switching time is related: to 
the applied field H by the equation: 


(H — Hy)i = S. 


Hy is approximately the same as the coercivity H,, and 
Sis a switching coefficient which depends upon the num- 
ber of domains of reverse magnetization nucleated and 
the speed at which the walls of these domains move 
through the specimen. 

The structure of the cores used in this investigation is 
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shown in Fig. 1. Two $-mil 4-79 Permalloy ribbons, 


% inch wide, are spot welded to stainless-steel straps 
coated with insulation. They are then wound around the 
bobbin and spot-welded to a stainless-steel band which 
constitutes the outer wrap of the core. By means of the 
attached current leads, current may be passed through 
the wraps of the core. The magnetic field produced by 
such a current will be everywhere transverse to the 
magnetization of the ribbon. These cores were annealed 
for one hour in a hydrogen atmosphere at various tem- 
peratures ranging from 425° to 900°C. The resulting 
coercivities of these cores are shown in Fig. 2. It is 
interesting to note that the lower the coercivity of the 
core, the more sensitive it becomes to strain. Cores 
with coercivities of 0.6 oersted or more have been 
completely unwrapped and rewrapped after annealing 
without any appreciable change in magnetic properties. 

The switching coefficient S also depends upon the 


coercivity of the core. In general, the higher the coerciv- 


ity of the core, the larger S becomes. Fig. 3 is a plot 
of S as a function of H, for cores from various anneals. 
On the other hand, as the coercivity is lowered, the 
transverse field effect on the switching behavior is 
diminished. For faster switching, therefore, using the 
coincidence of longitudinal and transverse field pulses, 
the optimum coercivity will be determined by a com- 
promise between these two effects. In most high-speed 
computer memories, writing is accomplished by the 
coincidence of two “half-select” pulses which are current 
pulses of half amplitude. The same compromise must 
be made, therefore, when writing is accomplished by this 
coincident-current technique since the switching speed 
increases as H, is raised but decreases with increasing S. 

A most disturbing feature of tape-wound cores of 
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Fig. 2—Coercivity as a function of anneal temperature. All anneals 
were made for one house in a hydrogen atmosphere. 


4-79 Permalloy is their tendency to “creep” into the 
demagnetized state or into the opposite remanent state 
when a large number of “half-select” pulses are applied. 
In many cases, creep can be brought about by applying 
pulses smaller than the 60-cycle coercivity H,. To 
minimize magnetic creep in the unselected cores, bi- 
directional writing pulses have been used. The relative 
amplitudes of the positive and negative pulses, however, 
must be carefully controlled. 

A block diagram of a dynamic core tester is shown in 
Fig. 4. “Read” pulses consist of 0.1-usec current pulses 
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Fig. 3—Switching coefficient S as a function of coercivity. 
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Fig. 4—Block diagram of a dynamic core tester. 


applied through the magnetic tape itself to produce a 
transverse magnetic field. The magnetic state of the 
core thus may be nondestructively sensed at any point 
in the test cycle. “Shaking” pulses are bursts of these 
same 0.1-usec pulses applied at the 2-mc rate. These 
shaking pulses produce the transverse field which 
switches a core when it coincides with a switch (longi- 
tudinal field) pulse. In a typical test program, a restore 
pulse is followed by a read, a write pulse, another read, 
and then a large number of disturb pulses alternating in 
direction. Another read pulse then shows how much 
change in magnetization (creep) has been brought about 
by the disturb pulses. This test program is shown in 
Figio: 

Fig. 6 shows the results obtained from one of the most 
stable cores tested with write and disturb pulses 7 usec 
long. The stable region is bounded by curves showing 
the maximum amount by which the second half of the 
bidirectional write pulse may be larger or smaller than 
the first half without causing creep. Also indicated are 
the switching fields necessary to switch the magnetic 
state of the core in 7 usec with (H,) and without (7,) 
a transverse shaking field. Also shown, for purposes of 
comparison, are the coercivity H, and Hmin, the field 
below which the hysteresis loop completely collapses. 
It is interesting to note that in this core, writing could 
be accomplished with a 0.49-oersted switch pulse and 
a burst of shaking pulses of 800 ma. Coincident-current 
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Fig. 5—Typical pulse program used to test for magnetic 
“creep” with the dynamic core tester. 
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Fig. 6—A stability plot obtained from one of the most stable cores. 
The two curves indicate the maximum difference in magnetic 
field pulses which can be tolerated in the “1” and “0” states of 
remanence. 


writing would require half-select (} H,,) pulses of about 
0.49 oersted. To prevent creep, the pulse amplitudes 
would have to be controlled to within plus or minus 4 
per cent; 

An average core as regards creep behavior is shown 
in Fig. 7. In this core pulses would have to be controlled 
to plus or minus 1 per cent or more. An interesting effect 
which is shown in Fig. 7 is that there is a preferred state 
of remanence into which the core will creep more readily 
than into the other state. Some such preference was 
demonstrated by nearly half the cores tested, and was 
independent of the orientation of the cores in the tester. 
Magnetic fields as large as 200 oersteds did not reverse 
the preferential direction, but they increased or dimin- 
ished the amount of preference, depending upon the 
direction in which they were applied. Meiklejohn and 
Bean? have observed a similar preferential direction 


2 W. H. Meiklejohn and C. P. Bean, “New magnetic anisotropy,” 
Phys. Rev., vol. 102, p. 1413; June, 1956. 
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g. 7—A stability plot for an “average” core. 
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Fig. 8—Outputs obtained when nondestructive readout interrogat- 
ing pulses are applied to cores of different coercivities. Such out- 
puts are measures of the ease of rotation of the remanent mag- 
netization. 


phenomenon which, they explain, is an interaction be- 
tween an antiferromagnetic and a ferromagnetic ma- 
terial. 

Another property of the cores which depends on 
coercivity is the ease with which the magnetization 
can be rotated by a transverse field. A convenient meas- 
ure of this property is the amplitude of the nondestruc- 
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tive readout signal obtained when a 0.1-ysec transverse 
magnetic pulse is applied. Fig. 8 is a plot of readout 
signals as functions of applied 0.1-usec pulses for cores 
of varying coercivity. Apparently, in cores of low co- 
ercivity the magnetization may be rotated reversibly 
with more ease than in cores of high coercivity. This 
same effect has been noted in ferrite memory cores. 


CONCLUSIONS 


Cores of 4-79 Molybdenum Permalloy can be made to 
have a wide variety of magnetic properties depending 
upon at what temperature they are annealed. An an- 
neal at about 575°C produces the best cores for memory 
application. 
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Cores annealed at relatively low temperatures are 
characterized by: 1) high coercivity, 2) large value of 
the switching coefficient S, 3) insensitivity to strains, 4) 
magnetization relatively difficult to rotate, and 5) 
switching time quite sensitive to an applied transverse 
field. The properties of cores annealed at high tempera- 
tures are the opposite. A slight preference to one mag- 
netic state has been noted in both types of core. 
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Formal Analysis and Synthesis of Bilateral 
Switching Networks” 


RAYMOND E. MILLERt 


Summary—Formal procedures for the analysis and synthesis of 
two-terminal combinational bilateral switching networks are pre- 
sented. A bilateral switching network is one which contains only ele- 
ments having the same switching transmission characteristic in both 
directions. 

Following the definitions for the terminology and notation, where 
some new terms are introduced, the definitions for a series-parallel 
network, a bridge element, and a bridge network are given. A 
condition, called the bridge condition, to test a given transmission 
function for possible bridge network realizations is presented. 

A stepwise decomposition procedure is developed which may be 
used for the analysis and synthesis of the series and parallel parts of 
the network. The steps are described both with linear graphs and 
connection matrices. The bridge condition partially formalizes bridge 
network synthesis. Redundant variables also are considered as an 
aid to network synthesis. Under certain conditions, the synthesis 
yields network realizations with the fewest possible number of ele- 
ments. 


I. INTRODUCTION 


HE application of Boolean algebra to the analysis 
‘De design of switching networks is quite well 
known. A Boolean function can be realized by 
an unlimited number of different switching networks. 
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The network chosen usually depends on some type of 
economy consideration. One type of economical switch- 
ing network design is one which uses the least possible 
number of elements; such a network is called a simplest 
switching network. The general problem of finding such 
a simplest switching network design to represent a 
given Boolean function has not yet been solved. The 
synthesis of simplest switching networks is compli- 
cated by the facts that many essentially different net- 
work configurations are possible to represent the func- 
tion, the simplest network is not necessarily unique, 
and that it is extremely difficult to show that a network 
is indeed a simplest one. 

This paper considers the problem of designing sim- 
plest switching networks fora special class of switching 
networks. An analysis and a synthesis procedure are 
presented for two-terminal combinational bilateral 
switching networks. A bilateral switching network is 
one which contains only elements having the same 
switching transmission characteristics in both directions. 
Relay contact networks are the most common example 
of bilateral switching networks, whereas electronic logi- 
cal networks do not normally fall into this class of 
switching networks. The network synthesis procedure 
given here is developed to yield simplest switching 
networks, and in certain cases it is proved to do so. 
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The process is capable of synthesizing bridge networks 
and nonplanar networks as well as series-parallel net- 
works. 

The basic postulates and theorems of Boolean algebra 
are assumed to be familiar to the reader.!:? The formu- 
las of Boolean algebra are to have the letter variables 
a, b,c, +++. The operation of complementation or nega- 
tion is denoted by a “—,” for example, 6; that of con- 
junction or intersection by juxtaposition, for example, 
xy; and that of alternation or union with a “\/,” for 
example, c\/d. 

Letter variables and their negations are called 
literals. A conjunction of literals such that no literal 
appears twice and no literal and its negation appear is 
called a clause. If a literal and its negation appear in a 
conjunction of literals, then the conjunction is called 
a zero value clause. A literal itself may be called a clause. 
An alternation of clauses is called an \V polynomial. A 
clause itself may be considered to be an \V polynomial. 
A simplest equivalent to an \/ polynomial is called a 
minimum \/ polynomial. Quine calls this a simplest 
normal equivalent. 

Two functions, f; and fe, are said to be matchable or 
to have a matching if and only if for each clause in fi 
the same clause appears in fz, and vice versa. 

Given an \/ polynomial f such that 
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then the right member is called an assembling of f if and 
only if upon expanding gigs - - - g», using only the dis- 
tributive and commutative laws, the resulting right 
member is matchable with the given f. A function is 
called fully assembled if and only if all of its clauses 
are used in the assembling process. It follows in this 
case that f’=0. 

A conjunction of some of the literals appearing in a 
clause is called a subclause of that clause. A clause may 
be considered to be a subclause of itself. A conjunction 
(without simplification) of subclauses 1s considered to be 
a clause or a subclause only if no literals appear more 
than once in the conjunction. In the following sections, 
all conjunctions of subclauses are considered to be 
clauses or subclauses unless specifically stated otherwise. 
Clauses and subclauses are denoted by capital letters 
and conjunctions of capital letters, for example, 4, 
ACE, and: SRE. 

The Boolean transmission functions will usually be 
denoted by the symbol f or 7, and are not to be con- 
fused with the symbols for a literal or a subclause. 

The minimum \/ polynomial is most generally used 
as the starting function for network synthesis. Various 


f= fig 


1G, Birkhoff and S. MacLane, “Algebra of Classes” in “A Survey 
of Modern Algebra,” The Macmillan Co., New York, N. Y., ch. 
11, pp. 311-325; 1948. 

2 R, Serrell, “Elements of Boolean.algebra for the study of in- 
formation handling systems,” Proc. IRE, vol. 41, pp. 1366-1380; 
October, 1953. 

3 W. V. Quine, “A way to simplify truth functions,” Amer. Math. 
Month., vol. 62, pp. 627-631; November, 1955. 
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procedures for obtaining a minimum \/ polynomial for 
a given Boolean function are well known.*~° 

The switching networks are assumed to be depicted - 
as weighted linear graphs. The weights are Boolean 
functions attached to the elements of the graph. An 
element consists of a line segment and its endpoints. A 
vertex is an endpoint of an element, and a linear graph 
is a collection of elements such that no two elements 
have a point in common which is not a vertex. A sub- 
graph is a graph containing a subset of the elements of 
the graph. A vertex and an element are incident with 
each other if the vertex is an endpoint of the element. 
The vertices of the graph may be classified as either 
internal vertices or terminal vertices. The transmission 
functions represent the transmission between terminal 
vertices. Switching networks represented by such 
weighted graphs are known as combinational switching 
networks. 

The switching network may also be represented by 
a connection matrix.°-? For a graph with p vertices 
the connection matrix is a pXp matrix. The matrix 
entries are Boolean functions. An entry c;; represents 
the connection from vertex 7 to vertex 7. The c;; entries 
are all defined to be the Boolean 1. If c;;=0, this means 
that no direct connection exists between vertices 7 and 
j. If ci; =1, this means that a short circuit exists between 
vertices 2 and j. Now if c,;;=0, 1, a literal, or an alterna- 
tion of literals for each c,;; of a connection matrix, then 
the matrix is called a primitive connection matrix. 

A network with two terminal vertices is called a two- 
terminal network. For convenience, the input vertex is 
labeled vertex 1, and the output vertex is labeled vertex 
2. The weights of elements incident to the input vertex 
thus appear as entries in the first row and column of the ' 
corresponding primitive connection matrix. The weights 
of elements incident to the output vertex appear as en- 
tries in the second row and column. 


4 Harvard Univ. Computation Lab. Staff, “Synthesis of Elec- 
tronic Computing Control Circuits” in “Annals of the Computation 
Oe Harvard University Press, Cambridge, Mass., no. 27; 

5M. Karnaugh, “The map method for synthesis of combinational 
logic circuits,” Trans. ATEE (Commun. and Electronics), vol. 72, pp. 

593-599; November, 1953. 

8 EK. J. McCluskey, Jr., “Minimization of Boolean functions,” 
Bell Sys. Tech. J., vol. 35, pp. 1417-1444; November, 1956. 

og, DOE Muller, “Minimization of \V Polynomial Subject to Sub- 
sidiary Conditions,” Digital Computer Lab., University of Illinois, 
Urbana, Il., Internal Rep. No. 51; October, 1953. 

8 J. P. Roth, “Algebraic Topological Methods for the Synthesis 
of Switching Circuits in Variables,” Inst. for Advanced Study, 
Princeton, N. J., Electronic Computer Project, Tech. Rep. No. 
56-02; April, 1956. 

* R. H. Urbano and R. K. Mueller, “A topological method for the 
determination of the minimal forms of a Boolean function,” IRE 
TRANS. ON ELECTRONIC CompuTERs, vol. EC-5, pp. 126-132; Sep- 
tember, 1956. 

0B. I. Aranovich, “The use of matrix methods in problems of the 
structural analysis of relay contact networks,” Avtomatika i Tele- 
mekhanika, vol. 10, no. 6, pp. 437-451; 1949. 

_ 1\F. E. Hohn and R. L. Schissler, “Boolean matrices and the de- 
sign of combinational relay switching circuits,” Bell Sys. Tech. J., 
vol. 34, pp. 177-202; January, 1955. 

_ 2 W. Semon, “Matrix Theory of Switching Networks,” Computa- 
tion Lab., Thesis, Harvard University, Cambridge, Mass.; 1954. 
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II. SERIES-PARALLEL AND BRIDGE NETWORK 
CLASSIFICATIONS 


The switching networks to be analyzed and synthe- 
sized are conveniently classified into two types, the 
series-parallel, and the bridge networks. The Boolean 
functions may be directly represented by series-parallel 
networks, and vice versa. Boolean functions also may be 
used to represent the transmission in a bridge network; 
however, the representation of a Boolean function by a 
bridge network is not direct. The definitions of the net- 
works show that a network is either a series-parallel or a 
‘bridge network, but never both, for a specified input 
vertex and output vertex. 
| Two equivalent definitions for series-parallel net- 
‘works follow.!3 
_ Definition 1: A network N is series-parallel with re- 
spect to the two terminals a and 0 if through each ele- 
‘ment of N there is at least one path from a to b not 
touching any junction twice, and no two of these paths 
‘pass through any element in opposite directions. 

Definition 2: A network is series-parallel if it is either 
a series or a parallel connection of two series-parallel net- 
works. A single element is a series-parallel network. 

To define a bridge network, a definition is given for a 
‘bridge element in a network. 

Definition 3: An element of a graph representing a 
‘switching network is a bridge element with respect to the 
assumed input and output vertices if and only if it 
belongs to a set of elements having the properties: 


Removal of this set from the graph reduces it to a 

| series-parallel network; 

_* When the set of elements is removed from the graph 
no internal vertex has fewer than two incident ele- 
ments; 

_ No proper subset of elements of this set when removed 

| from the graph will give a series-parallel network; 

_ No element of the set is incident to either the input 

| or the output vertex. 


Now, a bridge may be simply defined. 
_ Definition 4: A network is a bridge network if and only 
if it contains at least one bridge element in its graphical 
representation. 
Some examples of bridge networks are now given. 
- Example 1: Removal of element e of Fig. 1 yields a 
series-parallel network. Element e is a bridge element 
with respect to the terminal vertices 1 and 2. No other 
elements are bridge elements. 
_ Example 2: The bridge element sets for the network 
of Fig. 2 are: fe}, daft This example shows that more 
than one bridge element set may exist for a given graph, 
and also that the number of elements in each bridge 
element set need not be a constant for a given network. 
Some known facts for switching networks are par- 


8 J. Riordan and C. E. Shannon, “The number of two-terminal 
series-parallel networks,” J. Math. and Phys., vol. 21, pp. 83-93; 
August, 1942. 
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Fig. 2—A bridge network having two bridge element sets. 


ticularly useful in the synthesis procedures to be de- 
veloped. These are stated as three theorems. 

Theorem 1: A network with a minimum number of ele- 
ments which represents a Boolean transmission contains 
at least as many elements as the number of different 
literals appearing in the minimum \/ polynomial repre- 
sentation of the function. 

This follows from the definition of a minimum \V/ 
polynomial.¥ 

Theorem 2: If a switching network with bilateral ele- 
ments is represented by a connection matrix A, cross 
out the input vertex column and the output vertex 
row. The determinant formed from this resulting minor 
yields the transmission function. 

This determinant is called an output determinant for 
the network. 

The determinant of a Boolean matrix, An, is defined as 
the alternation of the 2! conjunctions of the entries of 
A,» in which each row subscript and each column sub- 
script is repeated once and only once in each conjunc- 
tion. 

This theorem is due to Aranovich.!? and Semon.!? 
An example demonstrates this theorem. 

Example 3: Consider the network graph and con- 
nection matrix C shown in Fig. 3. 

Now, letting 7;; denote the transmission from vertex 
4 to vertex 7, 


OnnO a 5 ae 
(Pip = CEG | Who Ak = oaet| * | = cava 
C 
i. "ee 4 


“RE. Miller, “Formal analysis and synthesis of bilateral switch- 
ing networks,” Doctoral dissertation, University of Illinois, Urbana, 
Ill.; 1957. The proof of this theorem and the other results ofthis 
paper appear in this dissertation. 
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Fig. 3—A network and its associated connection matrix. 


Theorem 3: Let A be the connection matrix for a 
network. Replace each nonzero nonunit entry by the 
corresponding symbol @yq. Form the determinant for the 
transmisssion function by theorem 2. Then an element 
represented by a;,=a;; is a bridge element if and only 
if each of a,; and a;; appears in a separate clause of the 
expanded \/ polynomial representation (7.e., not in the 
form a;;a;;), where the redundnat clauses are removed. 

Semon states this condition as a definition for a bridge 
element. It is stated here as a theorem and can be 
shown to follow from Definition 3. 


II]. THe BriDGE CONDITION 


It would be desirable to be able to test a given trans- 
mission function in some manner to determine whether 
the function might be representable as a bridge net- 
work with a saving in elements over a series-parallel 
network. A condition of this sort is given now and it is 
used in a later section for the synthesis of bridge net- 
works. 

Let A represent a subclause of an \V polynomial, 
equal to the function f. If A is representable as a bridge 
element, then the function f may be expressed as 


f = A(BC\V DE) \ zg, 
where 


1) A, B, C, D, and E represent subclauses or literals 
appearing in f or represent redundant subclauses 
or literals. 

2) ABC <the function f. 

3) ADE<the function f. 


4) gis the remaining part of f with the dBC and ADE 


related clauses removed. 
BD is a clause or a subclause of g, and 
5) « CE isa clause or a subclause of g, where zero must 
always be considered as a clause of g. 


Note that 5) may be restated as 


Ai = some clause of g, and 


CE = some clause of g. 


The symbol < 
includes. 

This condition can be shown to be necessary and 
sufficient for A to be representable as a bridge element. 
Some examples will demonstrate the condition. 


means is included in, and > means 
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Fig. 5—-The bridge network for Example 5. 


Example 4: 
f= abe V aden bd \/-ce 
= a(bc V de) VV bd V/ ce. 


Note this function satisfies the bridge condition. The 
bridge network that realizes this function is shown in 
Fig. 4. 

Example 5: Consider the network graph in Fig. 5. 
For this network 


f = xbd \/ xabc \/ xce \/ xead \/ yd \/ ybec \/ yac. 
= a(xbc \/ xed) \/ xbd \/ xec \/ yd \/ ybec \/ yac. 
Since 
(«b)(d) = xbd is a clause in g, and 


(c)(xe) = xec is a clause in g, 
the element a is seen to satisfy the bridge condition. 


IV. THE DECOMPOSITION PROCEDURE 


A systematic procedure for decomposing a _ two- 
terminal switching network is described which provides 
a method for classifying a given network as either a 
series-parallel or a bridge network. As described in the 
next section, the decomposition procedure also is useful 
for synthesis. 

The decomposition procedure is characterized in the 
form of a graphical analysis, and also as a partitioning 
and decomposition of the primitive connection matrix. 
The matrix representation provides a convenient repre- 
sentation for computer use. Let the connection matrix ~ 
be represented by C and the network graph by G. 
When elements are removed during decomposition, — 
causing a vertex to become isolated, the vertex is also 
considered removed from the network. 

The decomposition procedure is stated in five steps. 

Step (A)—Removal of Series Elements: If in G a ter- 
minal vertex has only one incident element, remove this 
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lement to form a new graph G’ with the terminal 
ertex of G’ being that vertex which was incident to 
he removed element. Reconsider G’ as a new G. Repeat 
intil no further elements may be removed. 

Step (B)—Removal of Parallel Elements: Remove all 
slements through which one and only one path exists 
( kegs terminal vertices. Let the new graph become 
1. new G. The network configuration for any one of the 
sets of elements in such a path is simply a series network 
rom the input vertex to the output vertex. 

Step (C)—Separation into Series Subnetworks: In G, 
f all paths from the input vertex to the output vertex 
Dass through a given internal vertex 7, then break G 
mto two subnetworks G’ and G’’, where G’ has input 
vertex 1 and output vertex r and G’’ has input vertex 7 
ind output vertex 2. Consider each G’ and G’’ as a new 
3, renaming the output vertex of G’ as 2, and the input 
vertex of G’’ as 1. Repeat this step until no further 
splitting occurs. 

Step (D)—Separation into Parallel Subnetworks: If all 
the internal vertices of G can be separated into two sets 
such that no element is incident to vertices in both sets, 
then break G into two subnetworks G’ and G’’, where 
each subnetwork consists of all the elements and ver- 
‘ices connected to the internal vertices of one of the sets 
of vertices. Each subnetwork has an input vertex 1 
und an output vertex 2 formed at the respective end- 
points of the elements incident to the terminal vertices 
of G. Let each G’ and G’’ represent a new G. Repeat 
‘his step until no further splitting occurs. 

Step (E)—Repetition: Repeat Steps (A) through (D) 
»n each G until no further reduction can be accomplished 
m any step. 

_ This procedure reduces the original network to a set 
of subnetworks or to the null set of elements. 

- The following properties connected with the decom- 
osition procedure can be shown to hold. 

Property 1: The decomposition procedure reduces 
‘he two-terminal network to a null set of elements if and 
only if the network is a series-parallel network. 

Thus, if the decomposition procedure leaves at least 
yne subnetwork which cannot be reduced further, then 
he original network is a bridge network. 

Property 2: A series network is reduced by Step (A) 
alone to the null set. 

_ Property 3: A parallel network is reduced by Step (B) 
lone to the null set. 

Property 4: Each subnetwork remaining after the 
-ompletion of the decomposition procedure is a bridge 
1etwork. 

Property 5: In each subnetwork remaining after the 
ompletion of the decomposition procedure there exists 
. “circuit” (i.e. a closed path) containing the input ver- 
ex and not the output vertex of that subnetwork, and 
. “circuit” containing the output vertex and not the 
nput vertex such that these circuits have at least one 
‘ommon element. 

Property 6: The decomposition procedure, when ap- 
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plied to a network, reduces the network to a unique set 
of subnetworks or to the null set of elements regardless 
of the order of the decomposition. 

The characterization of the network decomposition 
with the aid of connection matrices follows. Let a primi- 
tive connection matrix be 


1 Cho a 3 e Cin 

Chae 1h C23 C2n 
C= 

Cn1 Cn2 1 


where 


1) m is the number of vertices of the network; each 
vertex is assigned a row and column in C, 

2) C=C; for 447, and 

3) Each c¢,; for 147 is 0, a literal, or an alternation of 
literals. 


Step (A): If one and only one 4;40; j=2, 3,---,7n, 
say Cip#0, then remove row and column 1 from C. 
Permute row and column p, of C, to the first row and 
column of the remaining matrix of order 7—1 to obtain 
a new connection matrix C’. Let C’ be a new C and 
repeat the above process until no further reduction is 
possible. Repeat similarly on row and column 2. 

Step (B): Two separate inspection steps are required 
to accomplish this step in the matrix representation. 
They are called (1B) and (2B). 

Case (1B): If ¢1:=¢2140, replace them with zeros. This 
removes the single element paths from vertex 1 to 
wettex 2: 

Case (2B): If ¢1a40, where a=3, 4, -- -, ”, and row 
a has only one other nonzero entry, excluding cai and 
Caa, Call it ca» 40, and then inspect row 0b. If row 6 has 
only one other nonzero entry, excluding c_ and cp, call 
it Goa%0, and then inspect row d. If row d has only one 
other nonzero entry, excluding Ca and Caa, call it cae 40, 
where ea. Inspect row e. Continue this process, never 
repeating a vertex, until some vertex 7 is encountered 
such that c,.#0. This ends the sequence. The sequence 
iS Cia) Cad, Cod) Cde, * * * » Crg: Where in general for the c;,’s, 
no 7 equals any other z and no j equals any other 7. If 
such a sequence is obtainable, remove rows and columns 
a, b, d,e,- ++, and, from the connection matrix to ob- 
tain a new connection matrix C’. Let C’ be a new C and 
repeat (2B) until no further reduction is possible. 

A partitioning of C may be made as follows if and 
only if Case (2B) applies. We have 


CliinGe C13 
| 
— | a ———— 
Cs C21 C2 | C2 
| 
— —- | 5 
C31 C2 | C33 


236 


where rows and columns 3 through m have been per- 
muted in C in such a way that: 


» ome 
0 


2eGe = Urand Cs — (0), 
3) All rows of 
C22 | 


| C1 | C28 | 
| | 
have two and only two ¢;;4#0; 747, and, 
4) All columns of 


0 
| , since ¢1.=0 from step (1B), 


have two and only two ¢;;40; 747. 
Applying Case (2B) gives the new connection matrix 
C’, where 


Let C’ be a new C and repeat the partitioning process 
as often as possible. 

Step (C): Step (C) is applicable to the primitive con- 
nection matrix if and only if there exists a partitioning 
—after suitable permutations of the rows and columns, 
if necessary—such that 


Cu C2 : C33 
| 
| — 
& — C21 C22 C23 5 
| 
= | 
C31 | C2 C33 
| | 
where 
| | 
1) | CP priecCZo git §C28 | alone 
| | 
=the row for vertex 7, 
Cl 
[pate | Spal | | 
2) C22 —= C21 | C22 | C23 


=the column for vertex 7. 


3) C¥=0 and C#=0, 
| 
4) | Ci | Ci | C18 | 


contains the input vertex row, and 
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ee : C32 : 
| | 


C33 | 


contains the output vertex row. | 


The two subnetworks G’ and G’’ are then defined 
by the new primitive connection matrices C’ and C”’, 


respectively. 
| C22 


| 
| 
C’ = and C” = ——— 
| 
| 
| 


C82 


The 


column and the 


| 
[cn | 
| 


row then are permuted to the second column and row 
of C’, which is the normal position designated for the 
output vertex in this analysis. This gives the conven- 
tional connection matrix C for the reduced network 
represented by C’. The output vertex C’’ (1.e., the 
vertex 2 in C) is similarly transposed to the second row 
and column of C’’ to yield a new C for further analysis. 

Step (D): A partitioning of C, as shown below, exists 
if and only if Step (D) applies. 


a Ey a b 
il ! Cll Ce C13 

Ah ce egies 
eae C?1 C2 | (28 
b | C3l C2 | C33 


where 


1) Cus E at) 
C21 1 . 


The rows and columns denoted by a are from one 
set of vertices of Step (D), and the rows and 
columns denoted by 6 are from the other set of 
vertices of Step (D), 

6), C#=0and +e? =0: 


The primitive connection matrices for G’ and G’’, 
designated by C’ and C’’ respectively, are 


bo 
7 


1.061 130M 
, \ ? Ce ; | C13 
Cas | and C”’ = AEE: : 
| C21 : C22 1 C3t | C33 
| | 
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The omitted subnetwork designated by C# in 1) above 
need not be considered since it reduces to the null set of 
elements by Step (B) if ci: and cu are not already zero. 
Let C’ and C’’ each represent a new C for further de- 
composition. 

Step (E): Repeat Steps (A) through (D) on each C 
until no further reduction can be accomplished in any 
step. 

Some observations concerning the decomposition 
procedure and its application are in order. 


1) A subnetwork with three or fewer vertices will 
always reduce to the null set of elements since a bridge 
network requires four or more vertices. 

2) Step (B) of the decomposition procedure is a 
special case of an iteration of Steps (D) and (A). In 
theory it may be removed from the procedure. It is 
useful, however, because it reduces the number of itera- 
tions required to reduce a network, in many cases. It 
also allows network classification by Property 3. 

3) Steps (A), (C), and (D) are all required steps for 
a complete decomposition procedure of this type. No 
one of them may be removed from the procedure with- 
out destroying the decomposition process. 

4) The location of the zero entries in the connection 
matrix governs the partitionings which are possible in 
the matrix procedure. 


V. Two-TERMINAL SWITCHING 
NETWORK SYNTHESIS 


The starting point for the synthesis of a two-terminal 
switching network is considered to be the minimum 
\V polynomial representation for the transmission func- 
tion, where “don’t care” conditions may be included 
in the determination of the minimum \/ polynomial. 

The procedure to be developed attempts to synthe- 
size a network which attains the lower bound stated in 
Theorem 1, or when this is not possible, to form a net- 
work with the smallest possible number of additional 
elements. The obviously uneconomical network repre- 
sentations thus are excluded from consideration in this 
procedure. 

In Part A of this section, the steps of the decom- 
position procedure are restated as conditions for decom- 
posing the transmission function into subfunctions with 
fewer literals. These steps also specify the network syn- 
thesis using connection matrices. This decomposition 
synthesizes the parts of the function which certainly 
are represented most simply as series or parallel sub- 
networks. 

In Part B the bridge condition is applied to the func- 
tions which are not decomposable to investigate their 
synthesis into simplest bridge networks. 

The use of redundant elements for both series-parallel 
elements and bridge elements is discussed in Part C. 


A. Series-Parallel Network Elements 
Step (A): The network may be synthesized with a 
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single series element incident to a terminal vertex if and 
only if the minimum \/ polynomial, f, may be fully 
assembled to give 


f = af’ 


where x is some literal that appears in each clause of 
f. The output connection matrix for f is 


The partially expanded connection matrix for f may 
then be represented as 


fe Oe i. One 
Ct On Vos Ori C = 0 ke als 
Bs rene | L Pee: | 
which corresponds to the literal « being incident to the 
input vertex or the output vertex respectively. 
Further synthesis introduces more vertices into the 
network (z.e., more rows and columns into the connec- 
tion matrix). Then, in the matrix, zeros are inserted as 
the new entries in the row and column corresponding 
to the vertex to which » is incident. Only the function 


f’ need be considered for further synthesis. This means 
that only the submatrix 


aa 


‘requires further expansion. Let f’ be a new f and repeat 


this step as long as it applies. 

Step (B): If in f there exists an alternation of clauses, 
B,, for which all the literals in each of these clauses 
appear in no other clause of f, then each of these clauses 
should be represented as a separate series path from 
the input vertex to the output vertex for a simplest 
network representation. Here, f= B,\/f’; where f’ repre- 
sents the alternation of the clauses of the minimum 
\/ polynomial for f not in By. 

More paths in the final network realization may also 
be removable by Step (B) of the analysis which are not 
evident at this stage of synthesis. The output con- 
nection matrix for f is 


c=[? ) -[ 1 paki 
fine Br \/ f 1 ; 


The B, part of the network may be totally synthesized 
now, and could easily be expressed in primitive matrix 
form, and f’ may be treated as a new f for further de- 
composition. 

Step (C): A necessary and sufficient condition for 
Step (C) to be usable for synthesis is that some repre- 
sentation of the transmission function be factorable 
into a conjunction of two \/ polynomials f’ and f’’. This 
gives 


f — (f’) (f'2) 
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It should be noted that the minimum \/ polynomial 
form for f may not obviously be factorable, and some 
form for f using redundant literals or clauses may be 
required to obtain a factoring. The method to obtain a 
factoring is not always evident immediately. Some 
formal method would be desirable. 

For f=/f'f’’, the connection matrix 


tal altars 


may be expanded as 


{0 -- 
CaO 1 ft. 
he oe tk 


Now f/f’ and f’’ may each be considered as a new f with 
the submatrices 


1 7 1 VES: 
c= | i and on = ae 
/ i! IMA 1 


Step (C) should be repeated as often as possible. In 
general, for a factoring of f as a conjunction of k factors, 


f = (g1)(g2) > + + (ge), 


and by permuting vertex 2 to the last row and column, 
the following matrix structure is obtained by repeated 


expansion: 
E Oy ORS en 0 
ris 1h 5k OO) 
al Oe thr eo O 
Dr bs ete ee 00 
@) . . il Lk 
Qua eae el 


Each g; may be considered as a new f for further syn- 
thesis. 

The resulting network from Step (C) is not generally 
a minimum network representation unless the factoring 
is a fully assembled representation of the minimum \/ 
polynomial. 

Step (D): For f in the minimum \/ polynomial repre- 
sentation, if the clauses may be split into two sets of 
clauses which have no literals in common Step (D) ap- 
plies. The function may be written as 


iy = ih Wate 
where f’ and f’’ are the alternations of clauses in the 
first and second sets respectively. The f’ and f’’ may 
each be considered as a new f and Step (D) may be re- 
peated until it no longer applies. Assume that the origi- 


nal f for this step is thus split into m separate functions— 
call them gi, ge, - - * , ga—such that 


f= eV g2V +°-+NMV Bm. 


Each g; may be considered as a new f and further syn- 
thesis may be carried out to obtain a primitive con- 
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nection matrix for each g;. Let these be called C; and 


partition each C;; j=1, 2,---, m, thus, 
Ai? | Aip? 
Gee ee ear 
| 
Ao? | Aas! 
| 


where 


1 = ai aig! + * + Cin, 
* - . e) 
Ay = | , Avi = [Ani]? = ae 
Qo? 1 a3’ + + * don,? 
and A»? is the remaining submatrix. 


The primitive connection matrix representing the 
network for f then is 


o eer | 
VO Aine AG end ost ul 


A12™ 
j=l 
ae oe Wn eae 
Ao"! A 9911 0 | 0 
C= patie ———|———— 
0 | 0 
Aoi 0 0 | Ao?” 
where 


V Ani 

j=1 
is the alternation of the A1’’s and 7; to jm 1s some ar- 
rangement of the indexes 1, 2,---+, m, picked at the 
discretion of the designer. 

Step (£): Repeat Steps (A) through (D) until no 
further simplification can be accomplished by any step 
of the procedure. 

Theorem 4: If the minimum \/ polynomial represen- 
tation, and fully assembled factors are used in all the 
steps of repeated application (A), (B), (C), and (D), and 
if the connection matrix reduces to a primitive connec- 
tion matrix, then the network representation is a simplest 
network representation, in series-parallel form. 

This theorem follows from Theorem 1 and the struc- 
ture of decomposition. 

The following two theorems show that the minimum 
VV polynomial representation is invariant under certain 
decompositions. 

Theorem 5: Given a minimum \ polynomial, M, 
then the alternation of any subset of clauses of M is also 
a minimum \/ polynomial. 

Theorem 6: For a given minimum \ polynomial, M, 
if M can be fully assembled into a conjunction of VV 
polynomials, then each of these \V polynomials is also 
a minimum \V polynomial. 

The following examples demonstrate the use of the 
decomposition procedure for synthesis. 


Le 


EE ee 


a. 
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Example 6: Let T be the transmission function to be 
synthesized where 


T = abefpgqrs \/ abghpgrs \/ cdefpqrs \/ cdghpgqrs 
V alpqr \/ mnpoar. 


T is in a minimum \/-polynomial. Applying Step (A) 
successively to p, g, and r gives the matrix 


ieee eg eS 
Vee Oren Oo 0 
Perel ve 2 FO r| 

C8 i DP og: 0 
Oe a le ae 
eigen 


where the row and column numbers indicate the num- 
bers assigned to the vertices, and J’ =abefs\/abghs 
V cdefs\/ cdghs\/ al\/mn. Applying Step (B) to T’ gives 


T’ = (al \/ mn) \/ (abefs \V abghs \/ cdefs \/ cdghs) 


HB, \ Tt. 
The connection matrix for J’ thus may be expanded as 
noe oO oy 
5 6 Pe ae Noe 
BEDI Pherks idem 
Stat |: B10 
eto 9h OD “A 


Step (A) may be applied then to 7’’ followed by Step 
(C), giving 
T” = s(ab VV cd)(ef VV gh). 


A final connection matrix which realizes the network in 
a simplest form is thus 


he Ag il aa ee ell an My po) 9 

Cee TONS CORT Oe OO. 0 0 

Da bin Ger O Oa st tellin | 8 0 

Seep) Oh, 1g So, OO: 0 

Pe OVOre OD <g- 4k PoP AR 4) 0 
‘Cle SO CO tO 7 1 am ~O ab\/ cd 

Ope we Oi er aed ie. ~ 0) 0 

A og te OR CN eee: 7/ogen 0 imag 0 0 
Be san: 70 47 OF noir OD 1 ef \V gh 

9 lo 0 0 0 ab\Vcd 0 O ef\V gh 1 


The graph for this network is shown in Fig. 6. 
Example 7: Let f=abc\/ abc. Applying Step (A) gives 
f=c(ab\/ab), and the connection matrix 
1 0 C 
GCG =;0 1 ab \/ ab 
Geab\/, db J 


Miller: Formal Analysis and Synthesis of Bilateral Switching Networks 


239 


Fig. 6—The network for Example 6. 


eee) CI 
aT *) 


(b) 
Fig. 7—The network graphs for Example 7. 


The corresponding network graph is shown in Fig. 7(a). 
This is a simplest network realization in series-parallel 
form. The bridge network seen in Fig. 7(b) is also a 
simplest network realization for the function. 

In general, the bridge condition should be applied 
at each step of synthesis. It is desirable, however, to 
obtain the most simple subfunctions possible with the 
decomposition process before applying the more com- 
plex bridge network synthesis process. When this is 
done, series-parallel network realizations are favored 
unless simpler bridge network realizations are possible. 


B. Bridge Network Subfunctions 


Consider network synthesis using a subcondition of 
the bridge condition given in Section III, where the 
expression f=A(BC\/DE)\/g is an assembling of the 
minimum \/ polynomial, and BD and CE are clauses 
in g. Let f=fsV gi, where fp matches A(BC\/ DE) \/ BD 
\/ CE and g matches gi: V/BD\V/ CE. 

If the literals appearing in A, B, C, D, and E do not 
appear in gi, then f, may be synthesized into a simplest 
network realization, where the connection matrix for f 
is 


Wk Sasa ao 
ieee) ete 
Coat 
ED eee 
De Ole ae 


and the network graph is shown in Fig. 8. An inter- 
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Fig. 8—A partially synthesized bridge network realization. 


change of rows and columns 3 and 4, respectively, gives 
an equivalent network. 

Now, take the case that some literals appearing in 
the subclauses A, B, C, D, and E also appear in gi, 
gi1~0, and also that f is a final subfunction of the de- 
composition procedure. In this case, C, represents a 
partially synthesized network. Further assembling of 
the clauses of g; may allow further synthesis. Three such 
assemblings are illustrated. 


1): If g: contains a clause with the subclause A, and 
one of the subclauses B, C, D, or E, then attempt to 
factor f as 


F=hV A(FC, V DiEy) V FD, \V CiEi V £2 
where 
Ci, Di, and F, are some three of the subclauses B, 
Cp and 
C,£, is either BD or CE, 


g2 1s the alternation of the remaining clauses, and 
the function A FC\\/ FD\\/ gz is matchable with gy. 


Under these conditions, the subclause element not 
represented as Ci, Di, or &; has the element F placed in 
parallel with it to further synthesize the network. For 
example, let B=B;, C=Gi, D=Di, and H=,. Then 


f= A[(BV FICV DE] V (BV F)DV CEV g.. 


The connection matrix and network graph are shown 
in Fig. 9(a). 
2): Suppose g; contains clauses of the form HBC and 
HDE. Then f may be assembled as 
J = f8 VB ON DERN 2s 
= (AV H)(BCV DE) \V BD\/ CE\ gs. 
The connection matrix and network graph are illus- 


trated in Fig. 9(b). Element 7 is parallel with element 
A. 


3): Assume that f can be assembled as 
f = A(BC\/ DE) \/ BD\ CE V K(LC\V ME) 
V LM \/ AK(BM \V/ DL) VY ga, 


where, if AK has the value zero, the factoring AK(BM 
\VDL) is omitted. Such a function may be represented 
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by the connection matrix and network graph in Fig. 
O(c). 

Application of these three steps with the general 
bridge condition allows synthesis of “extended simple - 
bridge networks” of the general form shown in Fig. 10. 
The lines between vertices indicate some functional 
interconnection. 

Example 8: A bridge inside a bridge. Let 


f = ac \/ bd \/ aduv \/ advx \/ aduxy \/ advwy 
V bcuw \/ bevx \/ beuxy \/ bcvwy 
= O(ad \/ be) \/ ac \/ bd; 


where 


O = uw \V/ vx V uxy V vwy. 


The partially expanded connection matrix may then be 
written as 


1 3 a 
F210. 2at on 
Qaix0 Cuad 
Cy = 
She |G Page wk Ga 
Abed. Oss 


The function Q also satisfies the bridge condition, where 
O = y(ux \V/ vw) VV uw \V/ ox. 
An expanded connection matrix for Q is 
gr ASS 
Sl ie ease 
C= Aes Os Wiaontiey 
DE {ae “word 
ON ey 


The primitive connection matrix for f then becomes 


1 Se ot ee 
1 F Oren ts OO 
Peat) cites Ye ark M0 
Sta" (one lean Omer eo) 

C= 

Aye GFW) A apy 
Sy |) (0) ge gi al 
0 LOR Org Te vee 


The corresponding network graph is shown in Fig. 11. 

Other extensions of the simple bridge result from 
factorings similar to those in 1), 2), and 3). 

A combination of subnetworks, which are formed by 
applying the bridge condition synthesis to more than 
one part of a function, is useful in forming more com- 
plex bridge networks. For this, the terms “vertex equal- 
ity” and “element equality” are defined. 

Definition 5: A vertex equality exists between two ver- 
tices in a network if the two vertices may be combined 
into a single vertex without altering the transmission 
function of the network. 
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Cy ly Bl aD! PCA 6 
Eiaal. PAd MIeK 
i Mt Oia 

(c) 


Fig. 9—The bridge networks formed by some additional 
synthesis steps. 
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Fig. 10—The extended simple bridge network. Fig. 11—A “bridge inside a bridge” network. 
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Definition 6: An element equality exists between two 
elements of common weight in a network if the elements 
may be combined into a single element of the same 
weight in a network without altering the transmission 
function of the network. 

An example now demonstrates how these equality 
relations may be used to aid in the network synthesis. 

Example 9: Consider the transmission function: 


eo AbL i ABCK/ADER \ BDK BCL 
WCE \/ EL: 


There are two possible ways to apply the bridge con- 
dition to f. The first is 


f=, V 1 = A(BCK V DEK) V BDK V CEK 
V ABEL \/ BDCIW/ EL, 
where 
fo, = A(BCK V DEK) V BDK \/ CER, 
and 
gi = ABL\/ BDCLY EL, 


which gives the subnetwork for fi, shown in Fig. 12(a). 
The second is 


f=fe V go = C(ABK V BDL) V ABL\V/ BDK 
WECEIG\Y EL N\/VADEK, 
where 
fe = CABK \/ BDL) \/ ABL\/ BDK, 
and 
g2 = CEK \/ EL\/ ADEK, 


which gives the subnetwork for f,, shown in Fig. 12(b). 

Some interconnection of these two subnetworks to 
represent f may be possible. The proposed method is to 
try some interconnection and then test its transmission 
function by means of the output determinant described 
in Theorem 2, Section II. For this example the inter- 
connection shown in Fig. 12(c) produces a simplest 


network realization. The vertex equalities assumed for 


this interconnection are 3=3’, 4=4’, and 5=5’, and 
element equalities are assumed for each pair of ele- 
ments with like weight. 

When applying the general bridge condition to net- 
work synthesis the BD and CE ciauses in the bridge 
condition may be proper subclauses of some clauses in 
the minimum \/ polynomial. As stated in the following 
theorem, this condition does not lead to a necessarily 
simplest network. Other assembling of the function may 
give better results. 

Theorem 7: If the BD and CE clauses in the bridge 
condition, A(BC\/DE)\/ BD\V CE, are actually proper 
subclauses in the minimum \/ polynomial representa- 
tion for f, and if ABD\/ ACE\/f#f, then A is not repre- 
sentable as a single bridge element with this method of 
assembling. 
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(a) 
5’ 
‘ D K 
fee GC We > 
A 1 
4’ 
(b) 
3 ae 


(c) 


Fig. 12—A network interconnection. 


C. Additional Element Requirements; Redundant Clauses 
and Literals 


This part considers modification of the minimum \/ 
polynomial, or its subfunctions, which aid in the syn- 
thesis of a simplified network realization of the function. 
Any such modification of the function must leave the 
value of the function invariant. 

The modifications can be classified into the following 
four types. 

Type 1: A zero value clause is joined to the function, 
where the literals in this clause appear in the functional 
representation being modified. 

Type 2: A clause, B, which is contained in the func- 
tion (2.e., B<f) is joined to the function. 

Type 3: A clause in f is replaced by some other clause 
such that the functional value of f remains unchanged. 

Type 4: The minimum \ polynomial representation 
is replaced with some other minimum \ polynomial 
representation of the function. 

Although Type 2 includes the Type 1 modification, 
the Type 1 is listed separately since it is quite useful in 
practice. 

The following examples demonstrate the use of modi- 
fications on the functional forms to aid network synthe- 
sis. 

Example 10: Let 


f= abV ad\V cd\ bee. 
Type 2 = ab\ ade\/ ad \/ cd \/ bde 
e(bc \V ad) \/ ab \/ cd \/ ad = fy \/ ad. 


1958 


Fig. 13—A bridge network resulting from using redundancies. 


kes aunt 


Fig. 14—A series-parallel network formed by 
using redundancies. 


Types 1 and3 = e(bc V ad) V b(bc VV ad) VV ab \ cd 
(6 \/ e)(bc \/ ad) \/ ab V cd. 


I 


The resulting network graph is shown in Fig. 13. 
Example 11: Let 


f = axV bex \/ ab = x(a V bc) V abd. 


Type1 = x«(aV bc) V b(aV bc) = (bV x)(aV de). 


The resulting network graph, using Step (C) of the 
decomposition procedure, is shown in Fig. 14. 

These examples illustrate that redundant clauses and 
literals are useful in network synthesis. 

There usually are many ways to modify a given trans- 
mission function. One should attempt first to modify 
the function by using only literals which already appear 
in the minimum \/ polynomial representation, since 
such changes may lead to a network synthesis which 
gives a simplest network as described in Theorem 1, 
Section II. For such redundancies, Step (C) of the de- 
composition procedure, and the bridge condition may 
lead to further synthesis. 

The form of the bridge condition indicates the re- 
-dundancies which might apply to satisfy the bridge 
condition. Three possible forms are given now. 

Form 1: The function contains clauses of the type 


f = A(BC V —) V BHx V CK. 


The bridge condition factoring may be completed by 
joining the zero value clause AHxK& to f, giving 


FAL BON SP AAKS) \/ BH ENVCKS, 


which results in the network graph shown in Fig. 15(a). 
Form 2: The function contains clauses of the type 


f = Ax(BC\V—) V Bt V CH. 
Joining the zero value clause Ax#H to f gives 
f = Ax(BC \V 4H) V Be V CH. 


A resulting bridge network realization is seen in Fig. 
15(b). 
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Se dae 
a 


(a) 
oe 


Baorss c 
a al 


aN 
6 


(c) 


Fig. 15—Several bridge network forms resulting 
from redundancies. 


Form 3: Let 
f= Bt VV De BON CD 
=a BDO) BCS CD: 


By the Type 3 modification, replace Bx with BCx and 
Dx with DCx. Then, 


Sem BC DG) NT BC NCD, 


A resulting bridge network realization is shown in Fig. 
15(c). Forms 1 and 2 use the Type 1 modification. These 
redundant clauses correspond to “blocked paths” in the 
network. This blocking serves to prevent “sneak paths.” 
Form 3 uses no “blocked paths.” 

The following extension of Theorem 6 can be shown 
to hold. 

Theorem 8: If a minimum \/ polynomial representa- 
tion of f may be altered using redundancies in such a 
way that f may be fully assembled as 


Meet 4 en 


where g; is an \/ polynomial, and no fully assembled 
representation exists without using all of the redund- 
ancies applied, then all the g,;’s are minimum \ poly- 
nomials. 
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If no further network synthesis can be accomplished 
using only redundancies which appear in the minimum 
\/ polynomial representation, M, then redundant clauses 
using literals which do not appear in M may be tried. 
An upper bound on the number of different literals that 
should be used simultaneously is the difference between 
the number of elements required for the simplest known 
series-parallel network realization and the number of 
different literals appearing in M. The synthesis proce- 
dures using such redundancies are the same as previously 
described. 

Even though the synthesis procedure may lead to a 
simplest network realization, this in no way implies 
that this is a unique simplest network realization. In 
fact, the process was shown to yield a simplest realiza- 
tion only under very restricted conditions. 

The three main advantages claimed for the procedure 
are first, a complex transmission function may be split 
into a set of less complex functions which can be realized 
more simply than the original complex function; second, 
the necessarily series and parallel parts of the network 
are recognized as such in the transmission function, and 
third, some simplest bridge network configurations are 
recognized in the transmission function with the aid of 
the bridge condition. 
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VI. CONCLUSIONS 


The analysis of two-terminal networks is carried out 
by a decomposition procedure described in Section IV. 
This decomposition procedure is also restated and used 
for network synthesis for the necessarily series or paral- 
lel parts of a network as determined from the transmis- 
sion function for the network. A bridge condition also is 
developed to aid in the synthesis of two-terminal bridge 
networks. The minimum \/ polynomial is used as a 
starting functional representation for synthesis. In ad- 
dition, the use of redundant variables is considered as 
an aid to synthesis: 

Under certain conditions the synthesis procedure is 
shown to yield network realizations with the fewest pos- 
sible number of elements. 

Some extensions to multiterminal networks are pos- 
sible. More complex conditions seem to be required, 
however, before adequate synthesis procedures can be 
described to give minimal network realizations for 
multiterminal networks. A counting and classification 
procedure for all bridge network types with n vertices, 
for small ”, would also find ready application to aid in 
the formalization of bridge network synthesis. 


1s Miller, op. cit., pp. 71-78. 


A Transistor Pulse Generator for Digital Systems* 
DOUGLAS J. HAMILTON+ 


Summary—A design procedure is developed for a new transistor 
pulse generator circuit suitable for use as a building block in a digital 
system. The circuit produces a pulse whose shape is relatively inde- 
pendent of variations in transistor parameters and load current. 
Pulse durations in the range 0.5 microsecond to 20 microseconds 
and load currents of several hundred milliamperes may be obtained. 


INTRODUCTION 


HE requirements placed upon a pulse generator 
which is to be used as a building block in a digital 


computer or other pulse system may be general- 
ized as follows: 


1) The circuit is to be triggered by a waveform of 
somewhat arbitrary shape and repetition rate. 

2) The circuit must provide an output pulse whose 
amplitude and shape are relatively independent of 


* Manuscript received by the PGEC, May 1, 1958; revised 
manuscript received, June 18, 1958. 
+ Stanford Electronics Labs., Stanford, Calif. 


the input waveform, and are also independent of 
load variations over a considerable range. 

3) Because the circuit is to be used repetitively 
throughout the system, the output pulse must be 
reasonably independent of variations in transistor 
parameters. 


The blocking oscillator circuit appears attractive for 
such an application since it produces pulses whose shape 
is relatively independent of the input waveform. How- 
ever, the circuit has several disadvantages which must 
be overcome if it is to be useful from a system point of 
view. From some fundamental considerations of tran- 
sistor blocking oscillators a design procedure will be 
evolved for a new circuit which is capable of delivering 
pulses at repetition rates up to 500 ke, with durations in 
the range of 0.5 microsecond to over 20 microseconds 
having maximum variations in pulse duration of +10 per 
cent due to variations of transistor alpha from 0.95 
to 1.0 and load current from 20 ma to 200 ma. 
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SOME FUNDAMENTAL CONSIDERATIONS OF TRAN- 
SISTOR BLOCKING OSCILLATORS 


Two basic blocking oscillators circuit configurations, 
without bias or damping, are shown in Fig. 1. The cir- 
cuit of Fig. 1(a) uses collector-to-base feedback and is 
referred to as “base control,” while the circuit of Fig. 
1(b) uses collector-to-emitter feedback and is referred 
to as “emitter control.” A third configuration using emit- 
ter-to-base feedback is only slightly different from base 
control and is not discussed. For the purpose of analyz- 
ing the circuits the transformers are assumed to be ideal 
with a turns ratio 2 and a primary inductance L. 

As far as rise-times are concerned, there is no essential 
difference between the two circuits. For a given natural 
frequency corresponding to a growing transient in the 
case of emitter control it can be shown that there is 
another transformer turns ratio which will produce the 
same natural frequency in the case of base control.! It 
will, however, become apparent that emitter control 
is to be preferred where variations in pulse duration 
must be kept to a minimum. 

The model? which is used to represent the transistor 
in discussing the pulse duration has been chosen be- 
cause it retains only the properties which are essential 
to the understanding of the operation of the circuit; it 
is shown, together with its volt-ampere characteristic, 
in Fig. 2. For switching applications the model is a 
reasonable representation of an alloy junction transistor 
and has the following properties: 


1) The impedance between any two terminals is high 
when both junctions are reverse biased. 

2) When both junctions are forward biased the volt- 
age between any of its terminals is zero. 

3) When both junctions are forward biased the im- 
pedance between any of its terminals is zero. 

4) When the ratio of collector current to emitter cur- 
rent exceeds a value A, the collector junction be- 
comes reverse biased and the impedance between 
collector and base and between collector and emit- 
ter is high. 


It will be recognized that A is the dc alpha of the 
transistor. If base current is the independent variable, 
as it is in the case of base control, application of Kirch- 
hoff’s Current Law to the model shows that the col- 
lector junction will be reverse biased when 


Es A 
—> . 
Tan Aes A 


With the models for the circuit and the transistor 
thus defined, one is in a position to study the pulse 
duration. If in some manner the emitter junction is for- 
ward biased, regeneration will occur and will continue 


1J. G. Linvill, “Transistor electronics,” class notes for course 


EE227, Stanford University, Stanford, Calif.; 1957. 

2 This is similar to a model proposed by R. B. Adler, “A large 
signal equivalent circuit for transistor static characteristics,” Res. 
Lab. of Electronics, Mass. Inst. Tech., Cambridge; August 30, 1951. 
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Systems 


(a) (b) 


Fig. 1—Two basic blocking oscillator circuits. (a) Base control. 
(b) Emitter control. 


Ve 


Te = CONSTANT 


I, 
I.=Ale=A I, 
I-A 


Fig. 2—Transistor model and its V-I characteristic. 


until the collector junction becomes forward biased; 
it will remain forward biased until the primary induct- 
ance of the transformer causes the collector current to 
increase to a value [,=AJ,., at which time regeneration 
again occurs, continuing until both junctions are reverse 
biased. The pulse duration 7 will be defined as the 
amount of time that the collector junction is forward 
biased. 

Consider first the case of emitter control. For the 
duration of the pulse the collector junction is forward 
biased and the voltage across the primary winding of the 
transformer remains constant at £; hence the control 
current J,= E/nR, flowing in the emitter circuit remains 
constant. The collector current, however, is made up of 
two components; one is J,/n and the other is Jz, the 
current in the primary inductance of the transformer. 
Iz is initially zero but increases linearly with the time: 
I,=Et/L. It is sometimes convenient to visualize this 
process in terms of required collector current I,, and 
maximum available collector current Igmax,? as shown 
in Fig. 3. When the required collector current, which is 
I./n+Et/L, exceeds the maximum available collector 


3D. J. Hamilton, “Pulse Duration of Transistor Blocking Oscil- 
lators,” Hughes Aircraft Co., Culver City, Calif., Rep. M-122; 1956. 
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particular types of loading are shown in Fig. 4. If the 
load resistance is connected in shunt with the primary 
winding of the transformer as in Fig. 4(a), the required 


collector current is increased by an amount equal to the © 


load current and a correspondingly shorter time is re- 
quired for the collector current to increase to [,=ATe. 
If the load resistance is connected in series with the 
primary winding the voltage across the primary winding 
does not remain constant during the pulse, the current 
in the primary inductance does not increase linearly 
with time, and the emitter current does not remain con- 
stant. The resulting effect upon pulse duration is de- 


Fig. 4—Sketches showing the effect upon pulse duration of two types of loading* 
(a) Rx shunting primary. (b) Rz in series with primary. 


current, which is AJ,, the collector junction becomes 
reverse biased and the pulse ends. The time thus elapsed 
is the pulse duration and is easily calculated to be 


L 


NIN 


(A — 1/n). (1) 


T= 


Using similar reasoning for the case of base control 


one finds that 
If ( A ~) (2) 
eee Nt ee 


Eqs. (1) and (2) show that to minimize variations in 
pulse duration due to changes in A, the transistor dc 
alpha, emitter control should be used in preference to 
base control since variations in A will only be of the 
order of 5 per cent, while 4/1—A may change by 5 to 1 
in a normally encountered distribution of transistors. 
Thus, further discussion is confined to emitter control. 

There are several ways of connecting a load to the 
circuit, but in all cases a variation in load resistance 
produces-a variation in pulse duration which would 
generally be intolerable ina system. It has been shown 
that loading the circuit also deteriorates the rise-time 
and the trigger sensitivity.* 

Sketches of the effects upon pulse duration of two 


4 J. G. Linvill and R. H. Mattson, “Junction transistor blocking 
oscillators,” Proc. IRE, vol. 43, pp. 1632-1659; November, 1955. 


picted in Fig. 4(b). The effects of other types of loading 
may be similarly deduced. 

It is thus apparent that if the circuit is to be useful 
in a system whose requirements are those previously 
stated, modifications must be made to couple the load 
to the circuit in such a way that the immittance pre- 
sented to the circuit remains relatively constant with 
varying load. 


LIMITATIONS OF THE MODEL 


The property of the model which led to an easy cal- 
culation of the pulse duration is that with both junc- 
tions forward biased the voltage between any of the 


‘terminals is zero until J.=AT,. In reality this is in- 


correct because of junction potentials and the presence 
of 7’, the base spreading resistance. The latter is the 
major contributor to the voltage appearing between the 
base and collector and the base and emitter. If the model 
is modified to include rp’, analysis reveals that the pulse 


duration is given by 
1 
wie 
n 


£ 
Ss in lide 3 
¥ * MD [Re bitpih he mal 8) 


where 


1958 


Hamilton: A Transistor 


As an example consider a case in which 


L = 4.17 mh. 
n= 2 

A = 0.98 
R. = 100 ohms 
ry’ = 100 ohms 


Eq. (3) yields 7.’=9.4 microseconds. If 7,’ is assumed 
to be zero and (1) is used, one obtains 7,-=10 micro- 
seconds. In general, if R.>>r,’(1—A), and if accuracy 
of the order of 5 per cent is acceptable, (1) may be used 
in preference to (3). 

Clearly, to make the pulse duration independent of 
transistor parameter variations one should make R, 
>r'(1—A), where 7’ is the largest value of 7,’ that 


will be encountered ‘and A is the smallest value of A 


that will be encountered. This has the effect of maintain- 
ing the drop across R., and hence the emitter current, 
constant during the pulse. 

Several 8 microsecond blocking oscillators with R, 
=270 ohms and Raytheon 2N113 transistors were 
constructed and in all cases the measured pulse duration 
was within 2 per cent of the value calculated by (1). 


A NEw CiIrRcuIT 


Two problems must be solved in order to produce 
from the blocking oscillator a circuit useful in a system. 
Bias circuitry must be provided so that the circuit is 
energized only by a triggering waveform, and some 
means must be devised to couple the load to the circuit 
so that variations in load will not alter the pulse duration 
and rise time. These problems are solved by a new cir- 
cuit shown in Fig. 5. 

In this circuit 7; functions as a blocking oscillator, T» 
is used to couple the load to the circuit, and D, per- 
forms a dual function as a mechanism for providing 
cutoff bias and reverse base current to T>. 

The base-emitter circuit of T> is in the emitter loop of 
the blocking oscillator and all of the current flowing in 
the emitter of 7; flows in the base of 7». This current 
is made large enough to insure that the collector of T 
is forward biased during the pulse under all conditions 
of loading, which insures that the amplitude will be 
constant. This method of coupling the load has several 
advantages. First, a large current is available in the 
collector circuit of T2; second, the variation with load of 
the impedance between emitter and base of 7» will be 
small and can be masked by R, so that the pulse dura- 
tion is not affected; and finally, the impedance between 
emitter and base of J» is small enough to permit good 
rise-times to be obtained. A capacitor is placed in shunt 
with R, to eliminate the effects of R, during the rise 
time. 

Proper bias voltages to insure that both transistors 
are cut off until the application of a trigger are pro- 


Pulse Generator for Digital Systems 


247 


im N TRIGGER 
Fig. 5—New circuit which effects the desired improvements. 


vided by diodes D; and Dz, the former germanium, and 
the latter silicon. Both diodes are forward biased; D2 by 
(Iz—1;) and D,; by 1. Thus the emitter of 7; is reverse 
biased by an amount equal to the voltage across D,, and 
the emitter of 72 is reverse biased by an amount equal 
to the difference in voltages across D; and Ds, normally 
about 0.4 volt. When a trigger is applied D; is reverse 
biased, 7; turns on causing 72 to turn on and regener- 
ation occurs. 

It is the second function performed by D; which makes 
the circuit practical. During the pulse the collector of 
T, is forward biased and a large number of minority 
carriers are stored in the base region. Assume for the 
moment that D; is not present. At the end of the pulse 
the voltage across the transformer reverses and the emit- 
ter of 7, is driven negative. The base current of T> is 
reduced to zero and the minority carriers stored in 
the base are removed only by recombination. Because 
of the time required by this process, the duration of the 
pulse appearing across the load may be considerably 
longer than the pulse duration of the blocking oscilla- 
tor. With D, in the circuit the emitter of 7} is restricted 
to a reverse bias of only a few tenths of a volt. When the 
voltage across the transformer reverses, the base of JT» 
is driven positive and reverse current flows, sweeping 
out the minority carriers and reducing the storage time 
to a negligible value. Hence the duration of the pulse 
across the load is the same as the pulse duration of the 
blocking oscillator. Because all of the energy stored in 
the transformer during the pulse may not be used in 
sweeping out the minority carriers, additional damping is 
provided by Ds; and Rs. 

The ratio of reverse base current to forward base 
current is easily found. The forward base current Ip; is 
equal to J,, the emitter current of 7;. During the pulse 
the current built up in the primary inductance is 
I.(A—1/n). When the pulse ends this current must 
continue to flow; it appears in the secondary winding 
multiplied by ~ and in a direction opposite that of J,. 
Hence, J;,, reverse base current, is [.(nA —1). This as- 
sumes that R3/n? is much larger than the impedance be- 
tween emitter and base of JT». The ratio of reverse to for- 
ward base current at the end of the pulse is therefore 
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Hix 


= (nA — 1). 
of 


In calculating the pulse duration (1) must be modi- 
fied to include the effect of bias voltage and the base- 
emitter circuit of T>: 


IE 1 E- NV ve a= Ls) 
ee oe 
nR, nN E+ E, 


where V3, E;, and E are the magnitudes of the ferward 
base-emitter voltage of 72, the forward voltage of the 
silicon diode, and the supply voltage, respectively. 


DESIGN PROCEDURE 


A design procedure may now be established to be 
used in adapting the circuit to a particular system: 


1) Choose a value of base current Jy; for JT. large 
enough to insure that the collector will be forward 
biased during the pulse under all load conditions. 

2) Select R, so that with this current flowing the drop 
across R, will be large compared to variations in 
the emitter-base voltage of Z2 which are an- 
ticipated under all conditions of load and tran- 
sistor parameter variations. This amounts to mak- 
ing R, large enough to mask the variations in im- 
pedance between emitter and base of 72 so that 
I., and hence the pulse duration, remains constant. 

3) Select a value for C. The time constant of the emit- 
ter circuit of the blocking oscillator is approxi- 
mately 

PROT 


Rone 


where 7,’ is the base resistance of I>. C should be 
chosen so that this is about one fourth the pulse 
duration. Thus the pulse duration will not be 
influenced by the charging of the capacitor. 

4) Choose a turns ratio n. Here one wishes to satisfy 
two criteria. When and R, have been selected a 
supply voltage is automatically determined: E 
—nI,;R,.. It is necessary that 2 be chosen so that 
the maximum collector voltage is not exceeded, 
particularly at the end of the pulse when the 
transformer voltage reverses. Within this restric- 
tion one wishes to select, if possible, a value of 
n to obtain the best rise time. The selection of an 
optimum 7 will be discussed in the next section. 

5) Using (4) and the nominal value of Vz., select a 
value of primary inductance to produce the de- 
sired pulse duration. 


The 0.5 microsecond pulse generator shown in Fig. 6 
was designed to supply a load current of 20 ma to 200 
ma. Two transistors are used in the output to minimize 
the effects of a decrease in transistor current gain at high 
emitter currents. The use of two transistors has the 
additional advantage of reducing the parasitic impe- 
dance in the emitter loop of the blocking oscillator. As 
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Fig. 6—Pulse generator for 0.5 microsecond, 250 ke operation. 


can be seen from the waveforms in Fig. 7, the variation 
in pulse duration over the entire load range is less than 
10 per cent. The circuit has been operated at pulse 
repetition frequencies as high as 500 kc. 

In the tested circuits a 4-volt trigger was required for 
consistent operation, and no change in output rise-time 
was observed for triggers in excess of 8 volts. 


CHOOSING A TURNS RATIO 


The problem of choosing a turns ratio for the pulse 
generator which is to be used in a system can generally 
be reduced to the following terms: 


Given the maximum 7,’ and collector capacitance C,, 
and the minimum small-signal current gain a@ and 
radian cutoff frequency wo, which are to be expected 
from the transistors which will be available, approxi- 
mately what turns ratio should be used to obtain the 
minimum rise-time? 


Using the method of Linvill and Mattson‘ an analysis 
will be given here for the circuit of Fig. 5 which ex- 
presses the optimum turns ratio in terms of the prod- 
uct 1/C.wo. 

The circuit of Fig. 5 has certain natural frequencies 
which correspond to growing transients, and the regen- 
erative process may be thought of as a multiplication by 
these transients of the initial conditions established by 
a trigger source or other means. It is generally not neces- 
sary to calculate the actual rise time; in fact, to do so 
would be difficult because of the nonlinearities present in 
the regenerative process. One is interested in minimizing 
the rise time; this implies maximizing the growing 
transients. A simple model will be used which readily 
lends itself to the maximizing of the growing transients. 

If a transformer with good coupling is used, such as 
one with a ferramic cup-core, the leakage inductance 
may be neglected; the presence of leakage inductance 
changes the actual rise time but does not materially 
alter the optimum turns ratio.® During the transient 


5 J. G. Linvill, private communication. 
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Vertical: 5v/cm 


Horizontal: 1usec/cm 


Horizontal: 0.1usec/cm 


Load: 512 


(Trigger—4 volts, 0.2usec rise 
time) 


Horizontal: 0.1usec/cm 


Load: 510 


Fig. 7—Waveforms for the circuit of Fig. 6. 


period the base-emitter impedance of J» consists pri- 
marily of 7,’, and the base-emitter circuit will be re- 
placed by 7,’ in the model. The capacitor C will usually 
be sufficiently large that its impedance may be neglected; 
this can later be verified by calculating its impedance at 
the natural frequency and comparing it with 7’. The 
frequency dependence of alpha may be approximated by 


ao 


Ss 
1+— 


WO 


where s is the complex frequency variable; and the 
model is then as shown in Fig. 8. 

If one sets the determinant of the model equal to zero 
and solves for the natural frequency S corresponding 
to a growing transient, the following is obtained: 


1 ( NAO 1) 
SNES 
nC.= ee oe) 
in — in + i 


It can be seen from (5) that for a in the range 0.95— 
1.0 the actual value of a will not have much influence 
upon S; therefore a further simplification can be made 
by assuming that a)=1. It is also more convenient to 
maximize S/w, the ratio of the natural frequency to the 
radian cutoff frequency. Letting p=S/w» one obtains 


~( nN 1) 
BuNise Pp 


2n? — 2n +1 


/ = 
ry Coo = 
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Fig. 8—Model of the circuit during the transient period. 


from which 


at 4/14 4n 
jin (Qn ane eee 


The optimum value of 2 to maximize p is given by 


= 2+ —4/1+- 1+ kK . 
a cage 3 =( K ) (6) 


It will be noted that K is the product of the three tran- 
sistor parameters of interest: 75’C.wo. 

As an example, suppose that the available transistors 
Rave fi max = 100 ohms, Conax= 0:2 X10-™ farad, ena 
=62.8X10® radians per second (f=10 mc). The opti- 
mum turns ratio 1s Mope=3.55, and pmax =1 (the natural 
frequency is 10 mc). If C=0.001 X10~* farad, the capaci- 
tive reactance at this frequency is 15 ohms. 

It should be emphasized that the model used here is 
highly idealized and does not take into account the 
nonlinearities in emitter resistance and collector capaci- 
tance; the results are to be regarded as indicative only 
of the approximate turns ratio. It is also to be empha- 
sized that the criterion that the maximum collector 
voltage of the transistor should not be exceeded must 
always be satisfied in choosing the turns ratio. 


p= 


CONCLUSIONS 


The analysis in this paper has been aimed at present- 
ing a simple picture which illuminates the operation 
and some of the shortcomings of the blocking oscillator. 
The analysis has led to a new circuit, suitable for use as 
a building block in a digital system, which eliminates 
the most troublesome of these shortcomings without 
resorting to undue complexity. The design procedure 
presented has been useful in adapting the circuit to 
several different systems. While the circuit manifestly 
cannot qualify as a precision pulse generator, the varia- 
tions in pulse shape with load and transistor parameters 
are sufficiently small to permit its application in many 
systems whose requirements are not unduly stringent. 
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CORRECTIONS 


Douglas B. Netherwood, author of “Logical Machine 
Design: A Selected Bibliography,” which appeared on 
pages 155-178 in the June, 1958, issue of these TRANS- 
ACTIONS, has called the following corrections to the 
attention of the editor. 

On page 155, column 2, line 4, [23] should be read for 
[24]. On page 156, under Reading Suggestions, the 
references should be 29, 50, 61, 105, 143, 192, 234, 240, 
262, 298, 361, 381, and 456. In the Bibliography itself, 
cross-references shown as “(400)” should read “[397]” 
for the following entries: [23], [108], [221], [283] 
[286], [440], [442], and [449]. 


D. W. Davies, author of “Switching Functions of 
Three Variables,” which appeared on pages 265-275 of 
the December, 1957, issue of these TRANSACTIONS, has 
requested the editor to publish the following corrected 
figure. 


US fe Sti. Be EGS S Fle WH 2 es 


Fig. 5—Functions of classes d and e. 
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Correspondence 


Switching Circuits as Topological 
Models in Discrete Probability 
Theory* 


The purpose of this note is to call atten- 
tion to a substitution rule which can facilitate 
discussion in the related areas of information 
theory, switching theory, computers, and 
probability theory. The discrete probability 
theory and its notation is usually confusing 
to the beginner. It is believed that the 
scheme to be suggested will facilitate both 
presentation and understanding of the sub- 
ject. 

The central notion in the probability the- 
ory is that of an event Z with a certain prob- 
ability of occurrence Q(£) or nonoccurrence 
P(£)=1—Q(£). The event E can be as- 
signed a switching variable e which takes the 
value 1 when the event occurs and 0 when 
it does not occur. A joint event E involving 
two simple events £; and EF: can be related 
to the simple events in the most general way 
by writing the logical connection between 
these events in the canonic form of a switch- 
ing function of two variables! 


e = foer’es’ + frer’eo + foerer’ + frees (1) 


where the coefficients f; can be either zero or 
one in the switching function. The primed 
variables are logical negatives of the un- 
primed variables, thus e:/=1 is interpreted 
to mean that event £; did not occur. The 
sample space and events in that space may 
be interpreted effectively in terms of the 
topological model given for switching func- 
tions.? Simple and joint events correspond to 
collections of cells of various orders. 

Because of the correspondence between 
switching functions and switching circuits, 
these circuits may serve as_ topological 
models of an experiment or system. As an 
example of such a model, the joint event Z 
which occurs if and only if EZ; and Ez occur 
is represented by 


é€ = 1&2 


which, in turn, is represented by a switching 

- circuit consisting, for example, of two in- 
dependent relay contacts in series. In this 
representation, unity corresponds to closure 
and zero corresponds to an open contact or 
open network. Similarly an event £ which oc- 
curs if either or both of two events £; and Ez 
occurs is represented by 


e=ate 


which, in turn, corresponds to the parallel 
connection of two independent relay con- 
tacts. 

The fundamental relation which links 
the switching circuit model with the discrete 
probability theory is the following substi- 
tution rule which is offered without proof: 


* Received by the PGEC, March 31, 1958. 

1M. Phister, “Logical Design of Digital Com- 
puters,” John Wiley and Sons, Inc,. New York, N. Y., 
pp. 50-53; 1958. 4 

2R.H. Urbano and R. K. Mueller, “A topological 
method for the determination of minimal forms of a 
boolean function,” JRE TRANS. ON ELECTRONIC 
CompuTErs, vol. EC-5, pp. 126-131; September, 1956. 
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The probability Q(£) of occurrence of a 
complex event which can be represented in 
terms of a switching function of 1 variables 
in the complete product form 


M 
e= Li fbi (M=2"—1) (2) 


where ~; is a complete product, can be 
written from (2) by the following substitu- 
tions: 

1) An unprimed variable is replaced by 
the probability Q of occurrence of the asso- 
ciated event. 

2) A primed variable is replaced by the 
probability P=1—Q of non-occurrence of 
the associated event. 

3) The Boolean + and : go over into 
identical symbols in ordinary algebra. 

As an example of the application of 
the substitution rule, consider the following 
experiment. Five numbered dice are thrown 
and a payoff # occurs if sixes appear on one 
or more of the following combinations of the 
dice 


(1,2) (1,4,5) (3,4) (2, 3, 5) 
What is the probability of a payoff? 
The switching function is 


€ = e1€2 + e1€4€5 + C34 + Cr€3@5. BEF (3) 


Two possible switching circuits which yield 
(3) are shown in Fig. 1. In the series-parallel 
network the branches are not independent, 
while in the bridge network the branches are 
independent. Whether the branches are in- 
dependent or not, the substitution rule 
applies. 


él 2 
é3 4 
éy 4 @5 
é2 é3 &5 
O O 
ey €3 
© 
O 
€2 4 


Fig. 1—Switching circuit models of experiment. 


The canonic expansion of (3) yields 


’ 
E= €1€nl34l5+ Cr€nlse4es + C1€2esl4 Ct €1ln€s€4 C5 
+e, e2€3 "es@5+ €1€2€3' C45" + e)e2€3 "es! €5/ +€16263'€4' es 
+ €1€2! €3€4€5 + €' 2€3€4€5-+ €1/ €2! esses ter €2/ 63’ €4e5 


+ 61 €2€3€4/ 5+ €1/ €2€364€5" a. €1/ €2! 36405" ah €1€2 esses! 
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and the substitution rule gives 


Qe) = Q102030405+01020:04Ps-+010203P40s 
+010203P4P5+0102P30s05+Q0102P30iPs 
+0102P3P4Ps+0102P3Ps0s+01P20:0:05 
+P102030s05+ PiP2030105+ O1P2P30105 
+P1Q203PsOs+P1020:04P5+PiP20:0sPs 
+01P2Q:04Ps. 


(Since each of the Q’s is 2 and each of the 
P’s is % in this example, the numerical result 
is that the probability of a pay off is 
476/7776.) 

In applying the switching circuit repre- 
sentation to the law of conditional probabil- 
ity, four possible probabilities will be con- 
sidered: 


Q(A|B) which is the probability that A =1, 
given that B=1 

Q(A|B’) which is the probability that A =1, 
given that B’=1 (B=0) 

P(A|B) which is the probability that A =0, 
given that B=1 

P(A|B?’) which is the probability that A =0, 
given that B’=1 (B=0). 


The last two may be written in terms of oc- 
currence probabilities Q: 


P(A | B) = Q(4’|B) 
P(A| B’) = Q(4’| B’). 
Hence only the occurrence probabilities are 


considered. The law of conditional prob- 
ability states that for any x and y 


Q(«| ¥) = Oley) /O(y) 


and in this law xy may be interpreted as a 
circuit consisting of networks representing 
x and y connected in series. Thus the prob- 
ability of occurrence of x, when y is known 
to have occurred, is the probability of clos- 
ure of the network consisting of x and y in 
series, divided by the probability of closure 
of the network representing y. Thus for a 
conditional probability situation, there are 
two networks to which the substitution rule 
is applied. 

As an example, each of the conditional 
probabilities is examined with 


A=até 
[cs 2 ote 3. 
Hence 


AB = e + e163 = €1€2€3 + €1€'€3 + €1'€0€3 
+ er€2¢s' + e1'e263" 
AB’ = €1€2/ €3" 
A’B => €1' es! €3 
ABE = €1'€2/ €3" 


A = ej€2€3 + e1€0'€3 + €1'€2€3 + e1€2€3’ + e1/€2 €3 


ste €1€2'€3" 

B = ey€2€3 + e1€2'€3 + e1/€2€3 + eeres’ + e1' e263" 
+ e1'e0'e3 

B! = e62'¢3' + e1'¢0'e3'. 


Then 
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Q10203 + QiP203 + PiQ203s + Qi02P3 + PiQ2Ps 


Q(A | B) = Q(AB)/Q(B) = 


= QiQ3 + PiQ2 + Qi02P3 
Qs —- QoP3 


Q10203 + Q:P203 + PiQ203 + Qi02P3 + PiQePs + PiP20s 


, yp y OiP2Ps ca 
Q(A | B’) = Q(4B)/O(B’) = O.P.P; + P:P»Py = Q1 
Q(4'| B) = O(4'B)/0(8) = =e 
1—QiP2P3 — PiP2P; 1—P2Ps; Qi + QoP3 
PyPoP3 


Q(A’| BY) = Q(4'B)/O(B)) = 


Contributors 


Harvey L. Garner (S’51—A’51-M’57) 
was born in Lake, Colo., on December 23, 
1926. He received the B.S. and M.S. degrees 
from the University of Denver, Denver, 
Colo., and the Ph.D. degree from the Uni- 
versity of Michigan, Ann Arbor, Mich. 

From 1949 to 1951, he participated in 
the Cosmic Ray Research program at the 
University of Denver and the Inter-Uni- 
versity High Altitude Research Laboratory, 
Echo Lake, Colo. 

He was associated with the development 
and operation of the MIDSAC and MIDAC 
computers at the Engineering Research In- 
stitute of the University of Michigan during 
1951-1955. 

In September, 1955, he became an in- 
structor in the Department of Electrical 
Engineering at the University of Michigan. 
In this capacity he has been active as a 
co-supervisor of the MIC computer project 
and has been in charge of several special in- 
tensive summer computer courses. He is 
presently assistant professor of electrical 
engineering. 

Dr. Garner is a member of Sigma Xi, 
Sigma Pi Sigma, and the Association for 
Computing Machinery. 


2 
“e 


Roderick Gould (S’55—M’58) was born in 
Shanghai, China, on December 25, 1932. He 
received the B.A. degree in mathematics 
from Amherst College, Amherst, Mass., in 
1954, and the S.M. and Ph.D. degrees in ap- 
plied mathematics from Harvard Univer- 
sity, Cambridge, Mass., in 1955 and 1957, 
respectively. His graduate work was in the 
design and application of computing sys- 
tems, and his doctoral dissertation dealt 
with the application of graph theory to con- 
tact network synthesis. 

During 1957 and 1958, he was with the 
Comité d’Etude et d’Exploitation des Calcu- 
lateurs Electroniques, Brussels, Belgium. 
Dr. Gould died earlier this year. 

Dr. Gould was a member of the Associa- 


OPP, Pi Pabs a 


l- 


tion for Computing Machinery, the Society 
for Industrial and Applied Mathematics, the 
American Association for the Advancement 
of Science, Phi Beta Kappa, and Sigma Xi. 


Harry J. Gray, Jr. (S’45—-A’46-M’55) was 
born in St. Louis, Mo., on June 24, 1924. 
He received the B.S. degree in electrical en- 
gineering in 1944, the M.S. degree in elec- 
trical engineering in 1947, and the Ph.D. 
degree in 1953, all from the University of 
Pennsylvania, Philadelphia. 

After serving in the U. S. Navy as a 
radio specialist officer, he returned to the 
University of Pennsylvania in 1946. He 
worked on both the EDVAC and MSAC 
computers and then was associated with the 
development of a digital operational flight 
trainer (UDOFT). He contributed to the 
basic system organization of UDOFT, in- 
troduced the stability chart used to es- 
timate the effect of accumulated truncation 
errors in the numerical integration logical 
packages of UDOFT, and to the peripheral 
equipment. 

In 1954, he joined the Remington Rand 
Univac Division of Sperry-Rand Corpora- 
tion, Philadelphia, Pa., where he worked on 
magnetic and transistor circuits. He devel- 
oped the high-speed circuit system of LARC, 
and when he left in 1957 he held the position 
of staff engineer. 

At present he is associate professor of 
electrical engineering in the Moore School of 
the University of Pennsylvania. 

Dr. Gray is a member of the Association 
for Computing Machinery, Sigma Xi, Tau 
Beta Pi, and Eta Kappa Nu. 


\7 
‘* 


Douglas J. Hamilton (A’53) was born in 
Canton, Ohio, on December 6, 1930. He at- 
tended the College of Wooster, Wooster, 
Ohio, and Case Institute of Technology, 


COMPUTERS 


September 


It is not always desirable to use a canonic 
expansion and the substitution rule. When 
all switching elements in a network are in- 
dependent, it is usually easier to write the 


probabilities from inspection of the network ~ 


using series-parallel reduction. When the 
elements are not independent, the canonic 
expansion and substitution rule are simple, 
powerful tools. 


J. N. WARFIELD 
Purdue University 
Lafayette, Ind 


Cleveland, Ohio, receiving the B.S. degree 
from the latter institution in 1953. From 
1953 to 1957 he was employed by Hughes 
Aircraft Co., Culver City, Calif., where he 
was engaged in the development of tran- 
sistor circuitry for airborne digital com- 
puters. He received the M.S. degree from 
the University of California, Los Angeles, 
in 1956, and in 1957 joined the General Elec- 
tric Company Computer Laboratory, Palo 
Alto, Calif., where he was engaged in the 
development of transistor circuitry for mag- 
netic core and magnetic drum memories. 
Since January, 1958, he has been a research 
assistant at Stanford Electronics Labora- 
tories, Stanford University, Stanford, Calif., 
where he is studying for the doctorate. 

Mr. Hamilton is a member of Eta Kappa 
Nu, Tau Beta Pi, and an associate member 
of Sigma Xi. 


“e 


Vincent J. Korkowski was born on May 
17, 1929, in Brandon, Minn. He graduated 
from Northwestern Television and Electron- 
ics Institute, Minneapolis, Minn., in 1954. 

Since 1954, Mr. Korkowski has been with 
the Research Division of Remington Rand 
Univac, engaged primarily in research on 
magnetic storage devices and logical ele- 
ments. 


~~ 


Charles A. Krause (M’54) was born in 
Wahpeton, N. D., on October 18, 1926. He 
received the B.S. degree in electrical engi- 
neering from Iowa State College, Ames, 
in 1949 and the M.S. degree, also in elec- 
trical engineering, from Purdue Univer- 
sity, Lafayette, Ind., in 1950. He entered 
the U.S. Army in 1944 and attended South 
Dakota State College, Brookings, S. D., and 
Rutgers University, New Brunswick, N. J., 
under the A.S.T.P. program. He was later in 
charge of a telephone repeater and short- 
wave radio station of the Alaska Com- 
munication System (Signal Corps). 


1958 


From 1950 to 1953 he was a research 
engineer at Hughes Aircraft, Culver City, 
Calif, where he was engaged in the design, 
evaluation and prototyping of analog fire 
control system. From 1953 to 1955 he was 
a project engineer in the research depart- 
ment of the National Cash Register Com- 
pany, Electronics Division, Hawthorne, 
Calif., where he performed technical coor- 
dination of a research and development pro- 
gram which included the application to digi- 
tal systems of vacuum tubes, transistors, 
magnetic drums and heads, magnetic cores 
and ferro-electrics. He has been with the 
Norden Division, United Aircraft Corpora- 
tion, Gardena, Calif., since 1955 and is cur- 
rently head of the Special Equipment Sec- 
tion. His work there has included circuitry 
for machine tool control and for analog com- 
puting equipment, reliability studies, and 
the development of large-scale data han- 
dling systems. 

Mr. Krause is a member of Tau Beta Pi, 
Phi Kappa Phi, Pi Mu Epsilon, Eta Kappa 
Nu and Sigma Xi. 


+, 


Rodger R. Lowe (SM’58) was born in Los 
Angeles, Calif., on December 26, 1926. He 
received the B.A. degree in theoretical phys- 
ics from the University of California at Los 
Angeles in January, 1950. 

From 1950 to 1953 he performed de- 
sign and project engineering work in digital 
components and systems for the Electronic 
Engineering Company of California, Santa 
Ana, Calif. In 1954 and 1955 he served the 
National Cash Register Company, Elec- 
tronics Division, Hawthorne, Calif., as re- 
search project engineer in advanced digital 
components and techniques. 

From 1956 to the present time, Mr. Lowe 
has carried out special assignments in digital 
components and systems design and appli- 
cations for the Norden Division, United Air- 
craft Corporation, Gardena, Calif. 


* 


Raymond E. Miller (S’54—M’58) was 

born in Bay City, Mich., on October 9, 1928. 
“He received the B.S. degree in mechanical 

engineering in 1950 from the University of 
Wisconsin in Madison. From August, 1950, 
to May, 1951, he was employed by Inter- 
national Business Machines Corporation in 
Endicott and Poughkeepsie, N. Y. 

Until February, 1953, he served as a 
lieutenant with the U. S. Air Force. He then 
entered the University of Illinois Graduate 
College, and held a research assistantship in 
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the Digital Computer Laboratory from 
June, 1953, to September, 1956, and a re- 
search associateship from June to Septem- 
ber, 1957. He received the B.S. degree in 
electrical engineering in 1954, the M.S. de- 
gree in mathematics in 1955, and the Ph.D. 
degree in electrical engineering in 1957. He 
held the RCA Fellowship from September, 
1956, to June, 1957. 

At present, he is a staff engineer in the 
Information Research Department of the 
IBM Research Center, Yorktown Heights, 
NEY. 

Dr. Miller is a member of the Association 
for Computing Machinery, American Soci- 
ety of Mechanical Engineers, Pi Tau Sigma, 
Tau Beta Pi, Phi Kappa Phi, and Sigma Xi. 


Joseph Otterman (S’52—M’56) was born 
on April 12, 1925, in Warsaw, Poland. He 
attended the Hebrew Institute of Technol- 
ogy (Technion), Haifa, Israel, and was grad- 
uated in 1947 with the degree of Ingenieur 
in Electrical Engineering. During 1948-1950 
he served in the Israeli Army, ending his 
military service as a lieutenant in the ca- 
pacity of engineering reconnaissance officer 
of the Command of the South. From 1950 
to 1951 he worked in the Israeli Government 
Scientific Institute on the development of 
instrumentation systems for ballistic meas- 
urements. 

Coming to the University of Michigan, 
Ann Arbor, in 1951, he obtained the M.S.E. 
and Ph.D. degrees in electrical engineering 
in 1952 and 1955, respectively. 

Since 1955, he has been associated with 
the Engineering Research Institute, Uni- 
versity of Michigan, working on rocket in- 
vestigation of the upper atmosphere within 
the International Geophysical Year pro- 
gram, and on problems of simulation by ana- 
log and digital techniques. 

Dr. Otterman is a member of the Ameri- 
can Physical Society and Sigma Xi. 


William M. Overn was born on March 
20, 1924, in St. Paul, Minn. He served asa 
communications and radar officer with the 
U. S. Air Force in the Pacific area during 
World War II. In 1949 he received the 
B.E.E. degree in physics from the Univer- 
sity of Minnesota, Minneapolis, where he 
did graduate work until 1950. From 1951 
through 1952 he engaged in instrumentation 
work on the U. S. Corps of Engineers Under- 
ground Explosion test series as an employee 
of Engineering Research Associates, Inc. 
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Since 1953 Mr. Overn has been with the 
Research Division of Remington Rand 
Univac, engaged primarily in research on 
magnetic and ferroelectric storage devices 
and logical elements. 
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J. E. Robertson (S’48-A’50-—M’55) was 
born in Fairfax, Okla., on March 30, 1924. 
He received the B.S. degree in electrical en- 
gineering from Oklahoma A. and M. Col- 
lege, Stillwater, Okla., in 1947, the M.S. de- 
gree in 1948, and the Ph.D. degree in 1952, 
both from the University of Illinois. 

He has been associated with the Digital 
Computer Laboratory of the University of 
Illinois since 1950, and is now a research 
associate professor of electrical engineering. 

Dr. Robertson is a member of the Asso- 
ciation for Computing Machinery. 
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Jack L. Rosenfeld was born in Pitts- 
burgh, Pa., on June 6, 1935. In 1957 he re- 
ceived the S.B. and S.M. degrees from the 
Massachusetts Institute of Technology, 
Cambridge, Mass., where he is now studying 
for the Sc.D. degree in electrical engineering. 

He was employed by the Bell Telephone 
Laboratories as a cooperative course stu- 
dent, working on microwave tube develop- 
ment and electronic switching. He has also 
been associated with Texas Instruments, In- 
corporated, and Space Technology Labora- 
tories. 

Mr. Rosenfeld is a member of Tau Beta 
Pi, Eta Kappa Nu, and Sigma Xi. 
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Thomas D. Rossing (SM’56) was born 
in Madison, S. D., on March 27, 1929. He 
received the B.A. degree in physics and 
mathematics from Luther College, Decorah, 
Iowa, in 1950. Graduate study at Iowa State 
College, Ames, led to the M.S. degree in 1952 
and the Ph.D. degree in physics in 1954. 
During this period he conducted research on 
ultrasonic dispersion in gases. 

From 1954 to 1957, he was employed 
in the Research Division of Remington 
Rand Univac. His interests were in the 
field of ferromagnetic materials and devices. 
He has done research on devices for non- 
destructive readout and on vacuum-de- 
posited ferromagnetic films. Since 1957, he 
has been an associate professor of physics at 
St. Olaf College, Northfield, Minn. 

Dr. Rossing is a member of the American 
Physical Society and Sigma Pi Sigma. 
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SCIENCE EDUCATION SUBCOMMITTEE NEWSLETTER 


Vol. feiNo. 2 


PURPOSE 


SENEWS is a newsletter addressed to computer 
oriented members of IRE, ACM, and AIEE to help 
them promote interest and knowledge among high 
school age students; it provides a nationwide medium 
of communication pinpointed to this subject; it evolves 
from volunteer efforts of the JCC Science Education 
Subcommittee and relies on its readers for news mate- 
rial. 

Write to: C. W. Farr, Chairman, JCC Science Edu- 
cation Subcommittee, M.I.T. Lincoln Laboratory, 
Lexington 73, Mass.—or to your representative on the 
SENEWS Editorial Board: Richard W. Melville, IRE; 
George E. Forsythe, ACM; G. L. Hollander, AIEE. 


DANGEROUS VOLTAGE 


SENEWS bristles with stories of computer activity 
in high schools; but do not expect your local authorities 
to welcome you with open arms when you knock on the 
door and announce your intention to help them “dis- 
cover” computers. Dr. L. C. Van Atta, a leader in the 
energetic and successful industry-education program 
in southern California, has written us: “A word of warn- 
ing. It has been my experience that school systems tend 
to regard individual companies, specialized professional 
groups (petroleum engineers, aeronautical engineers, 
electrical engineers), and others as pressure groups un- 
less their offers of assistance to the schools are on a very 
broad basis. Naturally the schools and community are 
interested in the total education of the student, rather 
than in biasing the student toward any particular 
specialization.” 

Dr. Van Atta went on to point out that innovations 
(no matter how good) represent a disruption to estab- 
lished curricula. Since Sputnik the well-meaning offers 
of assistance have multiplied. In Los Angeles the super- 
intendent of schools has created a central committee for 
cooperation, and the school system has provided a full 
time Executive Secretary to coordinate the community 
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efforts and “arrive at a broad program in which al 
specialized subjects have their appropriate emphasis.” 

Avoid burning up the circuit; measure the input im- 
pedance of the school system in your community before 
applying your driving voltage. 


HAvE DeEsicn, WILL BUILD—CHAPTER Two 


In the last issue we reported how David Ecklein, 
a high school junior from Cedar Falls, lowa, was building 
a checker playing digital computer at home, using his 
own design, thirty-five hundred surplus tubes, and the 
help of his buddies in Tom Sawyer style. (Ed. note: 
Last minute correction of number of tubes was not in 
time for all publications.) 

David’s project has been interrupted for a while by 
summer employment at IBM’s research laboratory in 
Poughkeepsie, N. Y., where he will learn to program 
and operate one of the major computers. He is looking 
forward to talking with engineers about circuit designs, 
and learning about the well-known checker playing 
program of Dr. A. L. Samuel. Nobody who has met 17- 
year old David will be greatly surprised if he makes 
some significant improvements in the techniques before 
the summer is out. 

David Ecklein’s opportunity to work in a major 
computer laboratory came about through his contact 
with the JCC Science Education Subcommittee. It is, 
in David’s words, “beyond my fondest hopes. Not only 
will it provide experience but also a means to earn funds 
to carry out my project to conclusion.” 

WERNER BUCHHOLZ 


Junior HicH ScHooL COMPUTER PROGRAMMERS 


The thirteen year old in your home may not have 
learned to write the following coded program for the 
IBM 704 computer to solve (for 5 values each of a, b, 
and c) the equation P=(a+6)(a?+ab+6*)(c+1). But 
don’t be discouraged. We know of only one Junior 
High which has provided such instruction. 


a 


ee 
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TABLE I 
——— 
Opera-| Address, Tag, 
i tion Decrement Remarks 
| SBM | LXA | *HERE, 1 Load index register 
CLA | *A DATA+S5, 1 (a+b) 
FAD | *B DATAGS, 1 
STO *TEMP 
LDQ A DATA4S5, 1 a? 
FMP A DATA4S5, 1 
STO TEMP+1 
LDQ A DATA45, 1 ab 
FMP B DATA+S5, 1 
STO TEMP+2 
LDQ B DATA+5, 1 b? 
FMP B DATA+-S, 1 
STO TEMP+3 . 
CLA TEMP+41 (a2+ab +6?) 
FAD TEMP+2 
FAD TEMP+3 
STO TEMP+4 
CLA | *C DATA4S, 1 (e+1) 
FAD | *INFO 
STO TEMP+5 
LDQ TEMP (a+b) (a?+ab+b?) =(PROD.) 
FMP TEMP+4 
STO TEMP+6 
LDQ TEMP+5 (PROD.)(c+1) =(ANS.) 
FMP TEMP+6 
STO *X INFO-+5, 1 
TIX SBM +1, 1, 1 If contents of index register are 
greater than decrement (1), 
| decrease contents by decrement 
| (1) and transfer to SBM-+1. 
| Otherwise proceed to next in- 
struction (HTR) 
HTR SBM Halt and return to starting po- 


sition. 


* Note: 5 is in location 5 HERE; 5 values of “a” start at location 

A DATA, 5 b’s start at B DATA, 5 c’s start at C DATA: number 1 

% in location INFO; TEMP isa temporary storage register; X INFO 
storage location of final result. 


Harley Tillitt of the U. S. Naval Ordnance Test 
Station, China Lake, Calif., conducted an educational 
experiment with 24 selected eighth grade students at 
Burroughs Junior High School in November, 1957. 
The experiment involved lectures and demonstrations 
during 10 forty-minute instruction periods. Computer 
programmers will recognize the conventions developed 
by SHARE, an informal programming organization 
of IBM 704 users. 

“The success of the students in these experiments 
shows that it is feasible to introduce computer pro- 
gramming instruction below the college level,” con- 
cluded Tillitt in the report submitted to SENEWS. 
The problem coded here (Table I) is one of 14 problems 
assigned to eighth grade students during the experiment. 
The complete report came to SENEWS via George 
Forsythe. (Ed. note: If this report is published in more 
detail, SENEWS will announce it in a later issue.) 


Senews 


255 


SADSAC TI 


Samson Additive Digital Sequential Automatic Com- 
puter, SADSAC II, carried its designer-builder, Peter 
Samson, to first place in the Massachusetts State Science 
Fair, and fourth place in the 1958 National Science 
Fair. Peter lives in Lowell, Mass., and will enter M.I.T. 
this autumn with a four year National Merit Award 
scholarship. He visited Lincoln Laboratory recently on 
invitation from the Science Education Subcommittee— 
and when he left it was we who were inspired. 

First we reviewed specifications. SADSAC II is a 
relay computer, designed to add, subtract, and comple- 
ment four bit binary numbers; input is from two punched 
paper tape readers, and programming provides for 
jumping from one tape to the other; output is a solenoid 
actuated typewriter (actually a home rigged 1898 
Oliver). 

But then we explored motivation, and the “engineer- 
ing economics” of student computers. Peter’s early 
science interest found roots in his father’s electrical 
(hi-fi, etc.) equipment. High school algebra whetted his 
mathematics appetite. Attempting to build a machine 
to play ticktacktoe led him to design his own arithmetic 
circuitry, and the home grown SADSAC I computer was 
the inevitable result. After winning local honors with 
SADSAC I he went on to national fame with its succes- 
sor. 

Peter did not have help from a computer expert ad- 
visor; like most student builders of computers he dug the 
know-how out for himself, “and used common sense,” 
he hastens to add. He regarded pertinent literature as 
scarce but good; the hardware situation he found pretty 
sad. We asked why he used relays instead of tubes or 
transistors. He pointed out that relays are cheap and 
are ac energized at a few standardized voltages; tubes 
and transistors involve not only more money, but more 
complex circuitry and power supplies. He used dis- 
carded pinball machine components, and reliability 
was a nightmare. 

“Building SADSAC III this summer, Peter?” 

“No. I want to get in some study that I can’t find 
time for during the school year.” 


Book AND FILM 


“Local action kit,” prepared by the President’s Com- 
mittee on Scientists and Engineers, Washington 25, D.C. 

This kit is designed to: 1) tell how to bring together 
community organizations with the same basic objective 
into a single program, 2) promote self-evaluation and 
the recording of successful techniques and experiences 
useful to other communities throughout the nation, and 
3) establish a pattern of nationwide effort without 
sacrificing the incentives and benefits of creative local 
action. 

Here is what the kit contains: 1) “Guidebook for 
local action,” suggesting general techniques adaptable 
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TO OPERATE: 


MOVE HANDLE 
RIGHT @ LEFT 


Fig. 1—PAPAC-00, a 2-register, 1-bit, fixed-instruction binary digital computer (Rollin P. Mayer, July 14, 1958). 


to the needs and resources of your own community, 2) 
examples of tested projects in other communities, relat- 
ing how they began, how they are being carried out, 
their scope, financing, and results, 3) reference ma- 
terials such as information on scholarships, improve- 
ment of science curricula, youth activities, etc., and 
4) a bibliography of selected materials, visual aids, 


organizations, and publications helpful in your com- 
munity program. 

This kit was prepared for use in local programs for 
the improvement of science and mathematics education 
(not limited to computer education). It is available 
without charge to groups, not (for budgetary reasons) 
to individuals. 
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SWITCH 
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NOTE: EDGE VIEW SHOWING 
ARRANGEMENT OF PINS: 


BASE 
(SEE FIG.1) one PIN HERE. 


a CLUE HERE. 


fi LEAVE SPACE FOR 
FREE MOTION. 


HEAD OF PIN 
Fig. 2—Parts for PAPAC-00 (Rollin P. Mayer, July 14, 1958). 


Interested individuals can get free copies of “Local in Figs. 1 and 2. This computer was developed from the 
tion,” the monthly news bulletin. (Ed. note: Don’t model demonstrated in the Concord, Mass., High 
> misled by the price; this is good stuff.) School lectures on computers reported in the first 

issue of SENEWS:.! 

PAPAC-00, A Do-IT-YoursELF PAPER COMPUTER From the discussion below, the computer expert will 
In less than an hour you can build the simplified 


gital computer shown in Fig. 1, using only a pair of 1 SENEWS, IRE Trans. ON ELECTRONIC COMPUTERS, vol. EC-7» 
issors, three dozen common pins, and the parts shown pp. 186-187; June, 1958. 


258 


recognize that “PAPAC double zero” contains most of 
the units of a large-scale computer, but in simplified 
form. The control unit includes a counter and a system 
for controlling the parts of the computer according to 
the instruction being performed (in this model a simple 
fixed instruction is used; a large computer can draw 
from several instructions obtained from storage). The 
storage unit includes registers, bus, and selection switch; 
register contents are changed by hand rather than by 
the computer. The arithmetic unit can add. Input and 
output units have been eliminated by allowing the opera- 
tor to deal with the insides of the computer directly 
rather than by way of complicated equipment. Pro- 
prietary rights are held by the author. 

In operation, PAPAC-00 follows the same fixed in- 
struction over and over again. This instruction is: 
“Read the number out of the currently-selected storage 
register and add it to the adder, then get ready to use 
the next storage register for the next time.” The 
“counter” keeps track of which storage register to use 
next; since there are only two registers, numbered “0” 
and “1,” the counter alternates between them. The 
“switch” is controlled by the counter and allows only 
the selected register to be operated. Each “storage 
register” contains only a single binary “cell”; when the 
register is operated, the cell is forced against the “bus” 
if the cell is set to “1.” If a “1% has been read out in this 
way, the bus actuates the “adder,” preparing it to add 
the “1.” If the cell is set to “0,” the bus and adder are 
not operated, and “0” is added to the adder. Binary 
sums are as follows: 0+0=0, 0+1=1, 1+0=1, 
1+1=10. The adder forms these sums correctly except 
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that in the last case it forms a sum of “OQ” because it 
can handle only one digit The “control,” pushed back 
and forth by hand, performs this fixed instruction by 
operating the counter and switch, and by returning the 
bus to its “O” position (if it had read out a “1”) causing 
the sum to be formed in the adder. 

To assemble PAPAC-00, Fig. 1 should be used as the 
base, and the shapes of Fig. 2 should be fitted over it by 
following these steps: 

1) Punch a pinhole exactly through the intersection 
of each cross (+) in Figs. 1 and 2 (but not the dots in 
Fig 1). 

2) Cut out exactly on the lines, the parts in Fig. 2, 
in any order They are marked with a letter in a square 
box, from [A] to [Vv], and the next steps will be easier 
if you place each piece on the table in alphabetic order 
as you cut it out. 

3) Place a pin up through each hole with a circled 
number (from to N 

4) Taking each part of Fig. 2 in alphabetic order, 
place its wncircled number holes down over the cor- 
respondingly numbered pins. 

5) In first operating the computer you may find that 
some parts jam because the upper piece is down too 
far on the pins: pry such pieces up a little to provide 
space for free motion. 

6) The construction can be refined by cutting the 
pins and gluing the uppermost part to the remaining 
length. Caution: 

1) Don’t cut the stop pins too short. 

2) Glue only one moving part to the same pin. 

ROLLIN P. MAYER 
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PGEC News 


1959 WESTERN JOINT COMPUTER 
CONFERENCE 


Papers are being solicited for the 1959 
Western Joint Computer Conference, to be 
held at the Fairmont Hotel, San Francisco, 
Calif., on March 3-5, 1959. The theme of 
this conference will be “New Horizons with 
‘Computer Technology.” 

_ In keeping with the theme there is par- 
ticular emphasis on factual papers dealing 
with the newer applications of computer 
techniques, such as Information Retrieval, 
Operation Control, Pattern Analysis, De- 
cision Making, Computer Communications, 
Learning Concepts, and so forth, as well as 
on papers dealing with advances in com- 
‘puter component and systems design. 

It is also hoped that there will be two 
sessions of a speculative nature: A “Blue 
Sky Session” and “Philosophy and Re- 
sponsibility of Computers in Society.” 
Papers intended for the “Blue Sky Session” 
should deal with the extension of computer 
technology into areas not considered feasi- 
ble at present. They should indicate the ad- 
vantages of such extension, and also the area 
of research necessary to bring this applica- 
tion into the feasible range. Papers for the 
session on “Philosophy and Responsibility 
of Computers in Society” should deal with 
the philosophic and/or social implications of 
the widespread application of automatic 
computer techniques. The papers for these 
sessions should be of the type to invite seri- 
ous discussion. These two sessions will be 
definitely scheduled only if enough suitable 
papers are received. 
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Papers for the 1959 WJCC should be 
prepared based on a 20-minute delivery 
time. Selection of papers for presentation 
will be made from the complete text of the 
paper. There are no format requirements for 
these submission drafts. Three copies of the 
proposed paper should be submitted to 
Murray L. Lesser of the Technical Program 
Committee by October 1, 1958, and ad- 
dressed to Mr. Lesser at IBM Research Lab- 
oratory, San Jose, Calif. After review, final 
selection of papers will be made and the 
authors will be notified by December 1. Sub- 
mission of final texts of the selected papers, 
in the form required by the Publications 
Committee, should be made by February 1, 
1959. 

The 1959 WJCC makes the classical 
promise to be an _ exciting conference. 
Whether or not you are in a position to pre- 
sent a paper, you will find it well worthwhile 
to reserve March 3-5, 1959, to attend. 


1958 WESTERN JOINT COMPUTER 
CONFERENCE 


On May 6 through 8 the very successful 
1958 Western Joint Computer Conference, 
“Contrasts in Computers,” was held at the 
Ambassador Hotel, Los Angeles, Calif. Re- 
garded by the Conference Committee as the 
“150 per cent Conference” because every 
event—technical as well as social—was 
about half again as large as predicted, the 
Conference registered a total of over 1800 
people. It was the largest of the Western 
conferences held so far. An outstanding pro- 
gram was presented by the Technical Pro- 
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gram Committee. In the keynote session, 
which was devoted to exploring the possible 
social problems which may arise from wide- 
spread automation, a number of very in- 
triguing and significantly new computer 
applications in the social and behavioral 
sciences were suggested. Six invited panel 
discussions presented and explored several 
of the controversial questions of the present 
computer art, and six sessions of contributed 
papers described some of the recent ad- 
vances in computer hardware and computer 
applications. A “bonus feature” of the tech- 
nical program was the added technical ses- 
sion devoted to a paper describing the re- 
cently announced French Gamma-60 com- 
puter. 

Fifty exhibitors demonstrated recent 
components, equipment, techniques, and 
company progress in the 80 booths of the ex- 
hibit area. The social events of the Confer- 
ence—the women’s activities, the cocktail 
party and the luncheon—were all over-sub- 
scribed, and contributed to the over-all 
conference success and atmosphere. Follow- 
ing the Conference a one-day Symposium 
on “Small Automatic Computers and Input 
/Output Equipment” was sponsored by the 
Los Angeles Chapter of the ACM, and at- 
tended by 350 people. Eight papers from the 
manufacturers described their latest equip- 
ment. 

Proceedings of the Conference will be 
published in the late fall of 1958 and will be 
available to nonregistrants from the three 
co-sponsoring societies: the IRE, the ACM, 
and the AIEE. A Proceedings of the Sym- 
posium will be published at a later date by 
the ACM. 


by: 


Nore 


Note: Only considered are original papers not published or presented prior to 
the 1959 IRE National Convention; any necessary military or company clear- 


Call for Papers 


1959 IRE NATIONAL CONVENTION 
March 23-26, 1959 


Waldorf-Astoria Hotel and New York Coliseum, New York, N. Y. 


Prospective authors are requested to submit all of the following information 


October 24, 1958 


. 100-word abstract in triplicate, title of paper, name and address 
500-word summary in triplicate, title of paper, name and address 

3. An indication of the technical field in which the paper falls (e.g., ELEC- 
TRONIC COMPUTERS). 


ance of paper must be granted prior to submittal. 


Address all material to: Dr. George L. Haller, Chairman 


1959 Technical Program Committee 
The Institute of Radio Engineers, Inc. 
1 East 79 Street, New York 21, N. Y. 
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INFORMATION FOR AUTHORS 


The PGEC TrRansactTIONs is published quarterly and will bear date- 
lines of March, June, September, and December. Abstracts of papers 
appearing in the TRANSACTIONS will appear also in [RE PROCEEDINGS. 
The PGEC publication schedule requires about one month for review 
and correction of all accepted manuscripts. The professional IRE 
Editorial Staff requires an additional two months’ production time from 
receipt of manuscripts to completion of the printed journal. 


MANUSCRIPTS: Three copies of the manuscript should be sub- 
mitted. They should be typewritten (original and two carbon copies), 
and double spaced on only one side of each sheet. References should 
appear as footnotes, numbered consecutively, and include in the fol- 
lowing order the author’s name (including initials), title of reference 
work, journal name, volume, initial and final page numbers, and date 
of publication. Footnotes should be listed on a separate sheet and not 
inserted in text. Each paper must be accompanied by two copies of a 
summary not more than 200 words in length. Reviewing normally will 
require about four weeks from receipt of manuscript. 


ILLUSTRATIONS: Only original illustrations should be submitted; 
they will be returned if desired. Photostatic copies of originals are not 
acceptable, except where they are exceptionally clear, with sharp black 
and white contrasts. All line drawings (graphs, charts, diagrams, etc.) 
should be prepared on drafting cloth or white drawing paper in India ink. 
It is preferable that only the coordinate lines show in graphs. All 
lettering must be large enough to be legible when reduced 50 to 75 per 
cent in size. Photographs should be glossy prints. Figure numbers 
should be indicated on the back of each illustration. Figure numbers 
and captions should be listed on a separate sheet accompanying manu- 
script. All drawings, photographs, and other manuscript material should 
be not larger than 83 by 11 inches for ease in handling. 


Please submit all manuscripts to 

J.P, Nash, PGEC Editor: 
Missile Systems Division, Lockheed Aircraft Corp., 
3251 Hanover Street, Palo Alto, California. 


